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ABSTRACT 

Strategic  nuclear  targeting  studies  generally  include  more  target 
installations  than  there  are  weapons.  Hence,  a  weapon  is  not  assigned 
to  an  installation,  but  rather,  to  a  Desired  Ground  Zero  (DGZ).  The 
objective  of  this  study  was  to  investigate  optimal  DGZs  for  multiple 
nuclear  weapons  against  installations  in  a  target  complex.  To  accomplish 
this,  it  was  necessary  to  develop  the  target  Complex  Expected  Damage 
Function  (CEDF)  maximization  algorithm.  The  algorithm  locates  optimal 
DGZs  by  maximizing  the  CEDF;  the  CEDF  is  a  nonlinear  function  of  2m  vari¬ 
ables,  the  (X^,Y^)  DGZ  coordinates  for  each  of  the  m  weapons. 

The  algorithm  uses  two  CEDF  models  and  two  optimization  techniques. 
These  models  use  DIA  Physical  Vulnerability  System  probability  of  damage 
models.  The  CEP- Included  model  includes  each  weapon's  CEP;  the  simpler 
CEP-Excluded  model  assumes  each  weapon's  CEP  equals  0.  An  analytical 
expression  for  the  gradient  of  the  CEP-Excluded  model  was  calculated; 
the  algorithm  maximizes  this  CEDF  using  a  conjugate  gradient  with  restarts 
search  technique.  The  algorithm  maximizes  the  CEP- Included  CEDF  using  a 
direct  search  technique,  Powell's  method  of  conjugate  directions. 

This  investigation  characterized  three  factors  that  affect  the 
optimal  DGZ  locations  for  multiple  nuclear  weapons  in  a  target  complex. 

The  first  factor  wAs  gradient  symmetry;  this  symmetry  resulted  from  either 
a  geographically  symmetric  target  complex  or  collocated  weapons.  The 
second  factor  was  weapon  CEP.  Maximization  of  the  two  CEDF  models 
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produced  slightly  different  optimal  DGZs;  this  difference  depended  on  a 
weapon's  CEP  and  the  CEDF  model.  The  third  factor  was  the  initial  DGZ 
location  prior  to  CEDF  maximization.  The  algorithm  located  different 
CEDF  local  maximums  depending  on  the  initial  DGZ  condition.  However, 
the  investigation  revealed  that  the  most  successful  initial  XZ  condi¬ 
tion  is  to  use  the  coordinates  of  the  highest  valued  installations  as 
the  initial  DGZ  coordinates. 
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AN  INVESTIGATION  OF  OPTIMAL  AIMPOINTS 
FOR  MULTIPLE  NUCLEAR  WEAPONS  AGAINST 
INSTALLATIONS  IN  A  TARGET  COMPLEX 

I.  Introduction 

Effective  U.  S,  targeting  of  an  enemy's  resources  is  an  important 
part  of  U.  S.  military  air  power.  One  of  the  fundamental  objectives  of 
U.  S.  military  forces  is  to  sustain  deterrence  (Ref  8:  para  1-6). 
Deterrence  is  an  enemy's  state  of  mind  brought  about  by  the  existence 
of  U.  S.  military  power  or  the  enemy's  perception  of  U.  S.  resolve  to 
use  that  power.  Strategic  nuclear  targeting,  an  assignment  process, 
is  a  key  element  of  nuclear  deterrence.  The  nuclear  weapons  planner 
must  assign  a  weapon  system  to  a  specific  target.  Targeting  consists 
of  three  interacting  processes:  the  target  intelligence  process,  the 
threat  estimate  process,  and  the  operational  planning  process  (Ref  7:2-2). 
This  study  investigates  an  important  phase  of  the  target  intelligence 
process,  nuclear  weaponeer ing,  and  a  weaponeer ing  problem.  Lee  defines 
waaponeering  as  "the  process  that  determines  the  physical  vulnerabilities 
of  targets,  the  optimum  weapon  type,  the  number  of  weapons,  and  sometimes 
the  best  system  required  to  achieve  a  desired  level  of  damage  on  a  target 
or  a  target  system"  (Ref  18:122). 

Background 

Weapons  planners  allocate  weapons  to  Desired  Ground  Zeros  (DGZs) 

to  achieve  damage  to  installations  within  a  target  complex.  A  target 
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complex  is  a  geographical  area  that  includes  different  types  and  numbers 
of  target  elements  or  installations.  Nuclear  detonations  within  a 
complex  will  cause  insignificant  damage  to  installations  within  all 
adjacent  complexes  (Ref  20:6).  A  complex  may  contain  one  installation 
or  a  few  hundred  installations.  For  example,  a  50-square-mile  Air  Force 
base  may  be  a  target  complex.  Similarly,  the  runway,  the  maintenance 
facility,  the  parked  aircraft,  and  the  headquarters  command  post  are 
installations  of  this  target  complex.  However,  weapons  planners  do  not 
allocate  weapons  to  each  installation.  Instead,  they  allocate  weapons 
to  DGZs  and  plan  to  damage  more  than  one  installation  with  one  weapon. 

A  DGZ  is  a  point  on  the  surface  of  the  earth  at  or  vertically  below  the 
center  of  a  planned  weapon  explosion  (Ref  7:5-6).  It;  this  study,  a  DGZ 
refers  to  a  nuclear  weapon  detonation  at  a  specific  geographical 
location.  A  DGZ  may  be  located  directly  on  an  installation;  or,  if 
one  weapon  will  sufficiently  damage  two  or  more  installations,  then  the 
DGZ  may  be  located  between  the  installations. 

Weapons  analysts  use  the  concept  of  a  lethal  aimpoint  region  (LAIR) 
to  locate  DGZs  within  a  target  complex  (Ref  20:10;  21:4;  and  25:2-6). 

Ihe  LAIR  is  a  circular  area  whose  center  is  the  target  installation. 

It  represents  a  geographical  region  within  which  a  weapon  can  detonate 
and  achieve  at  least  a  minimum  probability  of  damage  (Pd)  to  a  target. 

Pd  is  the  probability  that  a  desired  level  of  damage  (severe,  moderate, 
light)  will  occur  to  a  target  (Ref  7:5-6).  "The  general  definitions 
of  the  three  damage  levels  are:  (1)  severe  damage  —  a  level  which 
requires  essentially  complete  reconstruction  or  replacement  of  one  or 
more  critical  major  elements  of  the  target,  plus  reconstruction,  repair, 
or  replacement  of  associated  structures  or  equipment.  Severe  damage 
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precludes  utilization  of  the  target  for  any  purpose,  (2)  moderate 
damage  —  a  level  which  requires  major  repairs  to  one  or  more  critical 
major  elements  of  the  target,  plus  major  reconstruction,  repair,  or 
replacement  of  associated  structures  or  equipment.  Moderate  damage 
precludes  effective  utilization  of  the  target  for  its  intended  purpose, 

(3)  light  damage  —  a  level  which  does  not  significantly  impair  the 
target  function,  but  requires  some  repairs  to  restore  the  target  to 
complete  usefulness"  (Ref  19:1-7). 

The  radius  of  the  LAIR  depends  on  specific  weapon  system  and  target 
parameters.  The  accuracy  of  the  missile  or  aircraft  system  that  delivers 
the  weapon  to  the  DGZ  affects  the  LAIR.  Also,  the  yield  of  the  nuclear 
weapon  affects  the  LAIR.  Yield  is  a  numeric  value  measured  in 
kilotons  (kt)  and  is  a  relative  indicator  of  the  explosive  energy  the 
weapon  releases  when  it  detonates.  This  explosive  energy  causes  damage 
to  installations.  A  nuclear  weapon  distributes  its  damage  energy  in 
several  ways  through  damage  mechanisms  or  weapon  effects.  For  ground 
targets,  the  most  prominent  mechanism  is  the  blast  effect.  The  primary 
elements  of  blast  are  overpressure  and  dynamic  pressure.  Overpressure 
creates  a  force  that  crushes  an  installation;  dynamic  pressure  creates 
a  force  from  the  resulting  high  wind  velocity  (Ref  11:80-82).  But 
thermal  effects,  cratering,  and  impulse  are  other  nuclear  weapon  effects 
that  may  contribute  to  target  damage.  The  occurrence  and  intensity  of 
these  weapon  effects  vary  for  different  weapon  yields. 

The  LAIR  also  depends  upon  target  characteristics,  specifically, 
the  vulnerability  of  the  target  to  blast  effects.  'Hie  Defense  Intelligence 
Agency  (DIA)  uses  a  Physical  Vulnerability  coding  system  to  quantify  a 
target's  susceptibility  to  blast  damage.  Each  installation  is 
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characterized  by  a  three-part  Vulnerability  Number  (VN).  The  first 
part  consists  of  a  two-digit  integer  reflecting  the  target's  relative 
hardness  in  terms  of  a  20-kt  weapon  and  a  specified  damage  level 
(severe,  moderate,  light).  The  second  part  is  a  letter  indicating 
whether  the  target  is  predominantly  sensitive  to  either  overpressure 
(L,M,N,0,P)  or  dynamic  pressure  (Q,R,S,T,U).  The  third  part  is  a 
K  factor.  This  factor  adjusts  the  target's  relative  hardness  for  weapon 
yields  other  than  20-kt  (Ref  6:34  and  19:1-7). 

In  this  paper,  four  factors  characterize  a  nuclear  weapon  —  yield, 
accuracy,  height  of  burst,  and  probability  of  arrival  (Pa).  Circular 
error  probable  (CEP)  is  a  numeric  value  measured  in  units  of  length 
that  represents  a  weapon's  delivery  accuracy.  A  500-foot  CEP  indicates 
a  weapon  has  a  50X  chance  of  being  delivered  within  500  feet  of  the 
target.  Similarly,  height  of  burst  is  the  weapon's  distance  above  the 
ground  when  the  weapon  detonates.  Pa  is  the  probability  that  a  delivery 
vehicle  (bomber,  missile)  and  its  weapon  arrive  at  the  target  and  the 
weapon  detonates  as  planned.  Pa  depends  upon  the  delivery  vehicle's 
pre- launch  survivability  (FLS),  weapon  and  weapon  system  reliability  (WSR), 
and  probability  to  penetrate  (PTP).  Each  of  these  factors  is  a 
probability  (Raf  7:5-7). 

The  weapons  analyst  plans  to  damage  installations  within  a  target 
complex  by  assigning  weapons  to  a  prioritized  list  of  DGZs.  In  addition 
to  Pd  and  Pa,  which  are  multiplied  together  to  calculate  an  installation's 
Damage  Expectancy  (DE),  the  value  of  each  installation  is  needed  to 
develop  the  prioritized  list.  The  value  of  an  installation  is  a  number 
that  represents  the  value  of  the  installation  relative  to  all  other 
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installations.  Most  value  rystems  cardinally  order  targets  over  a 
range  from  the  most  valued  target  (highest  value  number)  to  the  least 
valued  target  (Ref  7:6-19).  The  total  complex  expected  target  value 
damage  is  the  sun  of  each  installation's  value  multiplied  by  the 
installation's  cumulative  DE. 

There  is  a  shortfall  in  the  nuclear  weaponeering  process.  The 
prioritized  target  list  generally  has  more  DGZs  than  there  are  weapons 
available  to  assign  to  the  DGZs.  Hie  weapons  analyst  must  determine 
not  only  the  best  weapon-DGZ  combination  to  achieve  the  desired  attack 
objectives,  but  also  alternative  combinations  (Ref  7:5-6). 

The  method  that  strategic  nuclear  weapon  targeting  models  use  to 
address  this  problem  depends  on  the  specific  objective  of  the  model. 

One  objective  is  to  minimize  the  number  of  weapons  required  to  achieve 
at  least  a  minimum  acceptable  probability  of  damage  to  all  installations 
in  the  complex.  This  method  determines  the  minimum  number  of  weapons 
when  installation  Pds  are  prespecified.  A  different  objective  is  to 
achieve  the  maximum  total  expected  target  value  damage  for  the  complex. 
This  method  determines  the  Fd  to  each  installation  when  the  number  of 
weapons  available  is  prespecified. 

As  an  example,  suppose  a  preliminary  target  analysis  indicates  five 
DGZs  are  necessary  to  achieve  a  minimum  acceptable  Pd  for  each  instal¬ 
lation  in  a  complex.  However,  after  allocating  weapons  to  the  entire 
prioritized  DGZ  list  (all  complexes),  only  four  i/eapons  are  actually 
available  to  this  complex.  Should  the  four  weapons  be  targeted  against 
the  four  highest  expected  target  value  DGZs  or  should  an  attempt  be 
made  to  locate  four  new  DGZs,  perhaps  unrelated  to  the  five  potential 
DGZs?  The  former  choice  will  achieve  a  minimum  acceptable  Pd  on  some, 
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but  not  all,  of  the  installations  in  the  complex.  The  installations 
that  would  have  been  damaged  by  a  weapon  allocation  to  the  unassigned 
DGZ  probably  will  receive  insignificant  damage.  Conversely,  the  latter 
choice  may  increase  the  total  expected  target  value  damage  to  the 
complex  with  either  no  decrease  or  a  minimal  decrease  in  the  mininnjm 
acceptable  Pd  for  each  Installation. 

According  to  a  15  September  1983  literature  review,  AF  Studies  and 
Analysis,  Command  and  Control  Technical  Canter  (CCIC),  and  the  Single 
Integrated  Operational  Plan  (S10P)  Simulation  Branch,  Joint  Strategic 
Target  Planning  Staff  (JSTPS)  use  different  models  for  DGZ  optimization 
studies  (Ref  4;  13)  20}  22;  and  28).  Each  of  these  mathematical  models 
has  a  limitation.  Initially,  the  algorithms  generate  a  DGZ  list  for  a 
complex  using  an  unlimited  supply  of  weapons.  Then  the  algorithms 
assign  weapons  either  to  the  minimum  number  of  DGZa  required  to  achieve 
an  acceptable  level  of  damage  on  all  targets  or,  when  the  numbers  of 
weapons  are  constrained,  to  the  DGZs  that  achieve  the  best  total  expected 
target  value  damage  for  the  preplanned  DGZs.  The  second  situation, 
limited  weapon  supply,  is  more  realistic.  However,  the  development  of 
new  DGZs  in  the  constrained  weapons  case  to  maximize  total  expected 
target  value  damage  is  not  attempted.  In  some  algorithms,  DGZs  are 
relocated,  but  relocation  is  sequential.  One  DGZ  is  moved  until  its 
contribution  to  the  total  expected  target  value  damage  is  maximized, 
then  that  DGZ  is  assigned,  and  a  second  DGZ  is  sequentially  moved. 

Problem  Statement 

After  a  weapon  allocation  for  all  target  complexes  is  completed, 

not  all  complexes  may  be  allocated  enough  weapons  to  achieve  an 
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acceptable  Pd  for  all  installations. 

An  algorithm  is  needed  that  will  optimally  locate  DGZs  in  a  target 
complex  for  a  fixed  number  of  weapons ,  while  maximizing  the  total 
expected  target  value  damage  to  installations  within  the  complex. 

DGZ  Models 

Multiweapon  Optimizer  for  Strategic  Targets  (MOST),  Seiler,  and 
NUCWAVE  are  mathematical  models  that  Air  Force  agencies  use  for  stra¬ 
tegic  targeting  studies  (Ref  20;  21;  and  25).  The  models  locate  DGZs 
within  a  target  complex.  Each  of  these  models  depends  on  the  LAIR 
concept  and  uses  either  partial  enumeration,  or  linear  programming 
allocation,  or  sequential  allocation  to  determine  a  set  of  DGZs  for  a 
complex. 

MOST  determines  a  DGZ  list  in  two  phases.  Each  phase  satisfies  an 
associated  criteria.  These  phases  allow  MOST  to  achieve  its  objective  — 
determining  the  fewest  number  of  weapons  (DGZs)  required  to  achieve  at 
least  a  minimun  acceptable  Pd  for  each  installation  in  a  complex 
(Ref  21).  There  are  several  steps  in  the  first  phase.  Initially, 

MOST  generates  a  LAIR  for  each  installation.  These  LAIRs  satisfy  the 
criteria  to  achieve  a  minimum  acceptable  Pd  on  all  installations.  Next, 
MOST  compiles  subsets  of  DGZs  through  a  partial  enumeration  process; 
each  subset  contains  a  list  of  LAIR  intersections.  For  one  subset,  all 
installations  in  the  complex  must  be  included  in  at  least  one  intersection. 
Then  MOST  selects  the  subset  that  contains  the  fewest  number  of  LAIR 
intersections;  if  several  equivalent  subsets  require  the  fewest  number 
of  aimpoints,  then  the  subset  with  the  highest  total  expected  target 
value  damage  is  selected.  As  an  example,  consider  the  target  cotriplex 
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Figure  1.  LAIR  Intersections 


in  Figure  1.  Ihe  algorithm  would  select  the  DGZ  subset  that  contains 
the  LAIR  intersection  of  installations  B,  C,  and  D  as  one  DGZ  and  the 
LAIR  intersection  of  installations  A  and  E  as  the  second  DGZ  of  this 
target  complex  (the  shaded  regions  in  Figure  1). 

In  the  second  phase,  MOST  adjusts  the  final  DGZ  locations  within 
the  LAIR  intersection  regions  using  a  weighted  installation  value  system. 
This  process  maximizes  the  total  expected  target  value  damage  for  all 
DGZa.  If  installation  D  was  more  valuable  than  B  and  C  in  Figure  1, 
then  the  actual  DGZ  would  be  moved  proportionately  closer  to  instal¬ 
lation  D.  These  adjustments  to  final  DGZs  are  accomplished  sequentially. 
First,  the  DGZ  associated  with  the  greatest  number  of  LAIR  intersections 
would  be  maximized  (the  DGZ  associated  with  the  intersection  of  target 
LAIRs  B,  C,  and  D).  Then  the  DGZ  associated  with  the  second  greatest 
number  of  LAIR  intersections,  etc.  MOST  was  designed  to  find  the 
mininun  number  of  DGZs  for  Poseidon  ra-entry  vehicles,  irrespective 
of  weapon  supply  constraints.  If  there  are  not  enough  weapons  to 
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allocate  to  the  complex,  then  the  least  valued  DGZs  would  not  be 

assigned  weapons. 

The  objectives  and  purpose  of  the  Seiler  model  are  similar  to 
MOST.  Seiler  was  designed  to  study  the  assignment  of  nuclear  weapon 
missile  systems  (ICBM  and  SLBM)  to  installations  within  many  target 
complexes  (Ref  20).  Seiler  also  uses  two  phases  to  assign  weapons  to 
a  prioritized  list  of  DGZs.  In  the  first  phase,  generation  of  aim- 
points,  Seiler  creates  DGZs  using  the  LAIR  concept  and  a  tiered  DGZ 
system.  The  primary  tier  consists  of  the  minimum  number  of  DGZs  that 
are  required  to  achieve  a  minimum  acceptable  Pd  to  all  installations 
when  only  the  largest  yield  weapon  is  considered.  For  each  subset  of 
installations  contained  in  a  primary  tier  DGZ,  supplementary  DGZs  are 
created  for  the  next  largest  yield  weapon.  Supplementary  DGZs  are 
always  subsets  of  a  primary  tier  DGZ  or  a  higher-tiered  supplementary 
DG2.  Each  DGZ,  supplementary  or  primary,  achieves  a  minimum  acceptable 
Pd  on  a  subset  of  the  installations  in  a  target  complex  and  has  an  asso¬ 
ciated  DGZ  value.  This  value  depends  on  the  cumulative  total  expected 
target  value  damage  of  the  associated  installations. 

In  the  second  phase,  Seiler  uses  a  linear  programming  (IP)  algo¬ 
rithm  to  determine  an  optimal  (or  near  optimal)  assignment  of  weapons. 
The  LP  objective  is  to  maximize  the  total  complex  expected  target  value 
damage.  Seiler  accomplishes  this  assignment  using  the  primary  and 
supplementary  tiered  DGZs,  missile  delivery  vehicle  range  capabilities, 
and  constraints  on  the  number  of  primary  and  supplementary  tier  weapons 
available.  If  there  are  not  enough  weapons  available  to  allocate  to 
the  installations  in  the  complex,  then  lower  value  DGZs  (and  hence 


installations)  remain  untargeted  just  as  in  MOST. 

NUCWAVE  determines  the  number  and  the  location  of  DGZs  using  a 
different  approach  (Ref  25).  It  is  a  one-sided  nuclear  weapons  allo¬ 
cation  war  gaming  model.  The  user  can  select  one  of  two  strategies  — 

(1)  allocate  a  limited  number  of  weapons  to  DGZs  in  order  to  maximize 
the  total  expected  target  value  damage  to  all  target  complexes  con¬ 
sidered,  or  (2)  determine  the  minimum  number  of  DGZs  required  to  achieve 
a  minimum  acceptable  Pd  to  installations  within  all  target  complexes. 
NUCWAVE  generates  DGZs  using  the  LAIR  concept,  similar  to  MOST  and 
Seiler,  for  allocation  strategy  2.  Allocation  strategy  1  is  accomplished 
using  a  sequential  algorithm  and  will  be  discussed  later. 

The  NUCWAVE  algorithm  consists  of  three  phases,  irrespective  of 
the  allocation  strategy  chosen.  The  first  phase,  potential  allocation, 
uses  an  unlimited  supply  of  weapons  to  maximize  the  damage  attained  by 
each  weapon  until  a  sufficient  number  of  potential  DGZs  are  located  to 
satisfy  the  strategy  objective.  In  the  second  phase,  an  LP  weapon 
selection  program  uses  these  potential  DGZs  and  weapon  supply  con¬ 
straints  to  select  the  number  and  the  type  of  weapons  to  be  assigned 
to  each  complex.  In  the  final  phase,  real  allocation,  the  specific 
number  and  types  of  weapons  are  "optimally  allocated"  against  the  instal¬ 
lations  in  each  complex  (Ref  25:2-8).  If  allocation  strategy  2  has  been 
selected  and  the  weapon  selection  program  allocated  fewer  than  the 
required  number  of  weapons  to  a  complex,  then  the  lower  valued  DGZs  and 
their  associated  installations  would  not  be  targeted.  No  relocation  of 
the  DGZs  is  attempted  with  strategy  2. 

When  allocation  strategy  1,  maximize  total  expected  target  value 
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damage,  is  selected,  NUCWAVE  determines  DGZs  sequentially  in  both  the 
potential  and  real  allocation  phases.  NUCWAVE  starts  by  locating  the 
first  DGZ  at  the  highest-valued  installation  in  the  complex.  The  algo¬ 
rithm  then  moves  the  DGZ  to  a  location  that  maximizes  the  total  complex 
expected  target  value  damage.  The  algorithm  may  move  the  DGZ  closer  to 
several  installations,  thus  increasing  the  Pd  and  expected  target  value 
damage  for  these  installations.  Similarly,  the  algorithm  may  move  the 
DGZ  farther  away  from  other  installations,  thus  decreasing  the  Pd  and 
expected  target  value  damage  for  these  installations.  When  the  optimal 
location  is  determined,  two  steps  occur.  First,  a  weapon-DGZ  location, 
having  been  determined,  is  stored.  Next,  the  surviving  value  of  all 
installations  is  calculated  by  multiplying  the  previous  value  of  the 
installation  by  the  Pd  of  the  installation  from  the  current  weapon-DGZ 
combination.  NUCWAVE  then  selects  the  installation  with  highest  sur¬ 
viving  value  and  the  entire  process  is  repeated.  DGZs  are  sequentially 
determined  in  this  manner  until  the  specified  stopping  condition  is 
reached.  In  the  potential  allocation  phase,  the  stopping  condition  is 
that  a  user-specified  percent  of  the  total  expected  target  value  damage 
has  been  achieved;  in  the  real  phase,  the  condition  is  no  more  remaining 
weapons. 

After  a  weapon  allocation  is  made  for  the  entire  target  list,  only 
a  finite  nunber  of  weapons  may  be  assigned  to  a  target  complex.  Only 
NUCWAVE  attempts  to  locate  new  DGZs,  but  it  uses  a  sequential  optimi¬ 
zation  algorithm.  When  less  than  the  desired  number  of  weapons  are 
allocated  to  a  target  complex,  a  simultaneous  optimal  solution  speci¬ 
fying  the  location  of  the  final  DGZs  should  exist.  In  this  study, 
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optimal  means  the  best  location  of  DGZs  stjh  that  the  total  complex 
expected  target  value  damage  is  maximized. 

Objectives 

the  primary  objective  of  this  study  is  to  investigate  the  optimal 
DGZ  locations  within  a  target  complex.  In  order  to  accomplish  this, 
it  is  necessary  to  develop  an  algorithm.  This  algorithm  will  opti¬ 
mally  locate  the  DGZs  for  fixed  numbers  of  weapons  in  a  target  complex 
by  maximizing  the  expected  target  value  damage  to  all  installations. 

The  algorithm  will  not  be  restricted  to  one  type  of  weapon;  that  is, 
different  weapons  may  be  included  in  the  fixed  number  of  weapons. 

Also,  it  will  be  necessary  to  determine  the  sensitivity  of  the 
algorithm  to  two  factors  —  first,  the  mathematical  technique  used  to 
locate  the  optimal  DGZs;  second,  the  initial  starting  conditions  (lati¬ 
tude  and  longitude  coordinates)  for  the  DGZs. 

The  algorithm  will  consist  of  two  elements.  The  first  element  is 
a  mathematical  model  of  the  total  complex  expected  target  value  damage. 

The  second  element  is  an  optimization  technique  to  determine  the  maximum 
total  complex  expected  value  damage  and  to  locate  the  corresponding 
optimal  DGZs.  The  following  steps  are  an  outline  of  the  algorithm: 

1.  Specify  target  installation  parameters.  These  include 
installation  coordinates,  VN  numbers,  and  values. 

2.  Specify  weapon  parameters.  These  include  yield,  quantity, 

CEP,  Fa,  and  Haight  of  burst. 

3.  Specify  either  the  mathematical  form  or  an  acceptable  approx¬ 
imation  of  the  probability  of  damage  function  for  an  installation. 

4.  Determine  the  mathematical  form  of  the  Installation  Expected 
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Damage  Function  (IEDF).  This  function  represents  the  total 
expected  target  value  damage  to  an  installation  from  all 
weapons. 

5.  Specify  the  form  of  the  Complex  Expected  Damage  Function  (CEDF). 
This  function  is  a  summation  of  all  of  the  IEDFs. 

6.  Select  a  nonlinear  optimization  technique  to  maximize  the  CEDF 
and  to  locate  the  final  coordinates  of  the  DGZs. 

Scope  and  Assumptions 

This  study  will  develop  an  algorithm  subject  to  certain  restrictions 
that  optimally  locates  DGZs  in  a  target  complex.  Secondary  damage  will 
be  assumed  within  the  target  complex;  however,  secondary  damage  from 
weapons  detonated  in  adjacent  complexes  will  not  be  considered.  Also, 
the  algorithm  will  not  consider  target  avoidance  areas. 

Only  military/industrial  installations  that  can  be  modeled  as  point 
targets  will  be  considered.  Also,  since  blast  is  the  primary  damage 
mechanism  for  ground  targets,  other  nuclear  weapon  effects  will  not  be 
considered.  Each  installation's  susceptibility  to  overpressure  and  to 
dynamic  pressure  will  be  specified  with  VN  numbers.  Also,  the  mathe¬ 
matical  model  of  the  probability  of  damage  function  developed  by  the 
Defense  Intelligence  Agency  (DIA)  will  be  used  to  specify  the  instal¬ 
lation  expected  damage  function  (Ref  6). 

The  algorithm  will  consider  weapon  system  delivery  methods  and 
accuracy  since  they  will  affect  the  expected  target  value  damage. 

Delivery  methods  will  be  characterized  by  a  specified  Pa  for  each  weapon. 
However,  feaaible  delivery  constraints  will  not  be  considered,  for 
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example,  Multiple  Independent  Reentry  Vehicle  (MIRV)  footprinting  and  a 
weapon  delivery  system's  range  capability.  A  circular  normal  distri¬ 
bution  will  be  assumed  for  weapon  system  accuracy,  and  a  CEP  will  be 
specified  for  each  weapon.  Pa  and  CEP  are  a  function  of  range  for  some 
weapons,  but  in  this  study  they  are  prespecified  numbers. 

Different  optimization  techniques  and  initial  DGZ  conditions  will 
be  evaluated.  These  evaluations  will  characterize  the  properties,  capa¬ 
bilities,  and  any  limitations  of  the  algorithm. 

Overview 

This  paper  reports  the  methods  and  findings  of  a  study  that  inves¬ 
tigated  the  location  of  optimal  aimpoints  for  multiple  nuclear  weapons 
against  installations  in  a  target  complex.  A  CEDE  maximization  algorithm 
was  developed  to  optimally  locate  DGZs  for  these  weapons  by  maximizing 
the  Complex  Expected  Damage  Function  (CEDF).  The  algorithm  consists  of 
two  elements  --a  mathematical  model  of  the  CEDF  and  an  optimization 
technique.  Chapter  II  presents  the  mathematical  formulation  of  two 
CEDF  models  and  tha  gradient  for  one  of  these  models.  Chapter  III  pre¬ 
sents  an  overview  of  numerical  search  techniques;  it  also  discusses  the 
two  techniques  that  are  used  to  maximize  the  two  CEDF  models.  Chapter  IV 
contains  the  computerization  of  the  algorithm  and  the  verification  and 
validation  processes.  Chapter  V  discusses  the  algorithm's  convergence 
criteria,  and  specific  properties  of  symmetric  target  complexes,  and 
symmetric  CEDF  gradient  elements.  Chapter  VI  is  an  analysis  of  optimal 
DGZs  for  three,  four,  and  seven  installation  target  complexes.  It  also 
discusses  and  summarizes  the  three  conclusions  of  this  study.  Finally, 
Chapter  VII  presents  concluding  remarks  and  recommendations. 
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II.  Mathematical  Formulation  of  the  Complex 


Function  and  Gradient 


The  algorithm  determines  the  optimal  DGZ  locations  within  a  target 
complex  by  maximizing  the  Complex  Expected  Damage  Function  (CEDF). 
Initially,  the  CEDF  is  developed  from  prespecified  weapon  and  instal¬ 
lation  parameters;  then  the  CEDF  is  maximized  with  an  unconstrained, 
nonlinear  optimization  technique.  One  replication  of  the  algorithm 
determines  the  optimal  DGZ  locations  in  a  finite  number  of  iterative 
steps.  Each  iterative  step  finds  improved  DGZ  locations  and  an  asso¬ 
ciated  larger  CEDF  value  as  compared  to  the  previous  locations  and  CEDF 
value.  The  algorithm  iterates  until  no  significant  increase  in  the  CEDF 
is  possible.  This  chapter  explains  the  mathematical  formulation  of  the 
CEDF  and  its  gradient.  Chapter  III  presents  the  unconstrained,  non¬ 
linear  optimization  techniques  used  to  maximize  the  CEDF. 


Conceptual  Model 

The  CEDF  is  a  function  of  weapon  and  installation  parameters  and 
the  coordinates  of  the  DGZs.  The  conceptual  model  of  the  CEDF  is  shown 
in  Figure  2.  The  i  subscript  of  either  a  variable  or  a  parameter  refers 
to  one  of  the  m  weapons;  the  j  subscript  refers  to  one  of  the  n  instal¬ 
lations  in  the  target  complex.  All  parameters  are  constants  (either 


prespecified  or  calculated  values)  except  the  (X^,Y^)  DGZ  coordinates 
for  each  of  the  m  weapons.  The  basic  element  of  the  the  CEDF  is 


the  Fd.  .  —  the  probability  of  achieving  a  specified  level  of  damage 

l,  J 

to  installation  j  from  weapon  i.  Similarly,  the  expected  damage  to 


installation  j  from  weapon  i  is  DE.  .  —  the  product  of  Pd,  .  and 


the  probability  of  arrival  for  weapon  i,  Pa^.  The  cumulative  expected 

damage  to  installation  j  from  all  m  weapons  is  DE^ .  This  formulation 

for  determining  the  cumulative  damage  to  an  installation  from  multiple. 

bursts  is  similar  to  the  formulation  used  in  SIDAC  and  NUCWAVE 

(Ref  19:A-4  and  25:A-2).  For  each  installation,  the  Installation 

Expected  Damage  Function,  IEDF,  is  the  product  of  its  value  v,  and  DE,. 

J  J 

The  CEDF  is  an  unconstrained,  nonlinear  function;  it  is  the  summation 

of  n  lEDFs.  The  function  is  nonlinear  because  Pd.  .  is  a  nonlinear 

*»  J 

function.  Again,  all  CEDF  parameters  are  constants  except  the  (X^,Y^) 
DGZ  coordinates  for  each  of  the  m  weapons. 

In  order  to  determine  the  CEDF,  the  algorithm  requires  scenario 
dependent  inputs,  installation  and  weapon  parameters.  The  minimum 
necessary  installation  parameters  include: 

1.  The  nunber  of  installations  in  the  target  complex  -  n 

2.  The  coordinates  of  each  installation  -  (xj , y^ ) 

3.  A  VNTK  code  for  each  installation^  indicating  the  installation's 
susceptibility  to  blast  damage 

4.  A  value  from  a  relative  installation  value  system  -  v. 

The  minimum  necessary  weapon  parameters  include: 

1.  The  number  of  weapons  available  -  m 

2.  The  height  of  burst  for  each  weapon  -  HOB^ 

3.  The  yield  in  kilo tons  for  each  weapon  - 

4.  The  CEP  for  each  weapon  -  CEP^ 

5.  Hie  probability  of  arrival  for  each  weapon  -  Pa^ 

6.  The  initial  DGZ  locations  prior  to  optimization  -  (X^>Y^) 

The  assignment  of  specific  numeric  values  to  these  parameters  was 

not  a  critical  element  of  the  study.  Consequently,  several  hypothetical 
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target  complexes  were  used.  These  complexes  are  described  in  Chapter  IV. 

Two  mathematical  forms  of  the  CEDF  are  used;  hence,  there  are  two 
parallel  algorithms,  one  for  each  form  of  the  CEDF.  The  first  CEDF  is 
the  CEP-Excluded  version;  the  second  CEDF  is  the  CEP- Included  version. 
These  two  versions  are  explained  in  the  next  section.  The  only  differ¬ 
ence  between  the  two  CEDFs  is  their  respective  forms  of  the  probability 
of  damage  function.  The  Defense  Intelligence  Agency  (DIA)  developed 
these  Pd  models  to  provide  analytical  approximations  to  actual  blast 
damage  data.  The  CEP-Excluded  CEDF  uses  a  closed  form  analytical 
expression  of  an  installation's  Fd  function  that  is  independent  of 
weapon  delivery  system  accuracy,  that  is,  weapon  CEP  -  0.  The  CEP- 
Included  CEDF  uses  a  more  complicated  analytical  expression  of  an 
installation's  Pd  function  that  includes  weapon  CEP. 

The  two  CEDF  forms  are  used  for  three  reasons.  First,  a  closed 
form  analytical  expression  for  the  gradient  of  the  CEP-Included  CEDF 
expression  was  not  available;  hence,  this  CEDF  could  only  be  maximized 
with  an  optimization  technique  that  used  function  values.  However, 
since  gradient  optimization  techniques  are  generally  more  efficient 
than  function  value  techniques  (Ref  2:152;  5:321;  and  10:286),  a  second 
form  of  the  CEDF  is  desired.  Therefore,  an  analytical  expression  for 
the  gradient  of  the  CEP-Excluded  CEDE'  is  calculated.  This  CEDF  is 
maximized  using  gradient  optimization  techniques.  Chapter  III  explains 
different  optimization  schemes  and  the  optimization  techniques  used  to 
maximize  the  two  CEDFs.  The  second  reason  for  using  two  CEDFs  is 
verification.  The  results  of  the  algorithm  are  compared  to  insure 
that  they  provide  the  same  DGZ  locations  and  complex  expected  damage 
value.  The  last  reason  for  using  two  CEDFs  is  to  investigate  the 
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effect  that  the  assumption  of  CEP  -  0  has  on  the  location  of  the 

optimal  aimpoints. 


Probability  of  Damage  Models 

The  CEDF  is  an  unconstrained,  nonlinear  function  of  2m  variables  •« 
the  (X^Y^)  DGZ  coordinates  for  each  of  the  m  weapons. 


n  n  r  m 

CEDF  -  Z  IEDF.  -  Z  v,  *  1  -  TT(1  -  Pd,  .  *  Pa*) 

j-1  J  j-1  J  L  i“l  i,J 


(1) 


The  basic  element  of  the  CEDF  is  Fd^j  —  the  probability  of 
achieving  a  specified  level  of  damage  to  Installation  j  from  weapon  i. 

Pd^  j  is  a  function  of  two  independent  variables,  the  (X^,Y^)  coordi¬ 
nates  of  weapon  i.  Two  forms  of  the  Pd.  .  function  used  in  this  study 
are  part  of  the  D1A  Riyscial  Vulnerability  (FV)  System.  They  are  not 
independent  formulations.  These  formulations  are  described  in 
Mathematical  Background  and  Programing  Aida  for  the  Physical  Vulner¬ 
ability  System  for  Nuclear  Weapons  (Ref  6).  Therefore,  only  a  useable, 
but  limited,  description  will  be  presented  here. 

The  Pd.  .  depends  on  the  known  distance  s  between  DGZ.  and 
X»J  i 

installation  j.  The  coordinates  of  installation  j  are  (xj,y^).  The 
geometry  of  the  installation-weapon  interaction  is  shown  in  Figure  3. 

The  algorithm  uses  a  flat  earth  approximation  to  calculate  this  distance, 
that  is, 


s  - 


(2) 
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Figure  3.  Geometry  of  the  installation-weapon  interaction  (Ref  6:20). 

The  (X^,Y^)  and  (Xj,Yj)  coordinates  are  measured  in  feet  from  a  common 
origin.  Comparatively,  the  distance  r  is  the  distance  between  instal¬ 
lation  j  and  the  actual  weapon  impact  point.  When  CEP  *  0  ,  the 

impact  point  and  DGZj^  may  not  be  the  same  point,  and  s  will  not  equal  r. 

CEP-Excluded  Model.  If  the  CEP  of  the  weapon  delivery  system  can 
be  assumed  to  equal  0,  that  is,  the  actual  weapon  impact  point  is  the 
DGZif  then  the  distance  r  from  installation  j  to  the  impact  point  is 
known.  The  distance  r  can  be  calculated  from  the  (Xi,Yi)  DGZ  coordinates, 
the  (Xj>yj)  installation  coordinates,  and  Eq  (2).  The  distance  damage 
function,  i»  the  DIA  analytical  approximation  for  the  probability 

of  damage  function  when  weapon  CEP  -  0.  It  is  based  upon  actual  blast 
damage  data.  Pd(r)  is  the  complement  of  the  cumulative  log  normal 
distribution  function.  For  this  CEDF  version,  CEP-Excluded,  the  proba¬ 
bility  of  damage  function,  Pd^j,  is  the  distance  damage  function,  Pd(r). 
However,  it  will  be  referred  as  the  distance  damage  function  Pd(r)  to 
parallel  the  DIA  development.  The  shape  of  a  Pd(r)  function  is  shown 
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Figure  4.  A  probability  of  damage  function  P^(r) 

in  Figure  4.  The  independent  variable  is  the  distance  r  between  the 
installation  and  the  impact  point. 

The  distance  damage  function  ia  a  nonlinear  expression  in  integral 
form;  it  is  specified  by  the  location  and  dispersion  parameters,  <■*  and  -4  . 


The  parameter  t  is  a  dummy  variable  of  integration  for  the  normal  proba¬ 
bility  distribution.  '*Ihe  parameter  ot  ,  which  is  the  median  of  the  log 
normal  density  function,  is  the  distance  from  ground  zero  at  which  there 
is  a  50X  chance  of  achieving  a  specified  level  of  damage.  The  parameter  £ 
is  the  standard  daviation  of  ln(r)"  (Ref  6:6,7).  These  parameters  depend 
on  the  weapon  radius  (WR)  and  distance  damage  sigma  (cr^).  If  a  weapon 
is  datonatad  within  a  uniform  distribution  of  targets,  then  the  WR  is 
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the  radius  of  a  circle  centered  at  the  weapon  impact  point.  The  circle 
contains  as  man/  targets  undamaged  to  a  specified  level  inside  the 
circle  as  there  are  targets  damaged  to  a  specified  level  outside  the 
circle,  "o',  is  a  measure  of  the  variance  of  the  density  function.  A 
small  cru  Indicates  a  relatively  rapid  fall  off  of  the  damage  function; 
a  large  od  indicates  a  more  gradual  fall  off"  (Ref  6:11). 

Prior  to  1  September  1972,  the  analytical  approximation  of  actual 
blast  damage  data  was  the  circular  coverage  function  with  parameters  UR 
and  ,  However,  before  that  data,DIA  decided  that  the  distance  damage 
function  with  parameters  %  and  £  provided  a  better  fit  to  actual  blast 
damage  data.  Since  previously  measured  and  calculated  target  vulner¬ 
ability  data  depended  on  UR  and  a;, ,  DIA  developed  mathematical  trans¬ 
formations  to  determine  ot  and  p  from  UR  and  <rA  . 

P  ■  j-ln<l  -oj2) 

c*  •  URe  ^ 

With  these  transformations,  the  distance  damage  function  could  specify 
the  Pd  for  targets  characterized  by  the  Physical  Vulnerability  (FV) 
coding  system. 

Consequently,  the  probability  of  damage  to  installation  j  from 
weapon  i,  Pd^  j,  can  be  calculated  using  Eqs  (3)  and  (4),  after  crd  , 

UR,  and  r  have  been  determined.  UR  and  qj  are  parameters  that  are 
calculated  using  prespecified  user  values.  UR  depends  on  the  weapon's 
yield  and  HOB  and  the  installation's  VNTK  code.  Hence,  there  is  a 
unique  UR  for  each  weapon  i-installntion  j  interaction  —  WR,  , . 
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Likewise,  and  {3  depend  only  on  the  installation's  VNTK  code.  Hence, 

there  is  a  unique  and  6  for  each  installation  —  .  and  0;  . 

Appendix  A  presents  the  calculation  of  ^  and  WR.  The  independent 

variable  r  is  actually  a  function  of  two  independent  variables,  the 

(Xi,Y.)  DGZ  coordinates,  and  two  constants,  the  (x.,y.)  installation 
1  J  j 

coordinates. 

The  F^(r)  cannot  be  expressed  in  closed  form  in  terms  of  elementary 
functions;  however,  it  can  be  calculated  by  use  of  the  error 
function,  erf(u)  (Ref  6:21  and  1:298).  The  erf(u)  specifies  the  proba¬ 
bility  that  a  standard  normal  random  variable  observation  is  within  ±  u 
of  the  mean  value. 


Pd(r)  -  Wi>j(Xi,Yi)  -  0.5  +  0.5  erf  (^^*)  for  z(r)  >  0 


■  0.5  -  0.5  erf 


IkklL 

,  jX  , 


for  z(r)  <  0 


where 


z(r)  -  zCXpY^  - 


(5) 


(6) 


A  polynomial  function  of  the  independent  variable  u  can  approximate  erf(u) 
(Ref  14:185). 

In  summary,  if  the  CEP  of  the  weapon  can  be  assumed  to  equal  0, 

then  the  probability  of  damage  Pd,  .  to  installation  j  from  weapon  i 

J 

can  be  calculated  using  Eqs  (5)  and  (6).  Prespecified  target  and  weapon 
parameters  are  necessary  to  calculate  WR^  j,  jy  and  r(X^,Y^). 
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CEP-Included  ModtJ  ’f  the  CEP  of  the  weapon  delivery  system 
cannot  be  assumed  to  equal  Oi  that  is,  the  actual  impact  point  of  the 
weapon  is  unknown,  then  the  distance  r  from  the  impact  point  to 
installation  j  is  unknown.  The  geometry  of  the  installation-weapon 
interaction  is  shown  in  Figure  3.  The  distance  s  from  DGZ^  to  instal¬ 
lation  j  can  be  calculated  from  the  (X^,Y^)  DGZ  coordinates,  the 
(xj,yj)  installation  coordinates,  and  Eq  (2).  The  unknown  distance  ? 
from  DGZ^  to  the  actual  impact  point  is  a  function  of  s  and  the  inde¬ 
pendent  variables  r  and  -6-  • 

The  DIA  model  determines  the  probability  of  damage  to  installation  j 
in  the  following  way.  First,  for  each  possible  impact  point,  the  proba¬ 
bility  of  damage  is  multiplied  by  the  probability  that  the  weapon  arrives 
end  detonates  at  that  point  (Ref  6:19).  The  sum  of  these  products  for 
all  possible  impact  points  specifies  the  probability  of  achieving  the 
desired  level  of  damage  to  installation  j  from  weapon  i,  Pd^  y  This 
summation  is  a  multiple  Integral  over  the  area  that  contains  all  possible 


where  P^(r)  ■  distance  damage  function,  Eq  (3) 
O'2  -  CEP/1.1774 

*"|r2  +  s2  -  2  rs  cos# 


For  this  CEDF  model,  Pd^  j  and  P^(r)  are  not  the  same  function.  Pd^  ^ 
has  two  distinct,  yet  dependent,  terms:  Pd(r)  and  — e  2 o' 
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(l)  r  (r)  specifies  the  probability  of  damage  to  an  installation 
from  an  impact  point,  and 


(2) 


specifies  the  circular  normal  probability  of 


the  weapon  arriving  and  detonating  at  the  DGZ. 

A  closed  form  solution  to  Eq  (7)  does  not  exist;  however,  an 
analytical  approximation  does  (Ref  6:23). 


b 

(^f(r)dr 

a 


The  limits  of  integration,  a  and  b,  are  selected  such  that  when  r  <  a 
or  r  >  b,  f(r)  -  0.  They  ara  functions  of  s,  CEP,  WR,  and  o' . 

Appendix  B  presents  the  development  of  f(r),  the  determination  of 
a  and  b,  and  the  calculation  of  Fd^  j.  The  function  f(r)  has  two  dif- 
(•cent  form.  fo™  d.t*nd.  on  th.  disunc.  .  betw**,  DG^  and 

installation  J,  the  distance  r,  and  the  weapon's  CEP. 

This  integral  can  be  evaluated  using  Gauss-Lagendra  quadrature, 
a  numerical  integration  technique.  This  technique  approximates  a  definite 
integral  as  a  finite  series.  Each  term  in  the  series  is  a  weighted 
function  value. 


(b  -  a) 

2 


ww^  *  f(rk) 


where 


(b-.a)  +  (ajLb) 

2  Z2k  +  2 


(8) 

(9) 


Gauss-Legendre  quadrature  differs  from  the  more  common  trapezoidal 
numerical  integration.  In  Gauss-Legendre,  the  distances  between 


TABLE  I 


Steps  to  Calculate  Pd^  j 


Step 

Given  Parameters 
or  values 

Results 

1 

VNIK 

"v 

2 

3 

VNTtC,  Yieldi ,  hobi 

“i.j 

4 

V  Yi»  xj»  yj 

r(XjL,Y^),  s(X^,Y^) 

5 

h*  r 

Pd,  ,  (CEP-Excluded) 

5 

Ojj,  WRj^j,  s,  CEP^ 

a,  b 

6 

V®i, j ,  r,  s,  a,  b,  CEP^ 

Pd1>j  (CEP-Included) 

the  values  along  the  abscissa  are  not  equal.  The  values  of  the  quadra¬ 
ture  coefficients,  ww^,  and  the  base  points,  zzk,  can  be  determined  from 
the  N^1  Legendre  polynomial.  Gauss-  Legendre  quadrature  is  discussed  in 
more  detail  in  Appendix  C. 

In  summary,  if  the  CEP  of  the  weapon  cannot  be  assumed  to  equal  0, 
then  the  probability  of  demage  to  installation  j  from  weapon  i,  Pd^  j , 

>an  be  calculated  using  Eqs  (8)  and  (9).  Specific  target  and  weapon 
parameters  are  necessary  to  calculate  a,  b,  and  f(r, ). 

K 

Each  Pd,  ,  is  an  integral  part  of  a  Complex  Expected  Damage 

X,J 

Function  (CEDF).  Eqs  (5)  and  (8)  are  used  to  calculate  Pd.  ,  for  the 

J 

ChP-Exclud*  '  and  for  the  CEP- Included  CEDE  models.  Table  I  lists  a 
sutvrwy  of  the  steps  necessary  to  calculate  Pd^  y 
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The  CEDF  Model 


The  CEDF  is  an  unconstrained,  nonlinear  function  of  the  (X^,Y^)  DGZ 
coordinates  for  each  of  the  m  weapons. 

n  m 

CEDF  -  Vj  *  [l-TT(l-Wi  *  Pa±)  ]  (1) 

J-l  L  i-1 

The  i  subscript  refers  to  one  of  the  m  weapons;  the  j  subscript  refers 
to  one  of  the  n  installations.  The  2m  independent  variables  of  the  CEDF  are 

<X1-  V  *2'  V  -•  V  V 

The  2m  elements  of  the  DGZ  coordinate  vector,  &,  are  these  2m  variables 
in  a  revised  order. 

x  "  (Xp  *2’  •••»  Yl»  Y2»  ”*»  V 

The  (X^,Y^)  DGZ  coordinates  of  weapon  i  are  (x^x^).  Similarly,  the  2n 
parameters  specifying  the  (xj,yj)  coordinates  of  the  n  installations  are 

(xx.  yv  Xj,  y2>  x^,  yn) 

The  2n  elements  of  the  installation  coordinate  vector,  xx,  are  these  2n 
parameters  in  a  revised  order. 

xx  ■  (xj,,  X2»  ...»  Xp,  yj_,  y2»  •••»  yn) 
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DE1,1  -  Ml,l  *  Pal 


DE1,3  ’  Pdl,3  *  Pal 


DE2,1  *  M2,l  *  Pa2  DE2,3  -  Pd2,3  *  Pa2 

DE1  -  1  -  <1  -  DE^  ])(1  -  DEj  y)  DEj  -  1  -  (1  -  DE^  3>(1  -  DEj  3) 

>  (x*3!»6) 

(xx  „j  A  (X2  C>  >  3  * 


IEDFl  •v1*  DEl 


(xXjiXX^ 


<Va.'> 


l3’*V 


IEDF^  ■  *  DEg 


^xl* X3)  Pa^  ^X2,x4^  ^a2 


(xx2,xx5)  ^ 

DE1,2  “  **1,2  *  Pal 

DE2  2  ■  Pd2  2  ^ 

DE2  “  1  -’(1  ‘  DE1(2)(1  *  DE2,2> 

IEDF2  "  v2  *  DE2 


installation 

DGZ 


2  Figure  5.  A  representative  CEDF  with  three 

installations  and  two  weapons. 

The  (x.,y.)  coordinates  of  installation  j  are  (xx.,xx,.  )• 

J  J  J  J411 

Like  the  Pd,  , ,  the  damage  expectancy  for  installation  j  from  weapon  i 
is  a  function  of  two  independent  variables  DE^  jU^x^).  Similarly, 
like  the  CEDF,  the  cumulative  damage  expectancy  for  installation  j  from 
all  weapons  is  a  function  of  2m  independent  variables  DEj(x). 

A  pictorial  description  of  a  representative  CEDF  is  shown  in 
Figure  5.  The  target  complex  consists  of  three  installations  and  two 
weapons.  Each  installation  is  characterized  by  a  value  —  v^,  a  VNTK 
code,  and  (xx^xx^)  coordinates.  Each  weapon  is  characterized 
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by  a  yield  —  Y^,  height  of  burst  —  HOB^,  probability  of  arrival  —  Pa^, 
and  (x.,xi+2).  In  order  to  determine  all  of  the  DE^  .s,  six  Pd^  .s  are 
calculated  according  to  the  steps  in  Table  I,  one  for  each  combination 
of  i  -  1,  2  and  j  -  1,  2,  3.  For  the  weapon  coordinates  (x^jX^+2), 
i  -  1,  2 

CEDF(x)  -  vl  *  DEl  +  v2  *  DE2  +  v3  *  DE^  (1) 

However,  if  the  weapon  coordinates  are  changed  to  (x^1 ,x^+2' ), 
i  -  1,  2,  then  DEj  may  change  for  each  of  the  j  installations.  That  is, 
DEj(x)  may  not  equal  DE^ (x' )  for  all  j.  If  this  is  true,  then  the  CEDF(x' ) 
may  be  either  greater  than,  equal,  or  less  than  CEDF(x). 

In  order  to  maximize  the  CEDF,  it  is  necessary  to  locate  the  x* 

DGZ  coordinates  such  that 

CEDF(x*)  >  CEDF(;g)  for  all  x 

One  of  the  optimization  techniques  used  to  maximize  the  CEDF(x)  and  to 
locate  the  optimal  x*  DGZ  coordinates  required  the  gradient  of  the  CEDF(x). 

^  A* 


Gradient  of  the  CEP-Excluded  Model 

The  Complex  Expected  Damage  Function  (CEDF)  is  a  nonlinear  func¬ 
tion  of  2m  variables  —  the  (x^x^)  DGZ  coordinates  of  the  m  weapons. 


CEDF(x) 


(1 "  *  PV 


(i) 


All  parameters  of  the  CEDF(x)  are  constants  except  the  probability  of 

29 


■ L ■  ■ 


achieving  a  specified  level  of  damage  to  installation  j  from  weapon  i, 
Pd.  ..  Each  Pd.  .  is  a  function  of  two  independent  variables,  the 
(xi»xi+m)  coordinates  of  weapon  i.  Therefore,  to  calculate  the 
gradient  of  the  CEDF,  the  Pd^  .  must  be  differentiable  with  respect  to 
the  two  independent  variables.  A  closed  form  analytical  expression  for 
the  gradient  of  the  CEP-Included  Pd.  .,  Eq  (7),  was  not  available. 

This  was  one  of  the  reasons  for  formulating  the  second  version  of  the 
CEDF,  the  CEP-Excluded  model*  This  section  presents  the  calculation  of 
the  gradient  of  the  CEP-Excluded  CEDF. 

The  gradient  of  the  CEDF(x)  is  a  vector  of  2m  element.. . 


grad(CEDF)  ■  3CEDFfo).  .(^ 

35  V3x1  ,3x2 


where  CEDF(x)  -  v^DE^x)  +  v2»DE2(x)  +  ...  +  vn*DE^(x) 


(10  J 


and  DE 


;j(x)  -  1  -  (1  -  Pa^.Xl  -  Pa2Pd2)j)  ...  (1  -  Pa^j) 


Since  v 


J 


is  a  constant,  the 


element  of  grad(CEDF)  is  of  the  form 


d  CEDF 


+  v2* 


3°E2 

dxk 


+  ...  +  V  * 

n 


ajfn 

dxk 


j-1 


v.* 

J 


(11) 
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Each  of  the  n  terms  of  Eq  (11)  is  of  the  form 


I 


V.*  “  -  v,*| 

J  e)xk  J 


0(1  -  Pa, Pd,  ,)  v 

-<i  -  iVdi,j)*  -  vy 

Is 


3(1  -  PajPd.  ,) 

+  ..+  (1  -  Pa^Pd^^)*  ♦  - q  ^  — 4"***d  “ 

0(1  -  PaJPcL  ,)1 

+  ..+  (1  -  )*.••(!  "  | 


However,  since  each  Pd^  j  is  a  function  of  only  two  variables  (x^x^^) , 
for  all  i,  all  "  ^1^1, terms  equal  0  except  for  k  -  i  and  i  +  m. 

3xk 

Hence,  for  k  -  i  and  i  +  m 

ADE.  0(1  -  Pa<Pdi  ,) 

VjW  "  ’  VU  '  PalWl.j)*"' - d\  ~  ‘  PamPdm,j) 


»  Vj*(l  -  PajPd^  j)**’*  Pai 


3*k 


*...(1  -  P^Fd^) 


(12) 


Now  define 


m 

factor(j)  “  v.  TT  (1  -  P®!1^!  0 


and  rewrite  Eq  (12)  as 


f actor( j ) 

(1  -  * 


3Wi,J 
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for  k  ■  i  and  i  +  m  (13) 


Hie  gradient  of  the 
of  the  element. 
Eq  (13)  is  the  form 


of  a*T 


for  k  ■ 


CEDF  is  a  vector  of  2m  elements.  Eq  (11)  is  the  form 
Similarly,  each  element  is  a  summation  of  n  terms, 
of  each  of  these  n  terms.  Analytical  expressions 

i  and  i  +  m  are  needed  to  completely  specify  the 


gradient  of  the  CEDF(x). 


The  CEP-Excluded  version  of  the  Pd.  .  is 

J 


Pdi>j(xi,xi+m)  -  0.5  +  0.5  erf  for  z  >  0  (5) 


-  0.5  -  0.5  erf  (-i|U 

<rr  / 


for  z  <  0 


where 


z(r)  -  -x-  In 


WR 


iii 


e 


-Pj 


Ti 

3Mi  t 

First,  u  will  be  calculated 
Oxl 


i  ft-2 

-  -g-  ln(WR^^ j6  )  -  ln(r)  (6) 


Pj 


for  z  >  0,  that  is,  for  r  <  WR^ 


M 

Let  u  -  rrr  and  use  the  chain  rule 


JT 


where 


.05  3«f(u)/3u  \ 

3u  Ux^ 


(14) 


SjL.  1  31.1 
9  xi  JT  d  xi 


(15) 
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1/2 


(2) 


Now  since  r 


(x.  -  xx.)" 
i  J 


(x.  -  xx.,  y 

i+m  j+n 


and  z(r)  *  z(x^,x^4m)  from  Eq  (6),  and  again  using  the  chain  rule 


Slil  ■  (16) 

3X1  3r  3x^ 

where 

3Jzj  .  _1_M\ _ L 

dr  Pjr 


and 


(17) 


hence 


(18) 


The  derivative  of  erf(u)  from  the  Handbook  of  Mathematical  Functions  with 
Formulas ,  Graphs,  and  Mathematical  Tables  (Ref  1:298,801)  is 


3«rf(u)  .  B-u2 

3u  JtT 


(19) 


Hence,  combining  Eqs  (14),  (18),  and  (19)  specifies 


aWi.j  - 1  e“u^ 
2  W 
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for  z  >  0 


(20) 


where 


u  ■ 


A  similar  mathematical  development  was  used  to  calculate 


^  xi+m 


P?  ^J+n  "  xi+m) 


.2 


For  z  <  0,  that  is  r  >  W lks  -e  ^  , 

A»j 


then 


ffikl  -  -0.5  3-gSEM  ^iL 

3^  3  u  3 


(21) 


(22) 


(23) 


The  only  difference  between  this  development  and  the  previous  development 

for  z  >  0  is  the  sign  of  SJjEl  .  This  partial  derivative  is  positive 

3  r 


because 


r  >  WR 


and  | z |  is 


hence, 


3jzJ  „  1 

3r  fjr 


(24) 


Combining  Eqs  (IS),  (16),  (17),  and  (24)  yields 

^Xi  f2  \Pjr/\  r  / 


(25) 
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Now  Eqs  (19),  (23),  and  (25)  specify 


**LX 


i.  ^  i  m  p 

2  <|tT  J2  Wr/V 

J 


X.  -  XX 


for  z  <  0 


J2?^  <XX3  "  Xi) 


Hence,  compering  Eqs  (20)  and  (26),  >1-  is  the  seme  for  all  z. 


A  similar  development  indicates  -» . is  also  the  same  for  all  z . 

Oxi+m 


In  summary,  the  gradient  of  the  CEP- Excluded  CEDF(x)  is  a  vector 

elk 

of  2m  elements.  The  k  element  of  the  gradient  is 


3  CEDF  3-  *3^1 

3*k  p!  J  3xk 


where 


v  *  %T^  "  ) — — ^ - *  Pa,  *  for  k  -  1  and  i  +  m  (13) 

i  BXfc  (l-PalWi>j)  i  3^ 


Also,  for  k  -  i 


and  for 


whart 


and 


2!Sa. 

3xk 


(20) 


■  i  +  m 


a  Wi.i 

3-k 


W 


(22) 


u-  1*L 

JT 


(21) 


[(V 


XX j )2  + 


<x 


i-Hn 


(2) 
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III.  Optimization  of  the  Complex 
Expected  Damage  Function 

Two  approaches  are  available  to  maximize  a  function  of  n  variables — 
analytical  and  numerical  search.  An  analytical  approach  is  preferred  if 
the  roots  of  the  n  equations  defining  the  critical  points  of  the  function 
are  easily  determined  and  solved.  These  equations  are  the  first  partial 
derivatives  of  the  function  set  equal  to  0.  However,  if  these  analytical 
expressions  are  not  easily  determined  or  solved,  then  numerical  search 
techniques  are  necessary  to  determine  the  maximum  of  the  function. 

Numerical  search  techniques  require  an  organized,  efficient  exploration 
of  the  solution  space. 

Numerical  search  techniques  were  used  to  maximize  the  Complex  Expected 
Damage  Function  (CEDF)  because  of  the  complexity  of  the  CEDF.  This  chapter 
presents  a  general  methodology  and  overview  of  numerical  search  techniques 
that  are  used  to  maximize  unconstrained  functions.  Also,  this  chapter 
discusses  the  two  related  techniques  that  were  used  to  maximize  the  two 
versions  of  the  CEDF  —  the  CEP -Included  and  the  CEP-Excluded  models. 

The  primary  difference  between  the  two  CEDF  models  is  that  there  was  a 
closed  form  expression  for  the  gradient  of  the  CEP  Excluded  model. 
Therefore,  gradient  search  techniques  could  be  used  to  maximize  the  CEP 
Excluded  model.  The  CEP -Included  model  was  maximized  using  Powell's 
conjugate  directions  method  (Ref  23).  This  method  maximizes  a  function 
using  only  function  values.  The  CEP-Excluded  model  was  maximized  using 
a  conjugate  gradient  with  restarts  method  (Ref  24). 
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1 


Optimization 


Optimization  is  a  process  that  attains  the  best  or  most  effective 
results  for  a  problem,  while  satisfying  any  given  conditions  or  con¬ 
straints.  Optimization  can  be  either  maximization  or  minimization.  One 
part  of  this  study  was  to  maximize  the  CEDF(x),  a  nonlinear  function  of 
2m  independent  variables  —  the  (x^x^)  DGZ  coordinates  of  each  weapon. 
The  CEDF(jg)  is  an  unconstrained  function.  It  can  be  maximized  by  mini¬ 
mising  -CEDF(x).  that  is,  the  point  in  2m  spAce,  such  that  CEDF(x) 
is  a  maximum,  is  the  same  point  where  -CEDF(x)  is  a  minimum. 

In  this  chapter,  direct  references  to  maximizing  the  CEDF  are  not 
made.  Instead,  all  references  concerning  optimization  reference  mini¬ 
mizing  an  unconstrained,  nonlinear  function  of  n  variables,  f(x) ;  %  is 
an  n  element  vector  in  n-dimensional  space,  Rn.  The  gradient  of  f (5) 
is  Vf(j$);  the  Hessian  matrix  of  £(3)  is  H(^).  A  base  point  in  Rn 
is  #*■}  the  optimal  point  in  Rn  is  3*. 

There  is  an  important  difference  between  a  strict  local  minimum  and 
the  global  minimum  of  £(x).  The  following  two  definitions  are  extracted 
from  Avriel  (Ref  3:10),  A  real  valued  function  f(x)  with  domain  D  in 
Rn  has  a  strict  local  minimum  at  x*,  if  there  exists  a  number  6  such 
ttiat 

f(x*)  <  f(x)  for  all  %  €  D  (27) 

such  that  |  x  -  £*|  K  $  •  This  definition  states  that  x'v  is  a  local 
minimum  over  a  region  bounded  by  a  number  S  .  If  Eq  (27)  holds  for  al.l 
x  6  D,  that  is,  £  not  contained  within  a  bounded  region,  then  x*is  the 
global  minimum.  Optimization  techniques  locate  the  global  minimum  only 


38 


under  special  conditions.  That  is,  the  function  is  known  to  be  unimoda'l. 
Generally,  if  a  function  is  not  known  to  contain  a  global  minimum,  an 
accepted  procedure  is  to  search  D  from  a  number  of  initial,  separated 
base  points  to  determine  all  local  minimums.  Beveridge  and  Schechter 
state  "the  only  method  of  determining  the  global  optimum  is  the  direct 
comparison  of  the  function  values  at  various  local  optima"  (Ref  5:357). 

Numerical  Search  Techniques 

The  numerical  search  for  the  minimum  value  of  an  unconstrained 
function  f(x)  with  domain  D  in  Rn  is  a  sequential,  iterative  process. 

It  includes  the  successive  calculation  of  new  objective  function  values, 
f  ($*•),  and  the  comparison  of  these  values  with  the  best  value  that  has 

JU 

bean  obtained  so  far.  It  is  necessary  to  determine  £  by 

f(x*)  >  f(x2)  >  ...  >  f(^)  >  ...  >  £(&*)  for  all  £  €  D 

While  generating  the  sequence  of  x,S  each  unconstrained  numerical 
search  technique  must  consider  three  important  elements  —  the  search 
direction,  the  distance  to  move,  and  the  stopping  criteria.  From  a 
base  point  x*,  n  search  technique  must  select  (1)  a  direction  of 
movement  £  and  (2)  a  distance  to  move  t.  These  values  specify  the 
next  point  in  Rn 


+  td 


If  f(3S1+1)  <  fOc1),  then  3i+l 


is  a  bettor  estimate  of  the  local  minimum 
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than  x1.  The  stopping  criteria  fot  a  .’-oardi  tecnnique  depends  upon  the 
values  of  either  x^^x^),  or  Vf(x/+“),  Th  t  is,  if  either 

|&i+1  -  jg1!  4  &v  or  jf(xi+1)  -  tfx1)!  *  i 2 »  or  I  Vf(xi4l)|  <«S3, 

ju  4  4*1 

then  the  technique  stops  iterating,  and  x  ■■  x  is  the  optimal 
point  in  Rn  such  that  f(x,^)  is  a  minimum.  Numerical  search  techniques 
use  different  methods  to  determine  d  and  t. 

There  are  three  categories  of  numerical  search  techniques  (Ref  2:101), 
The  first  category  includes  direct  search  techniques.  These  techniques 
use  only  functional  values  to  locate  x^  from  x*.  The  second  category 
includes  gradient  or  first-order  search  techniques.  These  techniques 
use  f(x^)  and  Vf(x^)  to  determine  .  Generally,  gradient  methods 
are  more  efficient  and  preferred  to  direct  techniques  (Ref  2:152; 

5:321}  and  10:386).  However,  when  the  gradient  is  not  easily  obtained, 
direct  searches  ere  more  appropriate.  The  last  category  includes 
second-order  techniques.  These  techniques  use  Vf(jt^) ,  and 

tine  Hessian,  H(x^),  to  locate  x^,  Detailed  explanations  of  the  fol¬ 
lowing  techniques  can  be  found  in  moat  optimization  books  (Ref  2;  3; 

5;  and  10).  Hence,  only  a  brief  explanation  is  presented  here. 

If  an  unconstrained  objective  function  f(x)  is  not.  easily  differ¬ 
entiated,  than  a  direct  search  technique  is  necessary  to  minimize  f(x)„ 

These  techniques  use  two  stages,  on  exploratory  and  a  pattern,  to  move 
i  i+1 

from  £  to  £  ,  Two  older  techniques  are  the  Hnoke-Jeaves  pattern 

search  and  Rosanbrock's  method  of  rotating  directions.  In  the  explor¬ 
atory  stages,  Hooke-Jeeves  only  searches  along  the  axial  coordinate 

directions;  Rosenbrock  searches  along  a  set  of  mutually  orthogonal 

i  i+l 

directions  that  are  determined  from  x  and  £  * .  Both  of  those  techniques 
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use  a  fixed  step  .length  when  exploring  around  x,'*'.  The  exploratory 
function  evaluations  specify  the  direction  d  of  the  pattern  move.  A 
more  efficient  technique  is  Powell's  method  of  conjugate  directions 
(Ref  23).  In  the  exploratory  stage  ,  Powell's  method  searches  along 
conjugate  directions  that  are  determined  from  x*  and  n  -  1  of  the 
previous  n  exploratory  search  directions.  Conjugate  directions  are  a 
generalization  of  orthogonal  directions.  Also,  Powell's  method  does 
not  use  a  fixed  step  length.  Rather,  this  method  conducts  a  one- 
dimensional  search  in  each  of  the  conjugate  directions  from  A 
more  complete  description  of  conjugate  directions  and  Powell's  method 
of  conjugate  directions  is  presented  later  in  this  chapter. 

Gradient  search  techniques  are  separated  into  two  categories, 
either  those  techniques  that  follow  the  gradient  as  closely  as  possible 
(the  methods  of  steepest  descant)  or  those  techniques  that  use  the 

gradient  to  guide  the  search  (the  conjugate  gradient  methods). 

< 

Cauchy's  steepest  descant  method  uses  the  gradient  to  find  the  direction 
of  greatest  functional  dacrease  from  a  base  point.  The  greatest  de¬ 
crease  in  f(jj)  is  in  the  direction  of  the  largest  negative  gradient. 

That  is, 

d  -  -  VfCfc1) 

3i+1  -  a1  -  *  Vf<£1) 

The  steepest  descent  method  uses  a  one-dimensional  minimization  search 

in  the  direction  of  -  Vf(x^)  to  determine  the  step  length  t  and  to 

i+1 

subsequently  locate  £  , 

min  ^x1  -  t  VfCx1)) 


Conversely,  the  method  of  conjugate  gradients  locates  a  new  base  point 

jjC  by  searching  along  a  mutually  conjugate  direction  d  (Ref  9  and  24) . 
The  direction  d  is  determined  using  the  gradient  at  the  current  base 
point  and  the  previous  search  direction.  From  the  method  uses  a 
one-dimensional  minimization  search  in  this  direction  to  determine  the 
step  distance  t. 

min  fCj1  -  tg) 
t 

1+1 

This  one-dimensional  search  establishes  a  new  base  point  x  .  A  more 
complete  description  of  the  conjugate  gradient  method  is  presented  later 
in  this  chapter. 

If  first  and  second  partial  derivatives  of  f(x)  are  available,  then 
Newton's  method  could  be  used  to  minimize  the  function.  This  technique 
usee  the  function's  gradient  and  Hessian  to  specify  the  direction  and 
the  distance  of  the  maximum  decrease  in  f(x). 

«i+1  -  x1  -  Htx1)-1  VfCx1) 

Avrial  states,  "If  there  are  a  large  number  of  variables,  the  function 
and  derivative  evaluations  and  especially  the  matrix  inversions,  are 
time-consuming  and  expensive  operations"  (Ref  3:288). 

These  are  not  the  only  techniques  available  to  minimize,  uncon¬ 
strained,  nonlinear  functions.  However,  they  are  representative  of 
the  three  categories  of  techniques  —  direct,  gradient,  and  second- 
order.  A  detailed  presentation  and  summary  of  numerical  search 
techniques  for  each  category  is  provided  by  Gill,  Murray,  arid 
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Wright  (Ref  10).  In  addition  to  the  above  techniques,  several  vari¬ 
ations  are  available  to  minimize  unconstrained  nonlinear  functions. 

The  most  powerful  is  the  "variable  metric"  or  quasi-Newton  method. 

This  algorithm  differs  from  Newton's  method.  It  does  not  use  the 
Hessian  matrix,  Instead  of  calculating  H(x^),  the  technique  approx¬ 
imates  the  inverse  of  H(x^)  by  using  the  gradient  and  the  previous 
estimate  of  the  inverse.  There  are  other  variations  of  Newton's  method. 
Similarly,  finite  difference  techniques  are  variations  of  gradient 
methods;  they  use  function  values  to  approximate  V^Cx^).  Generally, 
it  is  not  possible  to  single  out  a  method  as  the  one  to  be  used  in 
every  case. 

Each  form  of  the  Complex  Expected  Damage  Function  (CEDF)  was  maxi¬ 
mized  using  only  one  technique.  Since  an  analytical  expression  for  the 
gradient  of  the  CEP-Included  model  was  not  available,  it  was  maximized 
using  a  direct  search  technique  —  Powell's  method  of  conjugate  direc¬ 
tions.  Conversely,  since  an  analytical  expression  for  the  gradient  of 
the  CEP-Excluded  model  was  calculated,  it  was  maximized  using  a  gradient 
search  technique  —  a  conjugate  gradient  with  restarts  method. 

Conjugate  Directions  and  Quadratic  Termination 

Conjugate  directions  are  a  generalization  of  orthogonal  directions. 
A  set  of  n  vectors  d^,  d^,  dn  in  Rn  are  orthogonal  if  their  inner 
product  is  0,  that  is, 

djTdi  -  0  for  all  i  ¥  j  (23) 

A  set  of  n  vectors  is  mutually  conjugate  with  respect  to  the  n  *  n 
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symmetric,  positive  definite  matrix  A  if 

d,TAd.  -  0  for  all  i  *  j  (29) 

Thus,  for  every  nxn  symmetric, positive  definite  matrix  there  is  at  least 
one  set  of  n  mutually  conjugate  directions.  If  the  matrix  A  is  the 
identity  matrix,  then  Eq  (29)  becomes  Eq  (28),  the  definition  of  ortho¬ 
gonal  directions. 

Powell's  method  of  conjugate  direction?;  and  conjugate  gradient 
methods  depend  upon  the  concept  of  quadratic  termination.  Powell  proved 
the  following  theorem  on  quadratic  termination  (Ref  23). 

1  9  rr\ 

Theorem:  If  <£  ,  d  , . . . ,  d,  ,  m  <  n  are  mutually  conjugate 
directions,  then  the  minimum  of  the  quadratic  function  f(x) 
is  a  point  in  m-dimensional  space,  Rm,  containing  x°,  the 
initial  point,  and  the  directions  d\  <^,  ...  _dra,  and  the 
minimum  of  f($)  may  be  found  by  searching  along  each  of  the 
directions  only  once.  The  required  minimum  is  the  point 


m 


i-1 


The  distances  t,  are  determined  by  minimizing  f(x)  in  each 
direction 

oj 

m 

min  f(x°  +  ]>  t^1) 

1  i-1 

where  f(x)  -  xTAx  +  bx  +  c 

Of  Of  fV  >v,v 

and,  A  is  a  symmetric,  positive  definite  Hessian  matrix. 
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Powell's  theorem  proved  that  the  minimum  of  a  quadratic  function  f(x) 
with  domain  D  in  Rn  and  a  symmetric,  positive  definite  Hessian  could  be 
located  in  n  steps.  Each  step  is  a  search  along  one  of  the  n  mutually 
conjugate  directions  d*.  However,  since  each  direction  d*  has  n 
component  directions  in  Rn,  each  step  requires  n  one-dimensional 
searches  to  minimize  f(x). 

CEDF  Optimization  Methods 

Powell's  theorem  is  the  basis  for  the  method  of  conjugate  direc¬ 
tions  and  conjugate  gradient  methods.  If  f(x)  is  quadratic,  then  the 
minimun  can  be  located  in  a  finite  (<n)  number  of  steps.  However,  even 
if  f(x)  is  not  quadratic,  the  concept-  of  quadratic  termination  can  still 
be  used  to  locate  the  minimum.  When  the  method  is  applied  to  non¬ 
quadratic  functions,  it  becomes  iterative  and  a  test  of  convergence  is 
necessary  to  determine  die  minimum  of  f(x).  This  section  presents  a 
brief  explanation  of  these  two  optimization  techniques.  Detailed 
explanations  of  them  are  available  in  Refs  9,  23,  and  24;  also,  most 
optimization  books  provide  complete  explanations  of  these  techniques. 

Powell's  Method  of  Conjugate  Directions .  This  section  presents 
an  algorithm  for  Powell's  method  of  conjugate  directions  (Ref  17  and  23). 
This  method  assumes  quadratic  convergence  of  f(x);  the  method  will  not 
locate  the  local  minimun  in  n  steps  unless  the  f(x)  is  quadratic. 

Instead,  the  method  iterates  f rom  x*  to  until  -^|  , 

In  this  development,  the  superscript  i  refers  to  the  iteration  and  the 
subscript  j  refers  to  one  of  the  n-dimensional  component  directions 
of  Rn.  The  starting  point  in  Rn  is  x°;  the  initial  search  directions  d| 
are  the  Rn  coordinate  directions,, 
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An  iteration  process  is  used  to  locate  x*  such  that  f(x*)  is  a 


local  minimum.  For  the  i1*1  iteration, 

1.  From  x*.  search  each  of  the  n  directions  d*  where  1-1,..., n. 

iV  '  /vj  J  '  9 

These  one-dimensional  searches  use  functional  values  to  locate 
a  minimum  in  each  direction.  A  quadratic  approximation  and  unimodal 
behavior  of  f(x)  is  assumed. 

2.  These  searches  locate  three  specific  points  in  R  —  x^,  the  last 
point;  xj,  an  expanded  point;  and  the  point  where  the  greatest 
decrease  in  f(^i)  occurred. 

3.  The  convergence  test  checks  ^  to  determine  if  f(x^)  is  a  local 
minimum.  If  x*  passes  the  convergence  condition  |x*  -  x^”1 j  <  &  , 

i 

then  x  -  x  .  If  not,  the  algorithm  continues. 

4.  The  modification  test  checks  the  decrease  in  f(x)  from  xj  to  x^* 


i+1 

These  functional  changes  specify  the  set  of  directions  for 

the  next  iteration.  The  same  mutually  conjugate  directions  may  be 
used  again  or  a  new  set  of  mutually  conjugate  directions  may  be 
determined* 


A  Conjugate  Gradient  Method..  This  section  presents  an  algorithm 

for  a  conjugate  gradient  with  restarts  method  (Ref  9  and  24),  Again, 

for  functions  which  are  not  quadratic,  the  method  will  not  locate  the 

i  1+1 

local  minimum  in  n  steps.  Instead,  the  method  iterates  from  x  to  x 
until  |  Vt(*1+1)|  <S, 


+  t 
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Ihe  direction  is  mutually  conjugate  to  the  previous  i  -  1  search 
directions; it  is  determined  using  the  previous  direction  d1”  and 
V £(&*■)»  Ihe  starting  point  in  Rn  is  x°;  the  initial  search  direction 
is  the  negative  of  the  gradient, -V  f(x°). 

An  iterative  process  is  used  to  locate  £  such  that  f(x  )  is  a 


local  minimum.  For  the  i^  iteration, 

1.  Calculate  Vf^) 

2.  From  jjc*,  use  a  one-dimensional  minimization  search  in  the  direction  d^ 
to  determine  the  step  length  t.  and  to  subsequently  locate  the  point 
xi+1- 


+  t 


min  fCx1  +  t^1) 
fci 

3 .  Calculate  V  f(jc,i+1) 

4.  The  convergence  test  checks  ^i+1  to  determine  if  f(£i+1)  is  a  local 

minimum.  If  ||  <$,  then  £  -  x;i+1.  If  not,  the 

algorithm  continues. 


5.  Compute  ^ 


nv^1H)ii2 

II V  t(&1)  II  2 


6.  Determine  the  next  mutually  conjugate  search  direction. 

d1+1  -  -V  £<s1+1) 

This  algorithm  locates  the  minimum  of  a  quadratic  function  with  a 


<7 


symmetric,  positive  definite  Hessian  matrix  in  n  or  less  iterations. 
However,  for  functions  that  are  not  quadratic,  the  minimum  will  gene¬ 
rally  not  be  determined  in  n  steps.  After  the  .  steps,  n  mutually 
orthogonal  directions  have  been  searched.  xn  may  or  may  not  have 
converged  rapidly  towards  x  . 

For  functions  with  slow  rates  of  convergence,  because  of  nearly 
i  i+1 

parallel  d  and  d  ,  Fletcher  and  Reeves  developed  the  restart  pro¬ 
cedure  (Ref  9).  After  every  n  +  1  steps,  the  method  reverts  to  the 
direction  of  steepest  descent,  the  largest  negative  gradient,  for  the 
next  search  direction.  That  is,  following  iteration  i  ■  n  +  1,  which 
located  x1*^,  the  direction  dn+^  would  not  be  specified  as  in  step  6 
above,  but  rather  d0*^  *  "Thus  the  whole  procedure  is 

restarted  from  the  currant  discarding  all  previous  experience  that 
would  normally  be  transmit ted  in  the  calculation  of  d*.  The  process 
remains  quadratically  convergent  provided  such  restarts  ore  not  more 
frequent  than  every  n  steps"  (Ref  9). 

The  CEDE  models  developed  in  Chapter  II  ore  maximized  with  these 
two  techniques.  Powell's  method  of  conjugate  directions  maximizes  the 
CEP-Included  CEDF  model;  a  conjugate  gradient  with  restarts  method 
maximizes  the  CEP-Excluded  CEDF  model.  These  methods  require  the 
vector  of  the  2m  independent  variables  x  and  the  function  CEDF(x); 
the  conjugate  gradient  technique  also  requires  the  gradient  of 
CEDF(£),  The  computerization  of  the  algorithm  that  maximizes  the  two 
CEDF(x)  models  is  presented  in  Chapter  IV. 

Greenwood  developed  a  similar  version  of  the  CEDF(x)  (Ref  12). 

His  algorithm  uses  a  different,  yet  related  approach  to  determine 
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optimal  DGZ  locations.  His  function  G(x)  depends  on  the  total  expected 
target  value  undamaged*  The  2m  first  partial  derivatives  of  G(x)  are 
set  equal  to  0.  Then  these  2m  nonlinear  equations  are  solved  itera- 
tively  to  yield  a  £  such  that  G(x )  is  a  minimum.  NUCWAVE  uses  a 
modified  Greenwood  technique  to  determine  optimal  DGZ  locations 
(Ref  25:4-3).  It  optimizes  one  weapon  at  a  time.  Hence,  it  iteratively 
solves  2  nonlinear  equations  to  determine  (X^,Y^)  —  the  coordinates  of 
weapon  i.  Then  it  repeats  the  process  for  the  next  weapon. 


IV.  Computerization ,  Verification,  and  Validation 
of  the  CEDF  Maximization  Algorithm 

The  Complex  Expected  Damage  Function  (CEDF)  maximization  algorithm 
includes  the  CEDF  models,  CEP-Included  and  CEP-Excluded,  and  the  opti¬ 
misation  techniques,  Powell's  method  of  conjugate  directions  and  the 
conjugate  gradient  with  restarts  method.  The  algorithm  determines 
optimal  DGZ  locations  for  a  finite  number  of  nuclear  weapons  against 
installations  in  a  target  complex  by  maximizing  the  CEDF. 

Evaluation  of  the  algorithm  consisted  of  three  related  stages  — 
construction,  verification,  and  validation.  These  stages  formed  an 
iterative  process  that  was  necessary  to  develop  user  confidence  in  the 
capability  of  the  algorithm.  Construction  is  the  formulation  and 
computerisation  of  a  model.  The  computerization  of  the  CEDF  maximi¬ 
sation  algorithm  used  a  modular  approach.  Smaller  segments  of  the  CEDF 
modal  were  developed  to  accomplish  lower  level  procedures.  These  seg¬ 
ments  became  subprograms  in  the  final  computer  code.  Verification  of 
the  CEDF  maximisation  algorithm  used  example  problems  to  insure  that  the 
results  of  each  subprogram  were  correct.  Validation  measures  the  rela¬ 
tive  agreement  between  the  model  and  the  system  modeled  (Ref  26:215). 
Validation  of  the  CEDF  maximization  algorithm  was  a  comparison  of  the 
results  from  the  algorithm  with  the  results  from  a  current  DGZ  model. 

This  chapter  presents  the  evaluation  of  the  CEDF  maximization  algorithm 
with  respect  to  these  three  stages. 

Computerization 

A  flow  chart  of  the  CEDF  maximization  algorithm  is  shown  in  Figure  6. 
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Figure  6*  Flowchart  of  the  CEDE  Maximization  Algorithm 
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Figure  7.  CEDE  Maximization  Algorithm  Subroutine  Hierarchy 


The  symmetry  of  Figure  6  illustrates  several  characteristics  of  the 
algorithm.  The  blocks  above  the  dashed  line  correspond  to  the  proce¬ 
dures  that  use  weapon  and  installation  parameters  to  develop  the  CEDF(^) 
and  the  gradient  of  the  CEDF(^).  The  lower  blocks  correspond  to  the 
optimization  techniques  that  were  used  to  maximize  the  respective  CEDF. 
The  left  blocks  correspond  to  the  CEP-Excluded  model;  alternately,  the 
right  blocks  correspond  to  the  CEP-Included  model.  The  upper  three  and 
lower  three  blocks  are  common  to  both  CEDF  models.  Each  block  is  a 
smaller  segment  of  the  CEDF  maximization  algorithm. 

The  computer  code  of  the  CEDF  maximization  algorithn  was  written 
using  FORTRAN  V.  Appendix  D  contains  a  listing  of  the  code  and  a  glos¬ 
sary  of  the  FORTRAN  variables.  The  computer  code  includes  a  driver 
module,  seven  subroutines  and  two  functions.  The  hierarchy  of  the  algo¬ 
rithm's  subprograms  is  shown  in  Figure  7.  All  program  variables,  inclu¬ 
ding  weapon  and  installation  parameters,  that  are  used  in  more  than  one 
subprogram,  are  stored  in  six  named  common  blocks.  Only  the  DGZ  coor¬ 
dinate  vector  £  is  transferred  between  subprograms  by  the  subprograms' 
argument  lists. 
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The  driver  module  OPTMZ  controls  the  CEDF  maximization  algorithm. 
OPTMZ  calls  the  five  highest  level  subroutines.  The  primary  functions 
of  these  subroutines  are:  (1)  INITLZ  inputs  user-specified  weapon 
and  installation  parameters,  (2)  WRADS  calculates  additional  instal¬ 
lation  and  weapon  parameters,  (3)  ZXCGR  is  a  conjugate  gradient  with 
restarts  subroutine  that  maximizes  the  CEP-Excluded  CEDF,  (4)  FWMIN 
is  a  conjugate  directions  subroutine  that  maximizes  the  CEP-Included 
CEDF,  and  (5)  OUTDGZ  outputs  the  final  DGZ  coordinates.  ZXCGR  calls 
GFUNCT,  a  subroutine  that  calculates  the  value  of  the  CEDF(x)  and  the 
gradient  of  the  CEDF(g).  PWMIN  calls  FUNCT,  a  subroutine  that  calculates 
the  value  of  the  CEDF(x).  GFUNCT  and  FUNCT  call  PDAM,  a  function  that 
calculates  the  probability  of  achieving  a  specified  level  of  damage  to 
installation  J  from  weapon  i,  Pd,  ..  PDAM,  in  turn,  calls  PDR,  a  func- 
tion  that  calculates  the  distance  damage  function,  P^Cr).  More  complete 
descriptions  of  these  subroutines  and  functions  are  given  below.  WRADS, 
PDAM,  and  PDR  are  modifications  of  subprograms  from  Mathematical  Back¬ 
ground  and  Programming  Aids  for  the  Hnyslcal  Vulnerability  System  for 
Nuclear  Weapons  (Ref  6). 

INITLZ.  This  subroutine  has  four  primary  functions.  First,  it 
reads  user-specified  weapon  and  installation  parameters  from  an  external 
file,  INDATA.  Appendix  E  contains  samples  of  an  input  file  and  an  output 
data  file.  For  each  weapon  i,  the  user  specifies  a  yield,  a  hob,  and 
a  Pa;  the  user  may  specify  initial  DGZ  coordinates.  Also,  for  each 
installation  j,  the  user  specifies  coordinates,  a  VNTK  code,  and  a 
value.  The  user  inputs  the  weapon  and  installation  latitude  and  longi¬ 
tude  in  degrees-minutes-seconds  and  the  direction  from  either  the  prime 

meridian  or  the  equator.  Positive  coordinates  are  east  cf  the  prime 
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meridian  and  north  of  the  equator. 

Second,  INITLZ  assigns  the  initial  coordinates  of  the  m  weapons, 
prior  to  maximization.  The  user  has  three  options:  (l)  provide 
independent  estimates  of  the  weapon  coordinates;  (2)  let  INITLZ  assign 
the  coordinates  of  the  m  highest  valued  installations  to  be  the  coor¬ 
dinates  of  the  m  weapons  in  decreasing  order  of  yield,  that  is,  the 
largest  yield  weapon  is  initially  located  at  the  highest  valued  instal¬ 
lation  ;  or  (3)  let  INITLZ  assign  the  coordinates  of  the  m  hardest 
installations  to  be  the  coordinates  of  the  m  weapons  in  decreasing  order 
of  yield,  that  is,  the  largest  yield  weapon  Is  initially  located  at  the 
installation  with  the  largest  VN  number. 

Third,  INITLZ  transforms  all  weapon  and  installation  degrees- 
minutes- seconds  into  coordinates  measured  in  feet  and  relative  to  a 
conroon  origin  in  an  XY  coordinate  system.  The  CEDF  maximization  algo¬ 
rithm  assumes  a  flat  earth  model  to  locate  all  coordinates.  Each  minute 
of  latitude  equals  6080  feet.  However,  one  minute  of  longitude  equals 
6080  feet  only  at  the  equator.  Whan  the  latitude  is  not  the  equator, 
one  minute  of  longitude  is  less  than  6080  feet  because  of  the  merging 
of  the  moridian  lines.  The  scale  factor  is  the  cosine  of  the  latitude. 

Lastly,  INITLZ  specifies  accuracy  requirements  for  ZXGGK  and  PWMIN. 
These  subroutines  need  prespecified  values  to  test  for  the  convergence 
of  $  to  the  maximum  value  of  CEOF(jc^). 

WRAPS.  This  subroutine  calculates  additional  weapon  and  installation 
parameters  from  tlie  user-specified  parameters.  For  each  installation  j, 
it  determines  a  distance  damage  sigma  (<rdj)  amd  a  Qy  Also,  for  each 
weapon  i-installation  j  combination,  it  calculates  a  weapon  radius,  WR^  j. 

These  calculations  are  described  in  Appendix  A. 
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PDR.  This  function  calculates  the  distance  damage  function  P^r), 

Slie  probability  of  achieving  a  specified  level  of  damage  to  installation  j 
from  weapon  i  whan  the  distance  r  between  Installation  j  and  weapon  i 
is  known.  See  Eq  (S). 

PDAM.  This  function  calculates  the  probability  of  achieving  a 
specified  level  of  damage  to  installation  j  from  weapon  i,  Pd^  y  For 
the  CEP-Excluded  CEDF  modal,  Pd|  j  is  the  distance  damage  function 
from  PDR.  For  the  CEP-Included  CEDF  model,  the  distance  r  between 
installation  j  and  the  impact  point  is  unknown,  and  Fd^  j  is  calculated 
using  Gauss -Legendre  quadrature  and  the  distance  damage  function. 

See  Eq  (8). 

GFUNCTi  This  subroutine  calculates  the  CEDF(x)  and  the  gradient  of 
the  CEDF(g) .  Onf  function  and  gradient  evaluation  requires  ra*n  calls 
to  function  PDAM.  These  cells  specify  Pd^  j  for  each  weapon  i-instal- 
lation  J  combination  using  the  CEP-Excluded  model.  CEDF(^c)  is  calculated 
using  Eq  (1).  Each  element  of  the  gradient  is  calculated  using  Eq  (11). 

FUNCT.  This  subroutine  also  calculates  the  CEDF(x)  using  Eq  (1), 
end  one  function  evaluation  requires  m*n  calls  to  function  PDAM.  How¬ 
ever,  these  calls  specify  Pd.  .  for  each  weapon  i-installation  j  combi- 

J 

nation  using  the  CEP-Included  CEDF  model. 

ZXQGR.  This  subroutine  minimizes  -CEDF(g)  for  the  CEP-Excluded 
modal.  It  ie  a  conjugate  gradient  with  restarts  routine  from  the 
International  Mathematical  and  Statistical  Libraries,  Inc.  (Ref  16:ZXGGR), 
ZXQGR  requires  function  and  gradient  evaluations  from  GFUNCT  and  the  DGZ 
coordinate  vector  g.  It  uses  two  control  parameters  —  DFPRED  and 
ACC.  DFPRED  specifies  an  estimate  of  the  expected  Increase  in  the  CEDF ; 

ACC  specifies  the  desired  accuracy  of  the  convergence  check.  This  check 
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requires  the  sum  of  the.  squares  of  the  gradient  elements  to  be  less 
than  ACC. 

RJMIN.  This  subroutine  minimizes  -CEDF(x)  for  the  CEP-Included 
model.  It  is  a  conjugate  directions  routine  from  Optimization  Tech¬ 
niques  with  FORTRAN  (Ref  17 i 331-343).  PWMIN  requires  function  eval¬ 
uations  from  FUNCT  and  the  DGZ  coordinate  vector  g.  It  also  uses  two 
control  parameters  —  ESCALE  and  E.  ESCALE  specifies  the  maximum  step 
size  multiplier  for  a  single  step  of  any  xkj  E  specifies  the  accuracy 
of  the  convergence  check.  This  check  requires  the  absolute  value  of 
the  differences  between  each  element  of  g  and  &  to  be  less  than  E. 

OUTDGZ.  This  subroutine  translates  the  XY  coordinates  of  the  final 
DGZs  from  feet  into  degrees -minu tea - s aconds  and  the  direction  from  either 
the  prime  meridian  or  the  equator.  Then  it  outputs  these  coordinates  to 
the  external  data  file,  TAPE 6. 

The  CEDF  maximization  algorithm  provides  three  sets  of  optimal 
DGZ  locations.  The  first  set  is  from  the  CEP-Excluded  CEDF  model  and 
ZXCGR  maximization;  the  second  set  is  from  the  CEP-Included  CEDF  model 
and  PWMIN  maximization.  The  last  sot  of  DGZ  locations  is  from  both 
CEDF  models  and  ZXOGR  and  PWMIN  maximization  —  a  mixed  technique. 

Verification 

The  verification  of  the  CEDF  maximization  algorithm  included  four 
phases.  Each  phase  verified  the  subprograms  of  the  algorithm 
using  example  problems.  For  each  computer  program,  the  results  of  each 
example  problem,  including  the  values  of  intermediate  variables,  were 
calculated  independently  of  the  respective  computer  program.  Then  the 
computer  program  solved  the  example  problem.  PRINT  statements  in  the 
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program  printed  values  of  most  FORTRAN  variables.  These  values  were 
compared  with  the  values  calculated  by  pencil  and  paper  to  verify  that 
the  computer  program  calculated  the  correct  values. 

The  first  phase  verified  five  lower  level  subprograms  —  PDR,  PDAM, 
WRADS,  INITLZ,  and  OUTDGZ.  Each  of  these  modules  was  coded  and  debugged 
as  a  small  FORTRAN  program.  Mathematical  Background  and  Programming 
Aids  for  the  Fhysical  Vulnerability  System  for  Nuclear  Weapons  includes 
example  problems.  Fifteen  of  these  problems  were  used  to  check  PDR, 

PDAM,  and  WRADS.  These  programs  calculated  the  same  values  as  the  exam¬ 
ple  problems.  The  outputs  of  1N1TLZ  for  several  example  problems  were 
compared  with  results  that  were  calculated  independently  of  the  computer 
program.  These  comparisons  indicated  INITLZ  was  properly  forming  the 
XY  coordinate  system  and  the  installation  and  DGZ  coordinate  vectors. 
Similarly,  the  outputs  of  OUTDGZ  for  several  test  cases  were  compared 
with  pencil  and  paper  calculated  results.  These  comparisons  indicated 
OUTDGZ  was  correctly  translating  the  final  DGZ  coordinate  vector  from 
feet  into  degreea-minutes-seconds  and  the  direction  from  either  the 
prime  meridian  or  the  equator.  These  five  subprograms  were  merged  into 
one  program  and  became  the  foundation  of  the  next  verification  phase. 

The  second  phase  verified  the  subroutine  FUNCT.  The  small 
programs,  PDAM  and  FDR,  became  FORTRAN  functions;  the  programs  WRADS, 
INITLZ,  and  OUTDGZ,  became  FORTRAN  subroutines.  FUNCT  calculates  the 
value  of  the  CEDF(^)  using  Eq  (1).  The  pencil  and  paper  calculated 
results  from  several  example  problems  were  compered  with  the  results 
from  FUNCT.  One  example  included  two  identical  installations.  Each 
installation’s  value  and  VNTK  code  were  15.0  and  15P2.  The  distance 
between  the  two  installations  was  6000  feet.  Ttoo  identical  weapons 
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Figura  8.  The  CEDF  for  a  two  weapon- two  installation  geometry. 

were  collocated  naifway  between  the  two  installations.  Each  weapon's 
yield,  height  of  burst,  and  CEP  were  100-kt,  1000  feet,  and  600  feet. 
The  independent  calculation  of  the  CEDF  for  this  complex  was  29.15;  the 
CEDF  value  from  FUNCT  was  29 . 1492 .  These  example  problems  indicated 
FUNCT  was  properly  calculating  the  CEDF(x). 

This  two  weapon-two  installation  complex  was  used  to  investigate 
the  results  of  moving  the  two  collocatad  DGZs.  The  geometry  of  this 
complex  is  shown  in  Figure  8.  Initially,  the  two  100-kt  weapons  were 
collocated  at  point  0,  and  the  CEDF  value  was  29.15.  The  CEDF  value 


decreased  as  the  two  DGZs  were  moved  in  opposite  directions  a  distance  d 
from  point  0.  A  graph  of  CEDF  versus  d  for  this  problem  is  also  included 
in  Figure  8.  This  example  indicated  the  existence  of  a  CEDF  maximum,  on 
the  line  considered,  as  the  distances  between  the  DGZs  and  the  instal¬ 
lations  varied. 

The  third  phase  verified  the  gradient  calculation  of  the  CEP- 
Excluded  CEDF  ;Ln  the  subroutine  GFUNCT.  GFUNCT  also  calculates  the 
CEDF(g)  using  Eq  (1).  The  gradient  vector  from  GFUNCT  was  checked  using 
two  weapon-installation  geometries.  Appendix  F  includes  tire  table  and 
calculations  used  to  verify  the  gradient  of  the  CEDF(x)  for  these  two 
examples . 

The  first  example  included  one  weapon  and  two  installations.  The 

first  installation's  value  and  VNTK  code  were  5000  and  11P2;  the  second 

installation's  value  and  VNTK  code  were  12000  and  15P2.  The  weapon's 

yield  and  height  of  burst  were  100-kt  and  1000  feet,  and  the  CEP  was 

0  feet.  Forty  values  of  the  CEDF(x)  were  calculated  for  different  DGZ 

locations.  The  x  direction  was  along  the  line  connecting  the  two 

installations.  These  40  values  were  then  plotted.  Figure  9  is  a  plot 

of  CEDF(x)  versus  x  for  this  example.  A  DGZ  between  the  two  installations 

was  selected  (x  *  63500)  and  the  gradient  calculated  using  two  methods. 

In  this  example,  the  gradient  liad  only  one  element  because  the 

y  variable  was  constant, and  only  the  x  variable  was  allowed  to  vary. 

The  gradient  values  for  the  two  calculation  methods  were  compared  with 

the  gradient  from  GFUNCT.  The  first  method  usad  a  difference  equation 
A  CEDF 

'  /5x  t0  aPProx^mate  fc^e  gradient.  For  the  DGZ  selected,  the  difference 
equation  approximation  of  the  gradient  was  3.939.  The  second  method  was 
pencil  and  paper  calculations  of  all  the  steps  necessary  to  determine 
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Figure  9.  The  CEDF  for  a  one  weapon- two  installation  geometry. 

the  gradient.  Chapter  II  presented  these  steps.  For  the  DGZ  selected, 
the  pencil  and  paper  calculation  of  the  gradient  was  3.9791.  Hie  value 
of  the  gradient  from  GFUNCT  for  the  DGZ  selected  was  3.9791.  These 
comparisons  indicated  the  subprogram  GFUNCT  was  properly  calculating 
the  gradient  of  the  CFl)F(x). 

The  second  example  included  two  weapons  and  three  installations. 

The  gradient  of  the  CEDF(x)  in  this  example  has  2m  or  f c- ;r  elements. 

Only  one  element  was  completely  checked  by  pencil  and  paper  calculations. 

A  location  for  each  DGZ  within  the  three-installation  complex  was  selected. 


Figure  10.  The  one  weapon* two  installation  geometry. 


Then  was  calculated  to  be  0.192548;  the  subsequent  value 

from  GFUNCT  was  0.19255136.  These  comparisons  indicated  GFUNCT  was 
correctly  forming  the  gradient  of  the  CEDF(g). 

The  last  phase  verified  the  CEDF  maximization  algorithm's  ability 
to  locate  a  local  maximum  of  the  CEDF.  All  subprograms,  the  subroutines 
ZXCGR  and  PWMIN,  and  the  driver  module  OPTMZ  were  merged  into  one  pro* 
gram  —  the  CEDF  maximization  algorithm.  The  two  installation-one 
weapon  complex  described  above  to  verify  GFUNCT  was  also  used  to  verify 
the  algorithm.  Figure  10  presents  this  complex,  several  initial  starting 
points  for  the  DGZ,  and  the  mean  location  £*  for  the  local  maximum  of 
the  CEDF(fl).  In  this  simple  example,  the  local  maximum  is  also  the 
global  maximum.  The  graph  of  the  CEDF(x)  versus  x  in  Figure  9  indicated 
the  maximum  CEDF  value  was  approximately  15000  for  65000  <  x  <  66000. 

The  CEDF  maximization  algorithm  was  run  with  seven  different 
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initial  DGZ.  From  the  initial  locations  1,  2,  3,  6,  and  7,  the  algorithm 
converged  to  a  maximum  CEDF  value  and  an  optimal  DGZ.  The  algorithm  did 
not  move  the  DGZ  from  the  initial  locations  4  and  5.  For  these  locations, 
the  ZXCGR  and  PWMIN  convergence  criteria  were  satisfied,  and  the  CEDF 
value  was  4950.  The  algorithm  did  not  move  the  DGZ  from  these  locations 
because  there  were  no  indications  of  a  CEDF  increase.  Chapter  VI  ex¬ 
plains  this  result  in  more  detail. 

The  mean  optimal  DGZ  location  and  CEDF  value  were  calculated  for 
the  other  five  initial  DGZ  locations.  For  the  CEP-Excluded  model  using 
ZXOGR,  the.  mean  location  of  the  optimal  DGZ  was  (65319,11184).  Ihe  stan¬ 
dard  deviation  for  x  was  8  feet;  for  y  it  was  15  feet.  The  mean  value  of 
the  CEDFQc)  was  15019;  the  standard  deviation  was  4.3.  For  the  CEP- 
Included  model  using  PWMIN,  the  mean  location  of  the  optimal  DGZ  was 
(65300,11172).  The  standard  deviation  for  x  was  58  feet;  for  y  it  was 
98  feet.  The  mean  value  of  the  CEDF(jj)  was  15006;  the  standard  devi¬ 
ation  was  16.6.  PWMIN  is  a  slower  optimization  routine;  hence,  less 
restrictive  convergence  criteria  were  established  for  PWMIN.  This  could 
account  for  PWMIN* s  smaller  CEDF  value  and  larger  standard  deviations 
for  x  and.  y. 

The  results  from  both  CEDF  models  were  compared  to  the  values  from 
Figure  9.  These  comparisons  indicated  that  the  CEDF  maximization  algo¬ 
rithm  located  an  optimal  DGZ  location  by  maximizing  the  CEDF(x)  for  this 
simple  two  installation-one  weapon  complex.  More  detailed  complexes 
are  considered  in  the  next  section  and  the  next  two  chapters. 

Validation 


Validation  measures  the  relative  agreement  between  the  model  and 
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the  system  modeled.  It  is  not  possible  to  make  comparisons  between  the 
CEDF  maximization  algorithm  and  the  real  world.  Similarly,  validation 
is  not  a  yes  or  a  no  answer;  it  is  a  qualitative,  relative  measure. 

The  CEDF  maximization  algorithm  results  for  two  example  problems  were 
compared  with  the  results  from  NUCWAVE  (Ref  29).  NUCWAVE  is  a  one-sided 
nuclear  weapons  allocation  war  gaming  model.  It  optimizes  the  damage 
to  a  set  of  targets  using  a  preselected  set  of  weapons. 

Two  of  the  primary  differences  between  the  CEDF  maximization  algo¬ 
rithm  and  NUCWAVE  are  NUCWAVE' s  starting  solution  and  optimization 
technique.  It  optimizes  sequentially  by  starting  the  largest  yield 
weapon  at  the  highest  valued  installation.  It  optimizes  over  the  (X,Y) 
coordinates  of  this  weapon.  Then  it  stores  the  final  coordinates  of 
this  weapon,  calculates  the  damage  of  all  affected  installations,  and 
determines  the  remaining  values  for  all  installations.  Then  it  opti¬ 
mizes  the  next  largest  yield  weapon  by  starting  it  at  the  highest 
remaining  valued  installation.  NUCWAVE  continues  to  iterate  through 
the  entire  weapon  set  until  no  further  movement  of  a  DGZ  results  in  an 
increase  in  the  total  expected  target  value  damage.  Chapter  I  includes 
a  description  of  NUCWAVE  methodology.  For  these  comparisons,  the  CEDF 
maximization  algorithm  assigned  the  initial  DGZ  locations  to  the  highest 
valued  installations. 

The  first  problem  included  one  weapon  and  two  installations.  This 
complex  was  very  similar  to  the  complex  in  Figure  10.  The  first  instal¬ 
lation's  coordinates  were  46°03'15"  N  -  45°10'00"  E;  its  VNTK  code  and 
value  were  11P2  and  5000.  The  second  installation's  coordinates  were 
46°03'25"  N  -  45°11'20"  E;  its  VNTK  code  and  value  were  15P2  and  12000. 

The  weapon's  yield,  height  of  burst,  and  CEP  were  100-kt,  1000  feet, 

63 


and  0  feet.  The  coordinates  of  NUCWAVE's  optimal  DGZ  were  46°03,22"  N  - 
45°10'58"  E;  the  total  expected  target  value  damage  was  16579  or  97.52% 
of  the  complex  value,  The  coordinates  of  the  CEDF  maximization  algo¬ 
rithm’s  optimal  DGZ  were  46°03'21"  N  -  45°10’50”  E;  the  CEDF  was  16788 
or  98.75%  of  the  complex  value.  This  represents  a  difference  of  approx¬ 
imately  575  feet  and  an  increase  in  CEDF  value  of  approximately  1%. 

Than  the  final  DGZ  coordinates  from  NUCWAVE  were  used  as  starting  coor¬ 
dinates  for  the  CEDF  maximization  algorithm.  The  CEDF  value  at  these 
coordinates  was  16647.  The  final  coordinates  for  this  run  of  the  algo¬ 
rithm  were  also  46°03'21M  N  -  45°10'50M  E.  These  results  Indicate  that 
the  CEDF  maximization  algorithm  achieves  comparable  results  with  an 
existing  model,  NUCWAVE. 

The  second  problem  included  two  weapons  and  five  installations. 

The  installations'  VNTK  codes  ranged  from  14P3  to  20P3;  the  instal¬ 
lations'  values  ranged  from  3000  to  12000.  The  total  complex  value 
was  33000.  Both  weapons'  yield,  height  of  burst,  and  CEP  were  100-kt, 
1000  feet,  and  0  feet.  The  coordinates  of  NUCWAVE's  optimal  DGZs  were 
46°01'58"  N  -  45°09'55"  E  and  46°00’48"  N  -  45°09'42"  E;  the  total 
expected  target  value  damage  was  28730  or  87.06%  of  the  complex  value. 
The  coordinates  of  the  CEDF  maximization  algorithm's  optimal  DGZs  were 
46°01'58"  N  -  45°09*54"  E  and  46°00'45"  N  -  45°09'38''  E;  the  CEDF 
value  was  29543  or  89.52%  of  the  complex  value.  This  represents  a 
difference  of  approximately  70  feet  in  the  first  DGZ  and  415  feet  in 
the  second  DGZ  and  an  increase  in  CEDF  value  of  approximately  3%. 

Again,  the  final  DGZ  coordinates  from  NUCWAVE  were  used  as  starting 
coordinates  for  the  CEDF  maximization  algorithm.  The  CEDF  value  at 
these  coordinates  was  29018.  The  final  coordinates  from  this  run 
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were  also  46°01'58"  N  -  45°09,54"  E  and  46°00*45"  N  -  45°09'38"  E; 
the  CEDF  value  was  also  29543. 

The  comparisons  between  the  results  from  the  CEDF  maximization 
algorithm  and  title  results  from  NUCWAVK  for  the  two  examples  indicate 
that  the  algorithm  correctly  determines  the  same  local  maximum  as 
NUCWAVE.  The  results  from  these  two  examples  do  not  validate  the 
algorithm,  but  because  the  DGZ  locations  were  consistent  between  NUCWAVE 
and  the  CEDF  maximization  algorithm,  the  algorithm's  results  are  not 
invalid.  Ttiaee  results  provide  the  user  confidence  in  the  capability 
of  tha  CEDF  maximization  algorithm. 
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V.  CEDF  Maximization  Algorithm  Properties 


The  Complex  Expected  Damage  Function  (CEDF)  maximization  algorithm 
determines  optimal  DGZ  coordinates  for  multiple  nuclear  weapons  against 
installations  in  a  target  complex  by  maximizing  the  CEDF.  Ihe  previous 
chapter  presented  the  computerization  and  evaluation  of  the  algorithm. 

The  flowchart  in  Figure  6  summarizes  the  algorithm's  modules.  It  also 
presents  the  algorithm's  three  CEDF  maximization  techniques:  (1)  ZXCGR, 
a  conjugate  gradient  with  restarts  optimization  method  that  maximizes  the 
CEP-Excluded  CEDF  model;  (2)  PWMIN,  Powell's  method  of  conjugate  direc¬ 
tions  that  maximizes  the  CEP-Included  CEDF  model;  and  (3)  a  mixed  tech¬ 
nique  that  uses  both  CEDF  models. 

The  ZXCGR  and  the  mixed  techniques  each  consist  of  two  stages.  The 
first  stage  of  the  ZXCGR  technique  has  a  less  restrictive  convergence 
criteria  than  the  second  stage,  and  its  DGZ  coordinates  are  used  as  ini¬ 
tial  DGZ  coordinates  for  the  second  stage  of  the  ZXCGR  maximization 
algorithm.  These  optimal  DGZ  coordinates  from  the  first  stage  of  the 
ZXOGR  algorithm  also  are  used  as  the  initial  DGZ  coordinates  for  the 
second  stage  of  the  mixed  technique.  Mixed  maximization  has  an  initial 
ZXCGR  stage  and  then  a  PWMIN  stage.  For  brevity  and  completeness,  the 
following  nomenclature  will  be  used  throughout  this  report.  ZXGGR  con¬ 
jugate  gradient  maximization  will  be  referred  to  as  ZXM.  Powell's 
method  of  conjugate  directions  will  be  referred  to  as  PWM.  Finally, 
the  mixed  technique  will  be  referred  to  as  MXM. 

This  chapter  contains  two  sections.  Each  section  presents  charac¬ 
teristics  of  the  algorithm's  three  CEDF  maximization  techniques.  Dif¬ 
ferences  and  similarities  between  ZXM,  PWM,  and  MXM  are  discussed. 
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Similarly,  sane  capabilities  and  limitations  are  presented.  The  first 
section  discusses  convergence  criteria,  installation  value  scaling,  user 
guidelines  and  user  cautions.  Also,  it  presents  comparisons  of  CEDF 
results  for  different  convergence  values.  The  second  section  presents 
characteristics  and  optimization  results  for  specific,  geometrically 
symmetric,  two  and  four -installation  target  complexes.  It  also  discusses 
the  effects  of  symmetric  gradient  elements  on  ZXM  maximization. 

Convergence  Criteria 

The  optimal  DGZ  coordinates  and  CEDF  values  from  the  CEDF  maximi¬ 
zation  algorithms,  ZXCGR  using  the  CEP-Excluded  model  and  PWMIN  using 
the  CEP- Included  model,  were  sensitive  to  the  convergence  control 
parameters. 

The  ZXCGR  control  parameters  are  ACC  and  DFFRED.  ACC  specifies 
the  desired  accuracy  of  the  convergence  check.  This  check  requires  the 
norm  of  the  gradient  to  be  less  than  ACC.  The  norm  of  the  gradient, 
IVFII,  is  the  sum  of  the  squares  of  the  gradient  elements.  When  ZXCGR 
locates  a  point  x*  in  2m  space,  such  that  the  norm  of  the  gradient  is 
less  than  the  prespecified  value  of  ACC,  the  optimization  routine  stops 
Iterating.  DFPRED  is  an  estimate  of  the  expected  increase  in  the  CEDF. 
ZXCGR  uses  it  to  determine  the  size  of  the  initial  change  in  each  x. 

The  values  of  the  installations  affected  the  choice  of  values  for 

ACC  and  DFPRED.  This  is  because  the  installation  values  directly  scale 

the  CEDF  and  the  magnitude  of  the  gradient  elements.  Eqs  (1)  and  (11) 

in  Chapter  II  present  this  relationship.  Most  of  the  example  problems 

in  this  report  used  installation  values  between  0  and  10000  (104).  If 
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installation  values  are  within  this  range,  then  the  most  versatile 
parameter  values  for  second  stage  ZXM  maximization  of  two  optimal  DGZs 
are  ACC  ■  0.001  and  DFPRED  -  1000.  If  the  installation  values  are 
not  within  this  range,  then  a  heuristic  guideline  is  suggested  to  assist 
the  user  in  estimating  reasonable  parameter  values.  For  either  case, 
the  first  stage  value  of  ACC  ■  0.01  determined  acceptable  DGZ  coor¬ 
dinates  for  a  wide  range  of  installation  values  and  number  of  weapons. 

A  general  guideline  to  determine  ACC  for  two  weapons  depends  on  the 
highest  valued  installation  in  the  target  complex.  The  highest  value 
is  rounded  up  to  the  largest  power  of  the  base  10.  'Ihen  ACC  equals  10”  ^ 
times  this  adjusted  value.  For  example,  for  a  four-installation  complex 
with  installation  values  between  2500  and  7000,  the  adjusted  value  would 
be  the  7000  rounded  up  to  10000  (10^).  Then  ACC  would  equal  10"^.  This 
heuristic  implies  that  a  smaller  ACC  is  needed  for  complexes  with  over¬ 
all  lower  valued  installations.  Intuitively,  this  makes  sense  because 
the  scaling  effect  of  smaller  installation  values  decreases  the  CEDF  and 
the  magnitudes  of  the  gradient  elements.  When  more  than  two  weapons  are 
used,  a  larger, ACC  value  is  needed  to  account  for  the  additional  gradient 
elements . 

As  an  example,  Table  tl  presents  the  results  of  two  ZXM  optimiza¬ 
tions  for  a  two  weapon-four  installation  target  complex.  Weapon  and 
Installation  parameters,  except  the  installation  values,  were  the  sane  for 
both  optimizations.  In  the  original  problem,  the  most  valuable  instal¬ 
lation's  value  was  7000 ,  and  an  ACC  of  0.001  was  used.  The  CEDF  maximum 
value  was  15436.  In  the  1/10  value  scaled  problem,  the  most  valuable 
installation's  value  was  reduced  by  a  factor  of  10  to  700,  and  a  smaller 
ACC  of  0.0001  was  used.  The  CEDF  maximum  value  was  1543.  Each  of  the 
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TABLE  II 


CEDF  Comparison  between  an  Original  Problem 
and  a  1/10  Value  Scaled  Problem 


Parameters 

Original 

Problem 

Scaled 

Problem 

Target 

VNTK 

Value 

Value 

1 

16P2 

3500 

350 

2 

22P2 

2500 

250 

3 

21P4 

5000 

500 

4 

19Q3 

7000 

700 

ACC 

0.0010 

0.00010 

H  VFjj  at  convergence 

0.0004 

0.00002 

CEDF  at 

convergence 

15436 

1543 

ZXM  CEDF  maximizations  located  essentially  the  same  coordinates  for  both 
DGZe.  Comparing  the  two  optimization  results,  the  coordinates  of  the 
first  DGZ  were  within  one  foot  of  each  other,  and  the  coordinates  of  the 
second  DGZ  were  within  five  feet  of  each  other.  Hence,  the  original 
and  the  scaled  optimization  problems  located  the  same  DGZ  coordinates 
without  regard  for  the  magnitude  of  each  installation's  value. 

Another  heuristic  is  suggested  for  estimating  the  value  of  DFPRED. 
For  the  m  weapons,  sum  the  values  of  the  m  highest  valued  installations. 
If  the  m  weapons  were  assigned  to  these  m  highest  valued  installations, 
then  this  sum  would  be  an  approximate  value  for  the  CEDF.  Next  deter¬ 
mine  the  total  value  of  all  the  installations.  Then  substract  the  value 
sum  of  the  m  highest  installations  from  the  complex's  total  value.  This 

difference  is  the  maximum  possible  CEDF  increase.  An  estimate  for  DFPRED 
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TABLE  III 


Comparison  of  ACC  Convergence  Criteria 


Parameter 

- ( - 

ACC 

0.0100 

0.0010 

0.0005 

11 7  F 1  at 
convergence 

0.0095 

0.0006 

0.0004 

Number  of 

function 

evaluations 

14 

18 

19 

CEDF  at 
convergence 

15433. 

15434. 

15437. 

DGZ  2  final 
coordinates 

(44840 J 23075) 

(44851,23090) 

(44850,23089) 

is  one-half  of  this  difference.  Again,  this  guideline  implies  that  for 

I 

lower  installation  values  DFFRED,  the  estimated  increase  in  the  CEDF, 
should  be  smaller.  j 

Generally,  by  decreasing  ACC,  ZXOGR  can  determine  better  estimates 
of  the  CEDF  maximum  and  its  respective  optimal  DGZ  coordinates.  For  the 
original  two  weapon-four  installation  complex  of  Table  II,  ZXOGR  was  used 
to  compare  the  CEDF  maximum  value  and  \V  F||  for  three  values  of  ACC.  : 

DFPRED  equaled  1000  for  these  three  examples.  Table  III  presents  the 
results  and  the  DGZ  2  optimal  coordinates  for  these  examples.  Ihe 
results  of  these  ZXM  maximizations  indicated  that,  by  decreasing  ACC, 

ZXM  can  determine  a  better  estimate  of  a  CEDF  local  maximum.  That  is, 

ZXM  can  achieve  a  larger  CEDF  value  and  a  smaller  ll  V  F  11  .  Ihe  final 
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coordinates  for  DGZ  1  were  the  same  for  the  three  cases.  Only  DGZ  2 
coordinates  were  different;  Table  III  indicates  this  difference  was 
barely  noticeable.  Since  the  final  DGZ  2  coordinates  were  within  17 
feet  of  each  other,  a  less  restrictive  ACC  is  acceptable.  That  is,  a 

_7 

value  of  ACC  smaller  than  10  times  the  adjusted  highest  installation 
value  is  unnecessary. 

The  results  of  a  similar  experiment  using  the  same  two  weapon-four 
installation  target  complex  indicated  that  the  value  of  DFFRED  also  did 
not  significantly  affect  the  CEDF  maximum  value  or  the  optimal  DGZ 
coordinates.  Five  values  of  DFPRED,  100,  1000,  2500,  5000,  and  6000, 
were  compared  using  a  constant  ACC  of  0.001. 

Occasionally,  ZXCGR  will  not  converge  satisfactorily  and  locate  an 
optimal  point  in  2m  space.  The  IMS!  subroutine  will  return  an  IER  ■  129 
error  message.  This  message  indicates  that  the  subroutine  abandoned  a 
line  search;  this  was  probably  because  of  conflicting  Information.  The 
gradient  may  indicate  that  a  point  is  not  optima].;  that  is,  IV  F||  >  ACC. 
However,  each  additional  iteration  may  be  on  either  side  of  the  optimal 
point  and  the  algorithm  i3  unable  to  terminate  satisfactorily.  For  most 
of  the  occurrences  of  this  error  message,  the  point  located  by  the  sub¬ 
routine  actually  was  a  good  estimate  to  the  local  CEDF  maximum.  Three 
options  are  available  to  the  user  when  title  algorithm  terminates  with 
this  error  message.  First,  select  another  ACC  value  and  rerun  the  same 
problem.  Second,  select  another  DFPRED  value  and  rerun  the  problem. 
Third,  compare  the  DGZ  locations  and  the  CEDF  maximum  value  with  the 
results  of  FWM  and  MXM.  Again,  for  most  occurrences  of  this  message, 

the  third  option  indicated  that  the  point  located  was  a  good  estimate  of 
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the  CEDE  maximum  value  and  the  respective  optimal  DGZ  coordinates. 

Just  as  ZXCGR  maximization  results  depended  on  the  values  of  the 
convergence  control  parameters,  ACC  and  DFFRED,  PWMIN  maximization  re¬ 
sults  depended  on  the  values  of  E  and  ESCALE.  E  specifies  the  desired 
accuracy  of  the  convergence  check.  This  check  requires  the  absolute 
value  of  the  differences  between  each  element  of  x  for  iteration  i  and 

/v 

each  element  of  x  for  iteration  i  -  1  to  be  less  than  E.  When  PWIN 

A/ 

locates  a  point  x*  in  2m  space,  such  that  all  element  differences  are 
less  than  the  prespeclfied  value  of  E,  the  optimization  routine  stops 
iterating.  ESCALE  is  the  maxitmm  step  size  multiplier  for  a  one-dimen¬ 
sional  search.  PWMIN  will  not  increment  each  x  by  more  than  ESCALE/E. 

Hie  effect  of  E  and  ESCALE  on  the  maxi  ram  CEDF  value  and  the  opti¬ 
mal  DGZ  locations  was  not  as  evident  as  the  ZXCGR  convergence  control 
variables.  Accordingly,  an  indepth  sensitivity  analysis  of  these  para¬ 
meters  was  not  accomplished.  Preliminary  investigations  indicated  that 
ESCALE/E  values  of  10000/0.1  were  the  most  effective  in  maximizing  the 
CEP-Included  CEDF  model.  E  values  of  1,  5,  and  10  often  resulted  in 
computer  runs  that  exceeded  60  seconds  of  computer  processing  (CP) 
time.  These  incomplete  runs  generally  aborted  after  the  third  or  fourth 
FWM  iteration.  Also,  ESCALE  values  of  1000  and  5000  were  examined. 

The  most  promising  values  of  ESCALE/E  were  5000/0.1  and  10000/0.1. 
These  two  combinations  were  used  for  more  than  143  CEDF  maximization 
algorithm  evaluations  using  3,  4,  5,  and  7 -installation  target  complexes. 
Eighty  of  83  runs  (96%)  using  the  ESCALE/E  values  of  5000/0.1  converged 
to  a  solution;  similarly,  54  of  £0  runs  (90%)  using  the  ESCALE/E  values  of 
10000/0.. I  converged  to  a  solution,  lue  other  nine  runs  were  termi¬ 
nated  because  cf  excessive  CF  time.  Twenty-eight  CEDF  maximization  runs 

72 


> 


were  identical  except  that  14  used  the  ESCALE/E  values  of  5000/0.1  and 
14  used  the  values  of  10000/0.1.  Differences  between  the  two  parameter 
pairs  for  two  criteria,  CEDF  maximum  value  and  CP  time,  were  evaluated. 
The  results  of  a  sign  test  indicated  that  there  was  no  difference  be¬ 
tween  the  parameter  pairs. 

Ihe  four  convergence  control,  parameters  need  to  be  specified  prior 
to  a  CEDF  maximization  algorithm  run.  Ihe  subroutine  INITLZ  initializes 
the  ACC  value  for  the  first  stage  of  ZXM  to  0.01.  The  user  provides  the 
ACC  value  for  the  second  stage  of  ZXM  through  the  external  file,  INDATA. 
Similarly,  the  user  provides  the  DFFRED  value  through  INDATA.  Appendix  E 
discusses  the  necessary  input  procedures.  The  subroutine  INITLZ  also 
initializes  the  values  of  ESCALE/E  to  5000/0.1.  If  the  user  desires 
different  PWM  convergence  control  parameters,  then  only  two  lines  of  the 
code  need  to  be  changed. 

The  norm  of  the  gradient ,  II  ^7  F II  ,  will  be  used  as  a  relative  indi¬ 
cator  of  convergence  for  all  ZXM  maximizations.  For  a  two- weapon  complex 
with  a  maximum  installation  value  of  10000,  W^J  Fll  -  0.001  implies  that 
the  mean  value  for  each  of  the  four  gradient  elements  is  approximately 
0.015.  That  is,  a  change  of  1000  feet  in  any  of  the  4  spatial  direc¬ 
tions  would  change  the  CEDF  by  only  15  value  points. 

Symmetry  Characteristics 

Two  simple  target  complexes  were  investigated  to  characterize  the 

CEDF  models  and  their  respective  maximization  techniques.  Initial  DGZ 

locations  were  selected  to  emphasize  special  features  of  ZXM  and  PWM 

optimization  techniques.  The  examples  included  either  symmetric  target 
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Figure  11.  A  symmetric  two-installation  complex, 

complexes  or  a  CEDF  with  symmetric  gradient  elements, 

A  two  weapon-two  installation  complex  was  analyzed  to  determine  the 
consequences  of  symmetric  gradient  elements  on  the  ability  of  ZXM  to  lo¬ 
cate  a  CEDF  maximum  value  and  optimal  DGZ  coordinates.  Figure  11  pre¬ 
sents  the  complex  geometry,  weapon  and  installation  parameters,  initial 
DGZ  locations,  and  the  ZXM  direction  of  DGZ  movement  for  five  cases. 

Each  of  these  cases  had  either  a.  geometrically  symmetric  weapon-instal¬ 
lation  complex  or  a  CEDF  with  symmetric  gradient  elements.  Tne  and  X2 
directions  were  along  the  line  segment  connecting  the  two  installations; 

the  and  x^  directions  were  perpendicular  to  this  line  segment.  For 
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each  of  these  cases,  FWM  converged  to  a  maximum  CEDF  value  of  9900. 

This  was  each  weapon's  Pa  times  the  total  value  of  both  installations. 

PWM  also  separated  collocated  DGZs.  The  results  of  the  five  cases  pro¬ 
vided  further  insight  into  the  capability  of  the  ZXM  maximization 
technique. 

CASE  1.  The  initial  DGZ  locations  for  two  identical  weapons  were 
collocated  halfway  between  two  identical  targets.  ZXM  neither  separated 
nor  moved  the  two  weapons.  This  was  because  all  four  of  the  CEDF  gra¬ 
dient  elements  equaled  0.  The  x3  and  x^  gradient  elements  were  0 
because  all  y  values  were  equal;  the  x^  and  x3  elements  were  0  because 
weapons  1  and  2  were  halfway  between  the  installations.  That  is,  one 
installation's  contributions  to  the  x^  and  x^  gradient  elements  can¬ 
celled  the  other  installation's  contributions.  ZXM  made  the  initial  DGZ 
coordinates  the  optimal  DGZ  coordinates  .with  a  CEDF  value  of  3465. 

CASE  II.  The  collocated  identical  weapons  in  CASE  I  were  separated. 
Weapon  1  was  moved  one  minute  of  longitude  west  (approximately  70  feet). 
With  this  move  from  the  complex's  geometric  center,  the  x^  gradient  ele¬ 
ment  was  no  longer  0;  ZXM  separated  the  two  DGZs  and  moved  weapon  1  towards 

installation  1  and  weapon  2  towards  Installation  2.  ZXM  converged  to  a 

-9 

maximum  CEDF  value  of  9900  and  two  optimal  DGZs  with  a  II V  F II  »  10  . 

CASE  III.  The  two  identical  weapons  were  again  collocated.  How¬ 
ever,  the  initial  DGZ  location  was  neither  halfway  between  the  two  iden¬ 
tical  installations  nor  along  the  line  segment  connecting  the  instal¬ 
lations.  The  gradient  elements  were  all  non-zero.  The  x^  and  ele¬ 
ments  were  equal  and  positive;  the  x3  and  x^  elements  were  equal  and  also 
positive.  In  this  case,  ZXM  did  not  separate  the  two  weapons  because 

their  respective  gradient  elements  were  the  same;  however,  it  did  move 
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the  two  weapons  together  towards  installation  1.  That  is,  ZXM  kept  the 
two  weapons  collocated  ..with  a  CEDF  value  of  5737  and  |iV  F  ||»  7.7  x  10”^. 

Next,  a  variation  of  this  case  was  examined.  The  same  initial  DGZ 
location  was  used,  but  the  yield  of  weapon  2  was  reduced  from  100-kt  to 
95-kt.  This  yield  reduction  altered  the  gradient  elements.  The  and 
X2  elements  were  no  longer  equal;  similarly,  the  x^  and  gradient 
elements  were  no  longer  equal.  The  two  different  weapons  separated 

from  the  same  initial  DGZ  location.  ZXM  moved  weapon  1  towards  instal-  j 

1 

lation  2  and  weapon  2  towards  installation  1.  ZXM  converged  to  a  maxi- 

«  ! 

mum  CEDE  value  of  9900  and  two  optimal  DGZs  with  11  ^  FH  ■  1.2  x  10  . 

1 

CASE  IV.  The  two  identical  weapons  were  separated  along  a  line 
segment  that  was  perpendicular  to  the  line  segment  connecting  the  two 
identical  installations.  Figure  11  displays  this  geometry.  The  two 
weapons  were  each  equidistant  from  the  two  installations.  Again,  because 
of  the  symmetry  of  the  target  complex,  the  identical  weapons,  and  the 

! 

identical  installations,  the  gradient  elements  were  symmetric.  The.  x. 

JL 

and  x^  gradient  elements  were  0.  That  is,  each  installation's  contri¬ 
butions  to  the  x^  and  x^  gradient  elements  negated  each  other.  The  x^ 
and  x^  gradient  elements  were  equal  in  magnitude,  but  opposite  in  direc-  j 

tion;  the  x^  element  was  positive  and  the  element  was  negative. 

Because  the  Xj  arid  x,,  gradient  elements  were  equal  to  0,  this  restricted 
the  DGZ  movements  to  only  the  x^  and  x^  directions.  With  the  and  x^ 
gradient  elements  with  equal  magnitude  but  opposite  direction,  ZXM 
moved  the  DGZs  directly  towards  each  other  to  point  P  in  Figure  11. 

This  point  was  along  the  line  segment  connecting  the  two  installations 

and  halfway  between  the  installations.  ZXM  converged  to  a  local  maximum  i 

CEDF  value  of  5729  and  one  collocated  DGZ  with  !i  V  F||  -  0.00028. 
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CASE  V.  A  variation  of  CASE  IV  was  examined.  The  initial  DGZ 
location  of  weapon  1  was  moved  one  minute  of  longitude  east  (approx¬ 
imately  70  feet).  The  x^  gradient  element  was  no  longer  0;  ZXM  moved 
weapon  1  towards  installation  2  and  weapon  2  towards  installation  1. 

ZXM  converged  to  a  maximum  CEDE  value  of  9895  and  two  optimal  DGZs 
with  8  VFil  -  0.00085. 

The  algorithm  located  CEDF  local  maximums  for  this  two  weapon- two 
installation  complex.  However,  the  algorithm's  optimal  DGZ  coordinates 
were  not  exactly  the  coordinates  of  the  two  installations.  For  this 
complex,  the  optimal  DGZ  locations  would  be  one  weapon  on  each  instal¬ 
lation,  since  the  weapons  and  the  installations  were  identical.  Never¬ 
theless,  these  examples  demonstrate  two  important  features,  an  opera¬ 
tional  characteristic  and  a  limitation,  of  the  CEP-Excluded  CEDF  model 
using  ZXM  maximization.  Both  of  these  are  a  result  of  the  synxnetry  of 
the  CEDF  gradient  elements.  These  cases  indicated  that  there  are  two 
types  of  gradient  symmetry.  There  is  symmetry  from  weapons  at  symme¬ 
trical,  initial  DGZ  locations  (CASES  I  and  IV)  and  from  identical  wea¬ 
pons  at  the  same  initial  DGZ  location  (CASE  III).  The  first  type  of 
symmetry  is  an  operational  characteristic;  the  second  type,  collocation 
of  initial  DGZ  locations,  is  a  minor  limitation.  This  limitation  means 
that  the  ZXM  maximization  algorithm  cannot  use  collocated  initial  DGZ 
locations.  The  PMW  maximization  algorithm  did  not  have  this  limitation. 

A  unique  one  weapon-four  installation  target  complex  also  exhibited 
symmetric  gradient  properties.  Figure  12  presents  U.e  complex  geometry, 
weapon  and  installation  parameters,  and  CEDF  values  for  different  po¬ 
tential  DGZ  locations. 
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Figure  12.  A  symmetric  four-installation  complex 
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The  geometric  shape  of  the  complex  was  not  a  square,  but  rather,  a 
rectangle.  The  north-south  distance  between  the  installations  was  8511 
feet;  the  east-west  distance  was  8448  feet.  There  were  four  local  maxi- 
mums  or  potential  DGZ  locations  for  this  target  complex;  one  at  each 
midpoint  of  the  four  line  segments  of  the  rectangular  perimeter.  The 
CEDF  value  at  each  local  maximum  was  7776.  The  CEDF  value  for  a  DGZ 
located  at  the  center  of  the  target  complex  was  3390.  The  CEDF  value 
for  a  DGZ  located  at  one  of  the  installations  was  4990.  Lastly,  for  a 
DGZ  located  halfway  between  an  installation  and  the  complex  center  along 
one  of  the  complex's  two  diagonals  at  one  of  the  points  L  in  Figure  12, 
the  CEDF  value  was  5660. 

CASE  VI.  The  initial  DGZ  was  the  geometric  center  of  the  target 
complex.  Figure  12  also  presents  this  location  and  the  optimal  DGZ 
location  for  PWM.  At  the  initial  DGZ  location,  the  damage  expectancy 
(DE)  for  each  installation  was  less  than  0.17.  FWM  moved  the  DGZ  in 
the  +x  direction  and  converged  to  a  maximum  CEDF  value  of  7776  and  an 
optimal  DGZ  location  between  installations  3  and  4.  The  DE  for  these 
installations  from  this  optimal  DGZ  was  approximately  0.78;  the  DE  for 
installations  1  and  2  was  approximately  0.001.  ZXM  did  not  move  the  DGZ; 
the  two  gradient  elements  were  0  because  the  initial  DGZ  was  at  the  geo¬ 
metric  center  of  the  complex. 

Next,  a  variation  of  this  example  was  examined.  The  initial  DGZ 
was  moved  one  minute  of  longitude  west  and  one  minute  of  latitude  north 
(approximately  120  feet).  PWM  moved  the  DGZ  in  the  +x  direction  as 
before.  However,  the  gradient  elements  were  >io  longer  0  because  the 
Initial.  DGZ  was  not.  at  the  geometric  center  of  the  complex.  ZXM 
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converged  to  a  CEDF  maximum  value  of  7776  and  an  optimal  DGZ  location 
between  installations  1  and  3  with  11 V  F II  -  0.0008.  A  possible  expla¬ 
nation  as  to  why  ZXM  moved  to  this  optimal  DGZ  is  presented  in  the 
next  case. 

CASF.  VII.  Tie  initial  DGZ  location  was  installation  1.  Figure  12 
also  presents  this  location  and  the  optimal  DGZ  locations.  PWM  again 
movid  the  DGZ  in  the  +x  direction  to  a  CEDF  maximum  value  of  7776  and 
an  optimal  DGZ  location  between  installations  1  and  3.  ZXM  did  not 
move  the  DGZ  towards  the  closest  local  maximum  as  it  did  in  CASE  VI. 
Instead,  ZXM  moved  the  DGZ  to  a  CEDF  maximum  value  of  7643  at  an  optimal 
DGZ  location  between  installations  3  and  4  with  II  7  Fll  -  0.0017. 

Investigation  of  the  first  20  iterations  of  ZXM  for  this  complex 
provided  a  plausible  explanation  as  to  why  ZXM  converged  to  this  optimal 
DGZ  instead  of  the  closest  DGZ.  Ihe  geometry  of  the  complex  was  not  a 
square,  but  rather,  a  rectangle.  Hence  at  the  Initial  DGZ,  the  two 
gradient  elements  were  not  exactly  equal.  The  x-j  gradient  element  was 
0.0352.  This  was  larger  than  the  X2  gradient  element  which  was  0.0310. 
Hence,  the  first  iteration's  search  direction  was  above  the  diagonal  of 
the  complex  along  the  line  U-U’  in  Figure.  12.  Figure  13  shows  an 
approximate  curve  of  CEDF  values,  using  seven  known  points,  along  U-U’. 
ZXM  located  the  point  U*.  ZXM  next  searched  along  the  line  segment  V-V' 
through  the  point  U*  and  the  two  potential  local  maximums  Vj*  and  V^*. 
Again,  Figure  13  shews  an  approximate  curve  of  CEDF  values,  using  five 
known  points,  along  V-V’.  From  U*,  ZXM  located  the  optimal  DGZ 
In  summary,  ZXM  does  not  always  move  the  initial  DGZ  towards  the  closest 
local  maximum. 
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Figure  13.  CEDF  values  along  line  segments  U-U'  and  V-V'  for  the 
symmetric  four -installation  complex. 


These  seven  cases  demonstrated  two  important  ZXM  maximization 
characteristics.  These  characteristics  depend  on  symmetry  of  the  CEDF 
gradient  elements.  This  symmetry  is  a  result  of  either  the  geometrical 
symmetry  of  the  target  complex  or  the  collocation  of  two  or  more  similar 
weapon  types.  This  second  characteristic,  a  minor  limitation,  prohibits 

i 

the  CEDF  maximization  algorithm  from  using  collocated  initial  DGZs.  The 
next  chapter  provides  a  description  of  initial  DGZ  locations  for  more 
typical,  nonsymmetric  target  complexes. 
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VI.  Algorithm  Results  for  Different  Initial  DGZ  Locations 


This  chapter  presents  results  of  the  CEDF  maximization  algorithm 
using  different  initial  DGZ  conditions.  The  important  results  are  the 
maximum  CEDF  value  and  the  optimal  DGZ  coordinates.  The  three  algorithm 
maximization  techniques  are:  ZXM,  conjugate  gradient  optimization  of 
the  CEP-Excluded  model;  PWM,  conjugate  directions  optimization  of  the 
CEP-Included  model;  and  MXM,  a  mixed  technique. 

Four  initial  DGZ  conditions  were  evaluated  using  three  target 
complexes.  However,  all  four  conditions  were  not  matched  with  each  of 
the  complexes.  The  four  initial  DGZ  conditions  for  m  weapons  against  a 
target  complex  were:  (1)  locating  the  weapons  at  the  m  highest  valued 
(HV)  installations,  (2)  locating  the  weapons  at  the  m  hardest  instal¬ 
lations,  (3)  locating  the  weapons  at  the  complex's  centroid,  and 
(4)  locating  the  weapons  at  m  pseudo-random  points. 

Intuitively,  the  most  logical  initial  DGZ  condition  was  the  highest 
valued  Installations,  and  the  least  logical  condition  was  random  loca¬ 
tions.  The  HV  condition  was  a  greedy  condition;  it  started  with  the 
maximum  damage  on  the  m  most  valuable  installations  and  then  searched 
for  other  DGZs  that  provided  an  increase  in  the  CEDF  value.  The  random 
locations  condition  was  not  canpletely  evaluated.  Instead,  for  a  two- 
weapon  complex,  six  pairs  of  initial  DGZs  were  evaluated,  and  the  six 
CEDF  values  and  optimal  DGZs  were  compared  to  each  other.  The  coordinates 
of  one  of  the  initial  DGZs  in  each  pair  were  fixed  and  common  to  all 
pairs.  The  coordinates  of  the  other  initial  DGZ  were  changed  for  each 
of  the  six  pairs  and  the  respective  six  runs  of  the  CEDF  maximization 
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algorithm.  The  results  of  these  pseudo-random  initial  DGZ  locations 
provided  additional  insight  concerning  different  initial  DGZ  locations. 

The  three  target  complexes  included  three,  four,  and  seven  instal¬ 
lations.  Ihe  CEDF  maximization  algorithm  located  optimal  DGZs  for  one, 
two,  or  three  weapons  against  these  complexes.  However,  each  complex 
was  not  matched  with  each  of  these  number  of  weapons.  That  is,  the 
three-target  complex  was  only  evaluated  using  one  and  two  weapons,  not 
three.  The  highest  valued  installation  in  any  of  the  complexes  was 
9000.  Hence,  for  the  convergence  control  parameters,  the  algorithm 
generally  used  values  of  ACC  ■  0.001,  DFPRED  ■  1000,  E  »  0,1,  and 
ESCALE  -  5000. 

Three  conclusions  were  made  from  the  results  of  these  examples. 
First,  the  algorithm  requires  some  indication  of  a  potential  increase 
in  CEDF  value  in  order  to  move  a  DGZ.  Second,  there  is  a  difference 
between  the  optimal  DGZ  coordinates  from  the  CEP-Excluded  model  using 
ZXM  maximization  and  those  from  the  CEP- Included  CEDF  model  using  PWM 
maximization.  This  difference  depends  on  a  weapon's  CEP  and  the  CEDF 
model  and  not  on  the  optimization  technique.  Third,  the  initial  DGZ 
coordinates  that  the  algorithm  uses  can  affect  the  maximum  CEDF  value 
and  the  optimal  DGZ  coordinates.  Statistical  evidence  of  these  con¬ 
clusions  is  not  presented.  Rather,  specific  examples  are  presented  that 
indicate  the  conclusions  are  not  invalid. 

A  Three- Installation  Complex 

CEDF  maximization  algorithm  results  were  analyzed  for  one  and  two 

weapons  against  a  three-installation  complex  using  different  initial  DGZ 

conditions.  Figure  14  shows  the  geometry  and  specific  parameters  of  the 
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target  complex.  The  total  available  target  value  for  the  complex,  reduced 
by  each  weapon's  Pa,  was  10890.  Figure  14  also  shows  the  optimal  DGZ 
coordinates  for  the  highest  valued  (HV)  initial  DGZ  condition  using  ZXM 
and  PWM  maximization.  DGZ  2  was  installation  3.  The  ZXM  algorithm  con¬ 
verged  to  a  maximum  CEDF  value  of  9812  or  90%  of  the  complex  value  and  to 
optimal  DGZ  coordinates  with  l7Fl  ■  0.00036.  Similarly,  PWM  converged 
to  a  maximum  CEDF  value  of  9223  or  85%  of  the  complex  value.  The  damage 
expectancy  (DE)  for  installations  1,  2,  and  3  from  ZXM  maximization  were 
0.96,  0.53,  and  0.99.  The  algorithm  did  not  move  DGZ  2  from  installation  3 
and  moved  DGZ  1  from  installation  1  towards  installation  2.  However,  the 
two  algorithms  located  the  optimal  DGZ  1  coordinates  480  feet  apart. 

This  difference  was  less  than  the  CEP  of  600  feet  and  initially  appeared 
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insignificant.  However,  the  differences  between  ZXM  and  PWM  maximum 
CEDF  values  and  DGZ  1  coordinates  were  important;  these  differences  do 
not  indicate  ZXM  is  a  better  algorithm.  These  differences  depended  on 
the  CEDF  model  and  are  discussed  in  more  detail  in  the  next  subsection. 

The  CEDF  maximization  algorithm  also  converged  to  a  local  CEDF 
maximum  for  the  centroid  initial  DGZ  condition*  ZXM  converged  to  a 
CEDF  value  of  6932  and  to  optimal  DGZ  coordinates  with  l|VF||  -  0,00003; 
PWM  converged  to  a  CEDF  value  of  6910.  However,  these  optimal  DGZ 
locations  were  not  the  same  locations  as  determined  using  the  HV  ini¬ 
tial  DGZ  condition.  Instead,  ZXM  and  PWM  moved  the  DGZs  towards  instal¬ 
lations  1  and  2  until  the  DE  for  each  installation  were  greater  than  0.99. 
The  final  DE  for  installation  3,  the  second  most  valuable  installation, 
was  less  than  0.001.  The  total  available  target  value  for  Installations 
1  and  2, reduced  by  each  weapon's  Pa, was  6930;  this  was  the  same  CEDF  . 
value  as  determined  by  ZXM  maximization.  These  optimal  coordinates, 
which  were  different  and  less  valuable  than  the  HV  Initial  DGZ  condition's 
optimal  coordinates,  were  also  identified  by  three  pairs  of  the  pseudo¬ 
random  initial  DGZ  condition. 

Six  pairs  of  the  pseudo-random  initial  DGZ  condition  were  also 
evaluated.  Figure  15  shows  the  initial  DGZ  locations  and  the  respective 
optimal  DGZ  locations  for  ZXM  maximization  for  two  cases.  It  also  shows 
the  approximate  weapons  radius  (WR)  for  each  installation.  For  all  six 
pairs,  PWM  maximization,  using  the  CEP- Included  CEDF  model,  converged  to 
the  same  optimal  DGZ  locations  as  determined  by  the  highest  valued  ini¬ 
tial  DGZ  condition.  The  pictorial  results  in  Figure  15  are  for  ZXM 
maximization  using  the  CEP-Excluded  CEDF  model. 
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Case  I.  For  the  initial  DGZ  locations  —  1-1’,  2-2' ,  and  6-6'  — 

ZXM  converged  to  a  maximum  CEDE  value  of  9810  or  90%  of  the  complex  value 
and  to  optimal  DGZ  coordinates  with  iv?\\  <  0.004*  this  was  the  same 
CEDF  local  maximum  that  the  HV  initial  DGZ  condition  located. 

Case  II.  For  the  initial  DGZ  locations  —  3-3',  4-4',  and  5-5'  — 
ZXM  converged  to  a  maximum  CEDF  value  of  6932  or  64%  of  the  complex  value 
and  to  optimal  DGZ  coordinates  with  "9  F  I  <  0.0004.  This  was  the  same 
CEDF  local  maximum  that  the  centre  id  initial  DGZ  condition  located.  Each 
of  the  optimal  DGZ  locations  are  slightly  different  locations;  however, 
each  of  the  locations  are  equivalent.  ZXM  moved  from  3-3',  4-4* ,  and 
5-5'  towards  installations  1  and  2  until  the  DE  for  these  installations 
was  greater  than  0.99. 

A  possible  explanation  exists  for  the  difference  between  the  two 
CEDF  local  maximuns  for  the  two  cases.  The  Case  I  initial  DGZ  locations 
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each  had  one  initial  DGZ  (1* ,  2',  and  6')  within  the  WR  of  one  of  the 
two  highest  valued  installations.  Ihe  other  initial  DGZ  (1,  2,  and  6) 
was  outside  the  WR  of  all  installations.  Alternately,  the  Case  II 
initial  DGZ  locations  had  neither  of  the  initial  DGZs  within  the  WR  of 
the  two  highest  valued  installations.  Hence,  for  Case  II,  the  algorithm 
moved  one  initial  DGZ  location  (3',  A',  and  5')  towards  installation  2 
and  the  other  initial  DGZ  location  (3,  4,  and  5)  towards  installation  1, 
the  mout  valuable  one.  Figure  4  shows  that  the  probability  of  achieving 
a  specified  level  of  damage  to  an  installation  at  the  WR  is  less  than  0.5. 
Using  WR  to  interpret  CETF  local  maximums  is  not  a  definitive  technique. 
However,  the  relationship  between  the  location  of  an  initial  DGZ  and  an 
installation's  WR  does  provide  insight  and  a  possible  explanation  for 
the  two  CEDF  local  maximums. 

In  summary,  these  results  point  out  the  first  of  three  conclusions 
of  this  study.  The  CEDF  maximization  algorithm  requires  some  indication 
of  a  potential  increase  in  the  CEDF  in  order  to  move  a  DGZ.  That  is,  if 
there  is  no  indication  of  a  CEDF  increase  in  the  direction  of  a  valued 
installation  and  there  is  an  indication  of  a  CEDF  increase  in  the  direc¬ 
tion  of  a  lesser  valued  installation,  then  the  algorithm  may  move  the 
DGZ  towards  the  lesser  valued  installation.  Eventually,  the  algorithm 
will  converge  to  a  less  valuable  CEDF  local  maximum. 

Using  the  CEDF  maximization  algorithm  to  evaluate  one  weapon  against 
this  three- installation  complex  produced  similar  results  for  three  ini¬ 
tial  DGZ  conditions.  For  the  highest  valued  condition,  ZXM  started  from 
installation  1  and  converged  to  a  CEDF  maximum  value  and  an  optimal  DGZ 
with  &  VFl  ■»  0,00030,  The  coordinates  of  this  DGZ  coincided  with  the 
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coordinates  of  one  of  the  two-weapon  HV  condition  DGZs.  Similarly,  FWM 
converged  to  the  same  optimal  DGZ  location  as  one  of  the  two-weapon  HV 
condition  DGZs.  This  result  indicated  that  the  two  weapons  against  the 
three  installations  were  not  dependent  but  rather,  unrelated  DGZs.  The 
fact  that  DGZ  2  never  moved  from  installation  3  for  the  two-weapon  exam¬ 
ple  also  indicated  that  the  two  DGZs  were  independent. 

A  comparison  between  the  one-weapon  centroid  initial  DGZ  condition 
results  and  the  *‘  ,,o-weapon  HV  initial  DGZ  condition  results  indicated  the 
sensitivity  of  the  gradient.  Both  of  these  examples  located  an  optimal 
DGZ  at  10' 09"  E  -  03*27"  N.  However,  when  the  optimal  coordinates  in 
feet  were  compared,  the  two  DGZs  were  approximately  30  feet  apart.  For 
one  DGZ,  the  1|V  F  |1  -  0.00031;  for  the  other  DGZ,  only  30  feet  away, 

||  VFft  -  0.00760. 

For  the  hardest  initial  DGZ  condition,  neither  ZXM  nor  FWM  moved  the 
one  weapon  initial  DGZ.  The  DGZ  started  at  installation  3  and  remained 
there.  The  8  V  Ftt  •  10"*^  at  this  point.  The  CEDF  value  for  ZXM  was 
3962  or  99 7.  of  the  value  of  installation  3;  the  CEDF  value  for  PWM  was 
3918  or  98%  of  the  value  of  installation  3.  The  difference  in  these 
CEDF  values  depended  on  the  CEDF  model  and  are  discussed  next. 

The  CEP  Effect 

CEDF  maximization  algorithm  results  for  the  three- installation 
complex  using  different  initial  DGZ  conditions  point  out  the  second 
conclusion  of  this  study.  Both  a  weapon's  CEP  and  the  CEDF  model  used 
to  maximize  the  CEDF  affect  the  maximum  CEDF  value  and  the  optimal  DGZ 
coordinates.  Specific  results  from  three  previous  three  installation 
examples  provide  evidence  to  support  this  conclusion. 
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First,  the  results  from  the  one  weapon-three  installation  complex 
using  the  hardest  installation  initial  XZ  condition  highlight  this  dif¬ 
ference  between  the  two  CEDF  models.  The  CEP- Excluded  CEDF  model,  using 
ZXM  maximization,  converged  to  a  CEDF  value  of  3962;  the  CEP- Included 
CEDF  nodel,  using  PWM  maximization,  converged  to  a  CEDF  value  of  3918. 

The  optimal  DGZ  coordinates  for  these  algorithms  were  within  1  foot  of 
each  other.  The  difference  in  CEDF  values  was  attributed  to  the  CEDF 
models.  The  Pdj^ j  for  the  CEP-Excluded  model  does  not  include  weapon 
accuracy  or  CEP.  This  probability  is  the  distance  damage  function  value, 
P(j(r).  The  Pd^  j  for  the  CEP-Included  model  does  include  weapon  CEP. 
Hence,  this  probability  is  less  than  P^(r).  Consequently,  the  PWM 
damage  expectancy  for  an  installation  is  less  than  the  ZXM  damage  ex¬ 
pectancy  for  the  same  installation. 

The  second  example  that  supports  the  conclusion  was  the  two  weapon- 
three  installation  complex  using  the  highest  valued  initial  DGZ  condition. 
Analysis  of  this  example's  results  provided  an  explanation  for  the  CEDF 
differences  between  ZXM  and  PWM.  Figure  14  shows  the  optimal  DGZ  coordi¬ 
nates  for  these  algorithms.  Only  XZ  1  coordinates  are  considered; 

XZ  2  coordinates  were  the  same  for  both  algorithms.  PWM  converged  to 
optimal  XZ  1  coordinates  approximately  480  feet  closer  to  installation  1 
than  ZXM.  Two  additional  initial  XZ  conditions  were  necessary  to  fur¬ 
ther  investigate  this  difference.  For  the  initial  XZ  coordinates,  the 
first  condition  used  the  EVJM  optimal  coordinates;  the  second  condition 
used  the  ZXM  optimal  coordinates.  Table  IV  presents  the  final  coordinates 
and  CEDF  values  from  these  initial  XZ  conditions.  Ihe  coordinates  are 
in  feet  and,  even  though  they  appear  different  within  the  three  optimi¬ 
zation  categories,  they  are  not.  The  ZXM  coordinates  were 
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TABLE  IV 


Comparison  of  ZXM,  FWM,  and  MXM  Optimal  DGZs 


Initial 

DGZ 

Condition 

ZXM  CEDF  ZXM 

Values  Final 

Start  End  Location 

PWM  CEDF  FWM 

Values  Final 

Start  End  Location 

MXM  CEDF  MXM 

Values  Final 

Start  End  Location 

Highest 

Value 

8985  9812 

858,1039 

8972  9220 

585,642 

9009  9223 

664,627 

PWM 

Optimal 

9480  9813 

837,1055 

9222  9222 

643,629 

9008  9218 

660,593 

ZXM 

Optimal 

9812  9812 

824,  976 

9073  9223 

627,632 

9009  9223 

639,638 

Mean 

Values 

9812 

840,1023 

9222 

618,634 

9221 

654,619 

10'09"  E  -  1Q'27"  N  and  the  FWM  and  MXM  coordinates  were  10'06"  E  -  03' 23"  N. 

Using  FWM  optimal  DGZ  coordinates  as  the  initial  DGZ  coordinates,  the 
algorithm  produced  three  results.  First,  the  PWM  optimal  DGZ  location  was 
the  initial  DGZ  location.  Second,  the  FWM  optimal  DGZ  coordinates,  the 
initial  coordinates,  were  not  optimal  for  ZXM.  At  these  initial  DGZ  coor¬ 
dinates,  the  ZXM  CEDF  value  was  9480.  ZXM  maximization  moved  the  DGZ  from 
10' 06"  E  -  03' 23"  N  back  to  10' 09"  E  -  03' 27"  N  and  a  maximum  CEDF  value 
of  9813.  ZXM  optimal  coordinates  are  initial  DGZ  coordinates  for  the  mixed 
technique,  MXM.  Third,  these  initial  MXM  coordinates  were  not  optimal 
for  MXM.  At  these  DGZ  coordinates,  the  MXM  CEDF  was  9008.  MXM  maximi¬ 
zation  moved  the  DGZ  from  10' 09"  E  -  03' 27"  N  back  to  10' 06"  E  -  03*23"  N, 
the  FWM  optimal  DGZ  coordinates,  and  a  maximum  CEDF  value  of  9218. 

The  algorithm  produced  three  similar  results  when  it  used  the  ZXM 
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optimal  DGZ  coordinates  as  the  initial  DGZ  coordinates.  First,  the  ZXM 
maximization  optimal  DGZ  location  was  the  initial  DGZ  location.  Second, 
the  ZXM  optimal  DGZ  coordinates,  the  initial  coordinates,  were  not  opti¬ 
mal  for  PWM.  At  these  initial  DGZ  coordinates,  the  IWM  CEDF  value  was 
9073.  PWM  maximization  moved  the  DGZ  from  this  initial  DGZ  back  to 
10' 06"  E  -  03' 23"  N  and  a  maximum  CEDF  value  of  9223.  Third,  the  ini¬ 
tial  coordinates,  10' 09"  E  -  03' 27"  N,  again  were  not  optimal  for  MXM. 

At  these  coordinates,  the  MXM  CEDF  value  was  9009.  MXM  maximization 
moved  the  DGZ  from  this  initial  DGZ  location  back  to  10*06"  E  -  03' 23"  N 
and  a  maximum  CEDF  value  of  9223. 

Finally,  analysis  of  a  third  three- installation  complex  provided 
further  insight  into  the  capability  of  the  CEDF  models.  In  the  previous 
example,  the  mixed  technique  moved  the  DGZs  from  the  ZXM  optimal  DGZ 
coordinates  to  the  PWM  optimal  DGZ  coordinates.  However,  this  readjust¬ 
ment  did  not  occur  in  all  examples.  For  instance,  ZXM  maximization  for 
the  pseudo-random  initial  DGZ  pairs,  3-3',  4-4',  and  5-5',  converged  to 
optimal  coordinates  that  were  different  from  the  PWM  optimal  coordinates. 
The  PWM  optimal  coordinates  were  the  HV  coordinates;  the  ZXM  optimal 
locations  were  near  installations  1  and  2.  The  mixed  technique  was  un¬ 
able  to  move  the  DGZs  from  these  ZXM  optimal  coordinates  to  the  PWM 
optimal  coordinates.  The  11  ^7  Fll  <  0.004  for  each  ZXM  local  maximum. 

These  CEDF  maximization  algorithm  results  indicated  that  each  CEDF 
modal  located  a  unique  set  of  optimal  DGZs.  This  occurred  because  of 
the  difference  in  Pd^  j  for  the  two  models.  Pd^ j  is  larger  for  the 
CEP-Excluded  model  than  it  is  for  the  CEP-Included  model. 
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Figure  16.  Multiple  local  optimal  DGZs  for  a  two  weapon  four  Installation  complex. 

Larger  Complexes 

For  a  four-installation  complex  using  the  four  initial  DGZ  conditions, 
the  CEDF  maximization  algorithm  produced  results  similar  to  the  results 
for  the  three-installation  complex.  Figure  16  shows  the  ge< me try  and 
specific  parameters  of  the  target  complex.  The  total  available  target 
value  for  the  complex,  reduced  by  each  weapon's  Pa,  was  17820.  Figure  16 
also  shows  the  optimal  DGZ  locations  for  several  initial  DGZ  conditions. 

Four  local  CEDF  maximum*  and  their  respective  optimal  DGZ  pairs  were 
located  for  this  complex j  the  HV  (highest  valued)  DGZs,  1-1',  2-2', 
and  3-3'. 

Again,  ZXM  and  FWM  located  their  highest  valued  DGZ  1  at  slightly 
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different  coordinates.  The  difference  between  the  two  locations  was  90 
feet.  ZXM  converged  to  a  CEDF  maximum  value  of  15434  or  87%  of  the  com¬ 
plex  value  and  to  optimal  DGZ  coordinates  with  ftVFll  -  0.00056.  A a 
with  the  three -installation  complex,  the  ZXM  optimal  DGZ  coordinates 
were  not  optimal  for  PWM.  The  mixed  technique  moved  DGZ  1  from  the  ZXM 
optimal  DGZ  to  the  PWM  optimal  DGZ  and  a  CEDF  maximum  value  of  15142  or 
85%  of  the  complex  value. 

The  three  remaining  local  CEDF  maxima  appear  to  be  related.  When 
the  initial  DGZ  condition  was  two  weapons  at  the  complex's  centroid, 

ZXM  converged  to  a  local  CEDF  maximum  value  cf  10921  and  to  the  optimal 
DGZ  pair  3-3'  with  ftV  Ffc  -  0.0007.  Similarly,  when  the  initial  DGZ 
condition  was  the  two  hardest  installations,  PWM  converged  to  a  local 
CEDF  maximum  value  of  10798  and  to  the  optimal  DGZ  pair  2-2 ' .  ZXM  and 
the  same  initial  DGZ  condition  produced  a  third  local  CEDF  maximum. 

When  the  initial  DGZ  condition  was  the  hardest  installations,  ZXM 
converged  to  a  CEDF  maximum  value  of  11431  and  to  the  optimal  DGZ  pair 
1-1'  with  llVFtt"  0.0027.  This  optimal  DGZ  pair  had  a  larger  CEDF 
value  than  the  pairs  2-2'  and  3-3'.  Yet,  it  had  a  smaller  CEDF  value 
than  the  pair  of  highest  valued  DGZs. 

The  local  CEDF  maximum  for  PWM,  when  the  initial  DGZ  condition  was 
two  weapons  at  the  complex's  centroid,  was  the  same  local  maximum  as 
determined  using  the  highest  valued  initial  condition.  Additionally, 
this  local  maximum  was  located  by  all  eight  of  the  pseudo-random  initial 
DGZ  conditions. 

Results  of  the  CEDF  maximization  algorithm  using  only  one  weapon 
against  this  four-installation  complex  were  examined.  ZXM  and  PWM, 
using  the  centroid  initial  DGZ  condition,  converged  to  the  highest 
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valued  DGZ  1  in  Figure  16.  This  is  one  of  the  two  op'imal  DGZ  loca¬ 
tions  determined  by  the  two  weapon  evaluation.  ZXM  converged  to  a  CEDF 
maximum  value  of  8547  or  78%  of  the  value  of  installations  1,  2,  and  3 
and  to  an  optimal  DGZ  location  with  |W  F|l  ■  0.0015.  FWM  achieved  a 
CEDF  maximum  value  of  8252.  Next,  the  algorithms  used  the  highest 
valued  initial  DGZ  condition  for  one  weapon.  Neither  algorithm  moved 
the  optimal  DGZ  from  the  initial  DGZ,  installation  4.  ZXM  and  PWM 
terminated  with  a  CEDF  maximun  value  of  6930. 

The  CEDF  maximization  algorithm's  results  for  the  three  and  four- 
installation  complexes  point  out  the  last  conclusion  of  this  study.  The 
most  likely  to  succeed  initial  DGZ  condition  is  to  use  the  coordinates 
of  the  m  highest  valued  installations  as  the  initial  DGZ  coordinates. 

For  all  examples  considered,  the  other  three  initial  DGZ  conditions 
located  at  least  one  local  CEDF  maximum  that  was  less  valuable  than  the 
local  CEDF  maximum  determined  from  the  highest  valued  initial  DGZ 
condition.  However,  there  is  always  an  exception.  The  CEDF  results  of 
the  simple  one  weapon-four  installation  complex  indicated  the  HV  initial 
DGZ  is  not  always  the  best.  For  this  reason,  the  CEDF  maximization  algo¬ 
rithm  does  not  include  a  decision  structure  to  determine  the  initial  DGZ 
condition  to  use.  Sometimes,  one  condition  may  be  more  likely  than 
another  to  succeed  and  to  achieve  the  most  valuable  local  CEDF  maximum. 

The  results  of  a  three  weapon-seven  installation  complex  were  ana¬ 
lyzed  to  further  define  the  CEP  effects  of  the  two  CEDF  models.  Only 
three  CEDF  maximization  algorithm  runs  were  made  with  this  complex.  For 
the  three  runs,  all  weapon  and  installation  parameters  remained  constant 
except  each  weapon's  CEP,  and  the  algorithm  used  the  highest  valued 
initial  DGZ  condition.  Using  each  weapon's  CEP  -  0  feet,  the  algorithm 
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converged  to  a  CEDF  maximum  value  and  optimal  XZ  coordinates  for  the 
three  weapons.  Each  weapon's  CEP  equaled  250  feet  for  the  algorithm's 
second  run.  For  this  example,  the  ZXM  optimal  XZ  coordinates  remained 
the  same,  as  they  should  have.  The  FWM  optimal  XZs  were  along  line 
segments  between  the  highest  valued  initial  XZs  and  the  optimal  XZs 
from  the  first  run  when  each  CEP  was  0  feet.  However,  because  each 
weapon's  CEP  >  0,  each  of  the  second  run  optimal  XZs  were  slightly 
closer  to  their  respective  highest  valued  initial  XZ.  Similarly,  each 
weapon's  CEP  was  400  feet  for  the  algorithm's  final  run.  Again,  the 
FWM  optimal  XZs  were  along  the  same  line  segments  as  the  optimal  XZs 
of  examples  one  and  two.  These  optimal  XZs  were  even  closer  to  their 
respective  highest  valued  initial  XZ.  Thus,  the  effects  of  each  wea¬ 
pon's  CEP  need  to  be  included  in  locating  optimal  XZ  coordinates  in  a 
target  complex. 
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VII.  Conclusions  and  Recommendations 


The  primary  objective  of  this  study  was  to  investigate  optimal  DGZ 
locations  within  a  target  complex.  In  order  to  accomplish  this,  it  was 
necessary  to  develop  the  Complex  Expected  Damage  Function  (CEDF)  maxi¬ 
mization  algorithm.  Ihe  algorithm  locates  optimal  DGZ  coordinates  for 
multiple  nuclear  weapons  against  installations  in  a  target  complex.  It 
does  this  by  maximizing  the  expected  target  value  damage  for  all  instal¬ 
lations.  Ihe  two  subobjectives  of  tills  study  were  to  determine  the  sen¬ 
sitivity  of  the  algorithm's  results,  the  maximum  CEDF  value  and  the  opti¬ 
mal  DGZ  coordinates,  to  two  factors:  first,  the  mathematical  technique 
used  to  locate  the  optimal  DGZs;  second,  the  initial  DGZ  locations  prior 
to  CEDF  maximization.  Ibis  chapter  discusses  these  objectives  and  their 
associated  conclusions. 

The  CEDF  maximization  algorithm  contains  two  related  algorithms, 
and  both  of  these  include  two  elements.  Ihe  first  element  is  the  CEDF, 
a  mathematical  model  of  the  total  complex  expected  target  value  damage. 
The  CEDE  is  an  unconstrained,  nonlinear  function  of  2m  variables  —  the 
(Xi,Yi)  coordinates  for  each  of  the  m  weapons.  There  is  a  CEDF  model 

for  each  of  the  related  algorithms.  The  basic  element  of  each  model  is 
Pd4  .  —  the  probability  of  achieving  a  specified  level  of  damage  to 
installation  j  from  weapon  i.  This  study  two  forms  of  the  Pd^  j 
function;  both  forms  depend  on  the  DIA  Fhysical  Vulnerability  system. 

The  major  difference  between  the  two  CEDF  models  is  their  respective 

Pd.  .  expressions.  The  CEP-Excluded  CEDF  model  assumes  each  weapon's 
if  J 

CEP  is  0  feet.  This  simplifying  assumption  results  in  two  conditions: 
a  less  complicated  mathematical  expression  for  the  CEDF  and  a  closed- 
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form  analytical  expression  for  the  gradient  of  the  CEP-Excluded  CEDF. 

The  CEP-Included  CEDF  model  includes  each  weapon's  CEP;  it  is  a  more 
complicated  expression. 

The  second  element  of  the  algorithm  is  a  nonlinear  optimization 
technique  that  maximizes  the  CEDF  models  and  locates  the  corresponding 
optimal  DGZs.  Since  an  analytical  expression  for  the  gradient  of  the 
CEP-Included  CEDF  model  was  not  available,  the  algorithm  maximizes  this 
CEDF  using  a  direct  search  technique  —  Powell's  method  of  conjugate 
directions,  PWM.  This  numerical  search  technique  requires  only  func¬ 
tion  evaluations  to  locate  a  local  maximum.  Conversely,  an  analytical 
expression  for  the  gradient  of  the  CEP-Excluded  model  was  calculated. 

The  algorithm  maximizes  this  CEDF  using  a  gradient  search  technique  — 
a  conjugate  gradient  with  restarts  method,  ZXM.  the  algorithm  also 
contains  a  third  CEDF  maximization  technique. —  a  mixed  technique,  MXM. 
This  technique  consists  of  two  stages,  an  initial  ZXM  stage  and  a  PWM 
stage.  The  optimal  DGZ  coordinates  from  the  ZXM  stage  become  the  initial 
D6Z  coordinates  for  the  PWM  stage. 

The  CEDF  maximization  algorithm  was  verified  and  validated  using 
two,  three,  and  five- installation  target  complexes.  The  CEDF  maximum 
value  and  optimal  DGZ  locations  for  two  example  problems  were  also  com¬ 
pared  with  results  from  NUCWAVE.  NUCWAVE  is  a  one-sided  nuclear  weapons 
allocation  war  gaming  model.  It  also  optimizes  the  damage  to  a  set  of 
installations  using  a  finite  number  of  weapons.  However,  NUCWAVE  deter¬ 
mines  a  sequential  optimal  solution;  it  optimizes  one  weapon  at  a  time 
until  no  increase  in  complex  damage  i3  possible.  The  comparisons  be¬ 
tween  the  results  from  the  CEDF  maximization  algorithm  and  from  NUCWAVE 
for  two  and  five- installation  target  complexes  indicated  that  the 
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algorithm  determines  the  same  local  maximum  as  NUCWAVE. 

A  symmetric  one  weapon-four  installation  complex  was  designed  to 
have  four  local  maximums.  CEDF  results  from  this  complex  were  analyzed 
and  pointed  out  two  ZXM  maximization  features,  an  operational  charac¬ 
teristic  and  a  limitation.  These  features  depend  on  the  two  types  of 
gradient  symmetry.  There  is  gradient  symmetry  from  either  geometrical 
symmetry  of  the  target  complex  or  collocation  of  two  or  more  similar 
weapon  types.  The  second  type  of  symmetry  is  a  limitation  and  prohibits 
ZXM  maximization  from  using  collocated  initial  DGZs. 

Further  analysis  of  three  and  four- installation  target  complexes 
indicated  the  presence  of  multiple  local  CEDF  maximums.  A  two  weapon- 
three  installation  target  complex  was  analyzed  using  CEDF  algorithm 
results.  There  were  two  distinct  local  CEDF  maximums  and  two  corres¬ 
ponding  pairs  of  optimal  DGZs.  Similarly,  a  two  weapon-four  installation 
complex  was  analyzed.  There  were  three  distinct  local  CEDF  maximums  and 
three  corresponding  pairs  of  optimal  DGZs.  The  CEDF  maximization  algo¬ 
rithm  located  these  local  maximums  using  different  initial  DGZ  conditions. 
For  both  complexes  analyzed,  one  local  maximum  was  definitely  the  highest 
valued  local  maximum  for  the  complex. 

Conclusions 


The  first  subobjective  was  to  determine  the  sensitivity  of  the  re¬ 
sults  of  the  CEDF  maximization  algorithm  to  the  mathematical  technique 
used  to  locate  the  optimal  DGZo.  Two  conclusions  of  the  study  empha¬ 
size  the  differences  in  CEDF  results  for  the  two  CEDF  models  and  their 
respective  optimization  techniques  —  ZXM  and  PWM. 

First,  the  algorithm  requires  some  indication  of  a  potential 
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increase  in  CEDF  value  in  order  to  move  a  DGZ.  That  is,  if  there  is  no 
indication  of  a  CEDF  increase  in  the  direction  of  a  valued  installation 
and  there  is  an  indication  of  a  CEDF  increase  in  the  direction  of  a 
lesser  valued  Installation,  then  the  algorithm  will  move  the  DGZ  towards 
the  lesser  valued  installation.  Eventually,  the  algorithm  may  converge 
to  a  less  valuable  CEDF  local  maximum. 

The  second  conclusion  is  that  a  weapon's  CEP  and  the  CEDF  model 
affect  the  CEDF  maximum  value  and  the  respective  optimal  DGZ  coordinates. 
All  three,  four,  five,  and  seven-installation  target  complexes  analyzed, 
that  used  weapons  with  CEP  >  0,  confirmed  this  conclusion.  ZXM  optimal 
DGZ  coordinates  were  not  optimal  for  PWM;  similarly,  PWM  optimal  DGZ 
coordinates  were  not  optimal  for  ZXM.  Results  indicated  that  each  CEDF 
model  located  a  unique  set  of  optimal  DGZs;  however,  the  distance  dif¬ 
ference  between  a  ZXM  and  a  PWM  optimal  DGZ  was  less  than  the  weapon's 
CEP.  This  difference  occurred  because  of  the  difference  in  Pd^  j  for 
the  two  models.  For  a  weapon  i-installation  j  interaction,  Pd^j  is 
larger  for  the  CEP-Excluded  model  than  it  is  for  the  CEP-Included  model. 
This  is  because  the  CEP-Included  Pd^  ^  is  reduced  by  a  factor  that 
depends  on  the  weapon's  CEP. 

The  second  subobjective  was  to  determine  the  sensitivity  of  the 
results  of  the  CEDF  maximization  algorithm  to  the  initial  DGZ  locations 
prior  to  optimization.  Four  initial  DGZ  conditions  were  evaluated  using 
three  and  four -installation  target  complexes.  The  four  initial  condi¬ 
tions  for  an  m -weapon  complex  were:  (1)  locating  the  weapons  at  the 
m  highest  valued  installations,  (2)  locating  the  weapons  at  the  m  hardest 
installations,  (3)  locating  the  weapons  at  the  complex's  centroid,  and 

(4)  locating  the  weapons  at  m  pseudo-random  points.  The  algorithm 
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using  these  initial  DGZ  conditions  located  more  than  one  local  CEDF 
maximum  for  three  and  four -installation  complexes.  Thus,  the  last  con¬ 
clusion  of  the  study  emphasizes  that  no  single  initial  DGZ  condition 
always  locates  the  most  valuable  local  CEDF  maximum.  Hence,  the  algo¬ 
rithm  does  not  include  a  decision  structure  to  determine  the  correct 
initial  DGZ  condition.  However,  this  conclusion  also  indicates  that  the 
most  likely  to  succeed  initial  DGZ  condition  is  to  use  the  coordinates 
of  the  m  highest  valued  installations  as  the  initial  DGZ  coordinates. 

This  investigation  characterized  three  factors  that  affect  the 
optimal  DGZ  locations  for  multiple  nuclear  weapons  in  a  target  complex, 
the  first  factor  was  gradient  synmetry;  this  symmetry  resulted  from  either 
a  geographically  symmetric  target  complex  or  collocated  weapons.  The 
second  factor  was  weapon  CEP.  Maximization  of  the  two  CEDF  models  pro¬ 
duced  slightly  different  optimal  DGZs;  this  difference  depended  on  a 
weapon's  CEP  and  the  CEDF  model,  the  third  factor  was  the  initial  DGZ 
location  prior  to  CEDF  maximization.  The  algorithm  located  different 
CEDF  local  maximums  depending  on  the  initial  DGZ  condition. 

Recommendations 

The  weapons  analyst  can  use  the  algorithm  to  solve  large  targeting 
problems  that  include  many  complexes  and  different  types  and  numbers  of 
weapons.  Ihe  algorithm  can  be  a  valuable  sensitivity  analysis  tool  to 
investigate  weapon  allocation  tradeoffs.  The  analyst  can  evaluate  the 
changes  in  total  complex  expected  target  value  damage  as  a  result  of  an 
increase  or  decrease  in  the  number  of  weapons  available  to  a  target 
complex.  Similarly,  the  analyst  can  estimate  the  effects  that  changes 
in  a  weapon's  yield,  CEP,  or  Pa  can  cause  to  the  optimal  DGZs. 

100 


*\j«  *w»twa/w*  -uk  wi  we  wuwt  jt  # 


w  .  w*  m' .  **v.  w  wv  m  ‘j  »■  ^  •  ■„  9  «,  »  -•  »  *  »>  *  »  •  *  *  e  *  <4  •  1 


The  CEDF  maximization  algorithm  does  have  some  limitations.  How¬ 
ever,  there  is  a  specific  improvement  or  recommendation  associated  with 
each  limitation.  The  following  reconmendations  would  provide  a  more 
capable  algorithm  for  strategic  targeting  studies: 

1.  Currently,  the  algorithm  accomplishes  only  two-dimensional 
location  of  optimal  DGZs;  the  user  provides  each  weapon's  height  of 
burst.  Optimization  of  each  weapon's  height  of  burst  could  be  added 
to  the  algorithm. 

2.  Currently,  the  algorithm  only  allows  military/industrial  tar¬ 
gets  that  are  modeled  as  point  targets.  The  algorithm  could  be  modified 
to  include  area  targets,  equivalent  area  targets,  and  target  avoidance 

areas. 

3.  In  an  analogous  manner,  the  algorithm  only  includes  blast  damage 
effects  for  these  point  targets.  Other  nuclear  weapon  damage  effects 
could  be  added  to  the  model. 

A.  Similarly,  other  optimization  techniques  could  be  used  to 
further  investigate  and  characterize  the  CEDF  local  maximumo  for  a 
target  complex. 

5.  User-specified  constraints  that  establish  a  minimum  acceptable 
Fd  for  some  or  all  installations  could  be  included.  This  addition  would 
provide  a  new  initial  DGZ  condition.  That  is,  locate  the  weapons  at 
the  installations  with  the  largest  minimum  Fd. 
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Appendix  A:  Determination  of  the  Distance 
Damage  Sigma  (qt)  and  the  Weapon  Radius  (WR) 

The  parameters  and  WR  are  necessary  to  calculate  i:he  probability 
of  achieving  a  specified  level  of  damage  to  installation  j  from  weapon  i. 

1.  Distance  damage  sigma*?,*  •  The  value  of  depends  on  the  T  factor 
of  an  installation's  VNIK  code.  Table  A-l  lists  the  T  factors  and 
their  associated  oj  values.  This  table  was  extracted  from  NUCWAVE 
Model  Methodology  Analysis  (Raf  25:3-4). 

2.  Weapon  Radius,  WR  ■  f (weapon  and  installation  parameters).  The  cal¬ 
culation  of  a  WR  dapends  on  the  concept  of  yield  scaling.  The  following 
information  on  yield  scaling  is  based  upon  Glass tone  and  Dolan's  presen¬ 
tation  in  The  Effects  of  Nuclear  Weapons  (Ref  11:100). 

"In  order  to  calculate  the  characteristic  properties  of  the  blast 
wave  from  an  axploslon  of  any  given  energy  if  those  of  another  energy 
are  known,  appropriate  scaling  laws  are  applied"  (Ref  11:100).  A 
1-kiloton  nuclear  explosion  is  the  reference  explosion  for  nuclear 
weapon  calculations.  Pressure  vs  range  data  have  been  tabulated  and 
graphed  for  the  1-kt  reference  oxplosion.  Also,  the  distance  scaling 
laws  use  the  cube  root  of  the  weapon's  yield  as  the  scaling  factor. 

That  is,  if  a  pressure  is  experienced  at  a  ground  distance  from  a 
1-kt  reference  explosion,  then  this  same  pressure  value  will  be  experi¬ 
enced  at  e  distance  d^  from  a  w-kiloton  explosion. 

"  ^(w)173 

The  needed  pressure  vs  range  (dw)  data  for  a  w-kiloton  explosion  can 
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TABLE  A-l 


Distance  Damage  Sigma  and  Target  Type 


T  Fad 

tor 

<* 

P  Targets 

Q  Targets 

L 

R 

0.1 

P 

S 

0.2 

M 

Q 

0.3 

N 

T 

0.4 

0 

U 

0.5 

be  determined  using  the  scaled  distance  d^.  Therefore,  to  determine 
either  the  overpressure  or  dynamic  pressure  from  an  explosion  of  w-kt, 
all  distances  need  to  be  transformed  to  the  1-kt  scaled  reference  frame. 
The  amount  of  pressure  an  installation  experiences  is  the  primary  deter¬ 
minant  of  the  installation's  probability  of  damage. 

Weapon  and  installation  parameters  are  needed  to  determine  the  WR. 
The  weapon  parameters  needed  are  the  yield  (Y)  and  the  height  of 
burst  (HOB).  The  HOB  and  the  subsequently  calculated  scaled  weapon 
radius  (SWR)  are  the  two  distances  that  need  to  be  yield  scaled.  HOB 
is  scaled  to  start  the  formulation;  after  the  SWR  is  calculated,  it  is 
inversely  scaled  to  specify  the  WR.  The  only  installation  parameter 
necessary  to  calculate  the  WR  is  the  VNTK  code.  The  following  presen¬ 
tation  is  based  upon  the  material  in  Mathematical  Background  and  Pro¬ 
gramming  Aids  for  the  Physical  Vulnerability  System  for  Nuclear  Weapons 
(Ref  6:57-61). 


WR  -  SWR 


1-tf 


(Y) 


1/3 


(A.l) 
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The  parameter  c  is  a  constant  that  equals  either  0,96  for  overpressure 
sensitive  installations  or  0.91  for  dynamic  pressure  sensitive  instal¬ 
lations.  The  SWR  is  calculated  from  the  scaled  height  of  burst  (SHOB) 
and  the  adjusted  VN  number  (VNfldj ) . 

SHOB  - 

(Y )1/3 


vNadj  -  VN  +  AVN 


An 


AVH-  5.485  ln(R)  for  overpressure  sensitive  installations 
-  2.742  ln(R)  for  dynamic  pressure  sensitive  installations 
iterative  procedure  is  used  to  calculate  the  VN  adjustment  factor  R 


R 


1 


K 

“Iff 


20 

Y 


1/3 


Re 


K  is  the  installation's  K  factor  and  the  exponent  e  equals  either  1/2  for 
overpressure  sensitive  Installations  or  1/3  for  dynamic  pressure  sensitive 
installations. 


SWR  -  expf^VNadj»SH0B^ 


The  function,  f(VNadj,SH0B),  is  a  polynomial  expression  whose  coefficients 
are  available  at  100-foot  increments  of  the  SHOB  between  0  and  900  feet. 
Hence,  the  SWR  is  derived  by  linearly  interpolating  between  a  low  SWR, 
that  is  calculated  from  a  low  SHOB,  and  a  high  SWR,  that  is  calculated 
from  a  high  SHOB.  For  example,  for  a  SHOB  of  632  feet,  a  SWR  low  is 
calculated  for  a  SHOB  of  600  feet  and  u  SWR  high  is  calculated  for  a  SHOB 
of  700  feet.  The  actual  SWR  is  a  linear  interpolation  of  the  high  and 
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low  SWRs.  Other  algorithms  use  a  table  look-up  with  parameters,  SHOB 
and  VN^j ,  to  specify  the  SWR  (Ref  19  and  25).  The  SWR  is  inversely 
yield  scaled  using  Eq  (A.l)  to  determine  the  WR. 
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Appendix  B:  Formulation  of  f(r)  and 
Calculation  of  the  Integration  Limits 


This  information  is  based  on  material  in  Mathematical  Background 
and  Programming  Aids  for  the  Physical  Vulnerability  System  for  Nuclear 
Weapons  (Ref  6:69-75). 

The  probability  of  achieving  a  specified  level  of  damage  to  instal¬ 
lation  j  from  weapon  i  depends  on  weapon  and  installation  parameters. 


(B.l) 


where  Pd(r)  -  distance  damage  function 
or  ■  CEP/1.1774 

y>(r,£)  -  (r2  +  s2  -  2rs  coa-©-)1//2 


A  closed  form  solution  to  Eq(B.l)  does  not  exist;  however,  an  analytical 


approximation  does. 


(B.2) 


-  (r2  +  s2  -  2rs  cosfr) 

«twc.  f(r)  -  P,(r)r  — («  ^  ”  L  do 

d  Lv, ac  j 

0 


(B.3) 


Eq(B,3)  can  be  rewritten  using  a  zeroth  order  modified  Bessel  Function, 


IQ(x)  (Ref  19:378). 


0 
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(3.4) 


£<r)  -  Pd(r>  ^  .  ** 


Wi,J  18 


calculated  using  normalized  distance  variables  and  Eq  (B.l). 


That  is,  r,  s,  and  WR  are  divided  by  t  ,  the  standard  deviation  of  the 
circular  normal  distribution. 

Let  r*  ■  ~  ■  the  normalized  distance  between  the  installation  and  the 

or 

impact  point 

dr'  ■  |dr  -  the  normalized  differential  element  of  r 

g 

a'  ■  =  *  the  normalized  distance  between  the  installation  and 

the  DGZ 

WR'  ■  5—  ■  the  normalized  weapon  radius 


Then  Eqs  (B.2)  and  (B.4)  become 


(f(r')dr' 


(B.5) 


-  (r')2  +  (s')2 


where  f(r')  -  Pd(r')r’  e 


Io(rV) 


(B.6) 


A  polynomial  approximation  of  the  zeroth  order  modified  Bessel  function , 
I0(rV),  specifies  a  value  of  f(r')  for  a  given  r’  (Ref  1:378). 

Eq  (B.5)  can  be  rewritten  as  a  definite  integral  with  limits  of  inte¬ 
gration,  a  and  b,  such  that  f(r')  -  0  for  r'  <  a  or  r'  >  b. 


Pdi>J  •  \f(r')dr' 


(B.7) 


Therefore,  it  is  necessary  to  determine  a,  the  smallest  possible  value 


Figure  B-l  Weapon-installation  geometry  to  determine  a  and  b 


of  r',  and  b,  the  largest  possible  value  of  r\  These  limits  depend  on 
weapon  accuracy  limitations  and  distance  damage  limitations.  The  weapon- 
installation  geometry  necessary  to  determine  the  limits  of  integration, 
a  and  b,  is  shown  in  Figure  B-l.  Ttoo  cases  are  examined. 

Case  I.  The  normalized  distance  between  the  installation  and  the 
DGZ  is  less  than  4  (0*s'*4).  The  distance  between  the  closest 
possible  impact  point  and  the  installation  is  r'  -  0.  The  dis¬ 
tance  between  the  farthest  possible  Impact  point  and  the  instal¬ 
lation  is  the  minimum  of  either  r'  ■s'  +  4  or  r'-R'- 
1.06*WR’exp(2.86* ).  The  point  at  r*  ■  s'  +  4  corresponds  to 
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the  maximum  distance  from  the  DGZ,  $>'  ,  that  a  weapon  could  be 
expected  to  impact.  The  probability  that  a  weapon  would  impact 
farther  than  ^*4 or  is  less  than  0.00005.  Similarly,  the  point 
at  r'  •  R*  corresponds  to  the  maximum  distance  from  the  instal¬ 
lation,  r',  that  the  weapon  could  detonate  and  expect  to  damage 
the  installation.  The  Pd(r)  for  r'  >  R'  is  less  than  0.0005. 

R'  may  be  either  greater  than,  equal,  or  less  than  s'  +  4. 

Case  II.  The  normalized  distance  between  the  installation  and  the 
DGZ  is  greater  than  4  (s'  >  4).  The  distance  between  the  closest 
possible  impact  point  and  the  installation  is  r'  *  s'  »  4.  The 
distance  between  the  farthest  possible  impact  point  and  the  instal¬ 
lation  is  again  the  minimum  of  either  r'  *  s'  +  4  or  r'-  R'. 

Therefore,  a  ■  max  (0,  s'  -  4)  and  b  -  min  (s'  +4,  R')  (Ref  6:73). 
Eq  (B.7)  can  be  evaluated  using  Gauss-Leg andra  Quadrature  between  the 
limits,  a  and  b, 


(B.8) 


where  r'. 

k 


(b  -  a)  (a  +  b) 
° u  zz,  +  : 

2  k  2 


and  f(t’,j5)  is  Eq  (B.6).  Gauss-Legandre,  the  quadrature  points,  zz^, 

end  the  coefficients,  ww^,  are  explained  in  Appendix  C. 

Eq  (B.8)  is  evaluated  to  determine  the  probability  of  achieving  a 
specified  level  of  damage  to  installation  j  from  weapon  i. 
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;ure  to  Integrate  f(r') 


Gauss -Legendre  quadrature  is  a  numerical  integration  technique  that 
approximates  a  definite  integral  as  a  finite  series  (Ref  27:125).  Each 
term  is  a  weighted  function  value. 

The  series  to  approximate  a  definite  integral  along  the  interval  [-l,l] 
by  Gauss -Legendre  is 

1  n 

\  f(y)dy  -  ^  w*k*f(yk)  +  Rn  (C.l) 

-1  k-1 

Where  wwk  ■  quadrature  coefficients 

yk  -  quadrature  base  points 

R  ■  remainder  (negligible) 
n 

n  *  number  of  quadrature  points 

Gauss-Legendre  integration  differs  from  symmetric,  trapezoidal  numerical 
integration.  In  Gauss-Legendre,  the  distances between  the  yk  base  points 
along  the  abscissa  are  not  equal.  Ihe  points  are  spaced  symmetrically, 
yet  unequally,  with  respect  to  the  midpoint  of  the  interval  [-l,l]  ,  the 
origin.  See  Figure  C-l.  This  method  is  more  efficient  than  equal 
spacing  trapezoidal  methods  because  it  requires  fewer  function  evaluations 
to  achieve  comparable  accuracy  (Ref  15:378).  Ihe  quadrature  coefficients 
wwk  are  positive  numbers  between  0  and  1;  they  are  weights  for  the 
f(yk)  values. 

Ihe  quadrature  coefficients  and  weights  are  calculated  from  the  nth 
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Legendre  polynomial,  Pn«  The  base  points,  yfe,  are  the  k  -  1,..,,10  zeros 
of  P  (y).  Similarly,  the  coefficients,  ww,  ,  are  calculated  from  P  (y) 

u  k  n 

and  the  y^  (Ref  1:888).  Table  C-l  lists  the  base  points  and  coefficients 
for  n  ■  10  (Ref  6:74  and  27:131). 

However,  to  calculate  Pd  the  interval  of  integration  is  not  j  -1,1*1 , 

i,j  J 

but  rather,  it  is  [a,b] .  For  this  interval  of  integration,  the  quadrature 
base  points,  coefficients,  and  limits  of  integration  specify  the  trans¬ 


formed  variables,  and  r^  (Ref  1:887).  Eq  (C.l)  becomes 


Vf(rk) 
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Quadrature  Base  Points  and  Coefficients 


n 

yk 

v,vk 

1,  10 

Q. 9739065285 

0.0666713433 

2,  9 

0.8650633667 

0.1494513442 

3,  8 

0.6794095683 

0.2190863625 

7 

0.4333953941 

0.2692667193 

5,  6 

0.148874339 

0.2955242247 

where  wk  - 


r  m  ^  *  a)..  +  (a  +  b) 


2  'k 


This  finite  series  is  used  to  calculate  Pd  .  For  each  k 

i»J 

an  and  associated  f(rk)  is  calculated. 


f(rfe)  «  f(r')  -  Pd(r')r"  exp 


<r')2  ♦  (s')2 

2 


I  (r\e ') 
o 


Then  f(rk)  is  lighted  by  w  and  sunned  to  form  Pd^  4. 


Appendix  D:  Computer  Code  of  the 
CEDF  Maximization  Algorithm 

This  appendix  presents  a  glossary  of  the  FORTRAN  variables  and  lists 
the  FORTRAN  V  source  code  of  the  CEDF  maximization  algorithm.  The  algo¬ 
rithm  has  the  capacity  to  handle  a  target  complex  with  up  to  ten  weapons 
and  tea  installations.  These  capacities  can  be  increased  by  changing 
the  array  dimension  variables,  MXM,  MXN,  MTM,  MIN,  MSQ,  and  M12.  Para¬ 
meter  statements  assign  values  to  these  variables;  the  glossary  describes 
these  variables. 

The  algorithm  has  several  options.  The  user  can  specify  one  of 
three  initial  DGZ  conditions  and  convergence  parameters  for  ZX06R. 
Chapters  IV  and  V  discuss  user  guidelines  in  selecting  a  particular 
option.  Also,  appendix  E  provides  instructions  in  how  to  create  an 
input  data  file. 

The  subroutine  INITLZ  assigns  values  to  the  two  PWMIN  maximization 
convergence  control  parameters,  E  and  ESCALE.  Chapter  V  discusses 
these  parameters.  Minor  code  changes  would  be  required  to  modify 
either  of  the  parameters. 

A  typical  two  weapon-four  installation  CEDF  maximization  problem 
requires  approximately  3.0  seconds  of  execution  time.  The  current  pro¬ 
gram  requires  50,000  words  of  core,  memory  on  a  Control  Data  Corporation 
(CDC)  6600  Cyber  computer. 

Glossary  of  Variables 

A[WRADS]  -  a  coefficient  that  is  used  to  adjust  the  VN  number. 

A[PDAM]  -  the  lower  limit  of  integration  for  the  Gauss-Legendre 

quadrature. 
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AA 

ACC 

ACC2 

B 

BB 

BESB 

and  BESC 

BEIA(j) 

CC(j) 

CEP(i) 

CEPA(i) 

CEPS(i) 

DFFRED 

DGLN(i,3) 

DGLT(i,3) 

DLNC(i) 

DLTC(i) 

DNCD(i) 

DTCD(i) 

E(i) 

EC 

EF 

ERF 


-  the  difference  between  weapon  i  and  installation  j 
x  coordinates  in  feet. 

-  the  ZXGGR  convergence  control  parameter. 

-  the  second  stage  ZXCGR  convergence  control  parameter. 

-  the  upper  limit  of  integration  for  the  Gauss -Legendre 
quadrature. 

•  the  difference  between  weapon  i  and  installation  j 
y  coordinates  in  feet. 

-  polynomial  approximations  of  the  modified  zeroth  order 
Bessel  function  that  are  used  to  determine  Pd^  ^ . 

-  the  beta  value  for  installation  j. 

-  an  intermediate  value  that  is  used  to  calculate  the 
gradient  of  the  CEDF. 

-  the  circular  error  probable  of  weapon  i  in  feet. 

-  the  adjusted  circular  error  probable  of  weapon  i  in  feet. 

-  intermediate  storage  of  the  circular  error  probable  of 
weapon  i  in  feet. 

-  an  estimate  of  tlie  expected  increase  in  the  CEDF. 

-  the  degree-minute-8econd  longitude  coordinate  for  weapon  i. 

-  the  degree-minute-second  latitude  coordinate  for  weapon  i. 

-  the  east/west  direction  from  the  prime  meridian  for 
weapon  i. 

-  the  north/south  direction  from  the  equator  for  weapon  1. 

-  the  longitude  coordinate  in  degrees  for  weapon  i. 

-  the  latitude  coordinate  in  degrees  for  weapon  i. 

-  the  PWMIN  convergence  control  parameter. 

-  the  value  of  CEDF. 

-  the  value  of  -CEDF  from  PWMIN. 

-  a  polynomial  approximation  of  the  error  function  that  is 
used  to  determine  a  distance  damage  function  value. 


ESCALE 

EV(j) 

EX 

F 

F(5,2) 

FACTOR(j) 

BN 

IV 

G(2) 

GR(i) 

GRAD(i) 

H 

HOB(i) 

I 

IER 

ILNC(j) 

ILTCCj) 

INASG(j) 

INC 

INCD(j) 

ISHOB(i,2) 

ISLN(j,3) 


the  maxinun  step  size  multiplier  for  a  single  step  of 
each  x. 

the  expected  target  value  damage  contribution  to  EC  from 
installation  j. 

an  exponent  that  is  used  to  calculate  the  VN  reduction 
factor. 

the  value  of  -CEDF  from  FUNCT  and  GFUNCT. 
the  ten  Gauss -Legendre  function  evaluations. 

PS(j)  *  value(j). 

the  sum  of  the  ten  Gauss-Leg endre  function  evaluations, 
the  value  of  -CEDF  from  ZXOGR. 

a  polynomial  expression  that  is  used  to  determine  an  SWR. 
the  2m  gradient  elements  of  the  CEP-Excluded  CEDF  model, 
same  as  GR(i). 

an  intermediate  value  that  is  used  to  determine  the 
appropriate  Bessel  Function  approximation,  BESB  or  BESC. 

the  height  of  burst  for  weapon  i  in  feet. 

generally,  the  subscript  of  a  weapon  array. 

an  error  code  from  ZXOGR. 

the  east/west  direction  from  the  prime  meridian  for 
installation  j. 

the  north/south  direction  from  the  equator  for  installation  j. 

a  logical  decision  variable  that  indicates  whether  an 
installation's  coordinates  have  been  assigned  to  a  DGZ. 

the  user-specified  indicator  variable  that  controls  the 
assignment  of  initial  DGZ  locations. 

the  longitude  coordinate  in  degrees  for  installation  j . 

the  low  and  high  scaled  heights  of  burst  for  weapon  i 
in  feet. 

the  degree-minute-second  longitude  coordinate  for 
installation  j. 
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ISLT(j,3)  -  the  degree-minute-second  latitude  coordinate  for 
installation  j. 

ITCD(J)  -  the  latitude  coordinate  in  degrees  for  installation  j. 

ITERS  -  the  number  of  ZXCGR  calls  of  function  GFUNCT. 

J  -  generally,  the  subscript  of  an  installation  array. 

JJ  -  the  argument  of  the  modified  zeroth  order  Bessel  functions, 

BESB  and  BESG. 

K(j)  -  the  K  factors  for  installation  j. 

KK  -  the  incremental  contribution  to  grad(i)  or  grad(i-Hm) 

from  installation  j. 

KO  thru  K7  -  the  eight  coefficients  of  the  polynomial  G(2). 

LNGMN  -  the  length  of  one  minute  of  longitude  in  feet. 

M  -  the  number  of  weapons  for  the  target  complex. 

MAXFN  -  the  maximum  number  of  function  evaluations  ZXCGR  is 
authorized. 

MAX1T  -  the  maximum  number  of  iterations  PWMIN  is  authorized. 

MSQ  -  the  dimension  variable  for  a  PWMIN  work  array  [-2m  *  (2m  +  3)]. 

MTM  -  the  dimension  variable  for  all  2m-element  weapon  arrays. 

MIN  -  the  dimension  variable  for  all  2n-element  installation 

arrays. 

MXM  -  the  dimension  variable  for  all  m-element  weapon  arrays . 

MXN  -  the  dimension  variable  for  all  n-element  installation  arrays. 

M12  -  the  dimension  variable  for  a  ZXCGR  work  array  [-  12  *  m]. 

N  -  the  number  of  installations  in  the  target  complex. 

NON  -  an  intermediate  value  in  calculating  the  gradient  of  the 

CEDF. 

NS  -  the  standard  deviation  scaled  distance  between  weapon  i 

and  installation  j. 

NW  -  the  dimension  variable  for  a  PWMIN  work  array. 

NWR  -  the  standard  deviation  scaled  weapon  radius. 
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l*la  \ 


I'V  v— » 


N2 

ONCE 

ORLN 

ORLT 

Mi) 

PDAM 

PDMG(i,j) 

PDR 

PP 

PS(j) 

R[PDAM] 

R[PDR] 

REM 

RR 

RI 

R2 

S 

SHOB(i) 

SIG 

SIGMA(j) 

SWR 

T(J) 


-  the  dimension  of  the  xy  weapon  coordinate  vector  x. 

-  a  logical  decision  variable  that  controls  the  algorithm 
so  that  ZXOGR  runs  a  second  time  with  a  smaller  ACC. 

-  the  longitude  coordinate  in  degrees  for  the  origin  of 
the  XY  coordinate  system. 

-  the  latitude  coordinate  in  degrees  for  the  origin  of 
the  XY  coordinate  system. 

-  the  probability  of  arrival  for  weapon  i. 

-  the  probability  of  achieving  a  specified  level  of  damage 
to  installation  j  from  weapon  i. 

-  Pd.  i  -  the  probability  of  achieving  a  specified  level  of 
damage  to  installation  j  from  weapon  i. 

-  the  distance  damage  function  value. 

-  the  distance  damage  function  value  used  to  calculate 
Pd^  j  for  the  CEP- Included  CEDF  model. 

-  the  probebility  of  not  achieving  a  specified  level  of 
damage  to  installation  j. 

-  the  standard  deviation  scaled  distance  between  the  impact 
point  and  the  installation. 

-  a  distance,  either  S  or  R,  from  the  subroutine  PDAM. 

-  a  remainder  that  is  used  to  translate  the  final  DGZ  coor¬ 
dinates  from  feet  into  degree-minute-second  coordinates. 

-  the  flat  earth  distance  between  weapon  1  and  installation  j. 

-  an  intermediate  value  that  is  used  to  calculate  the  VN 
reduction  factor. 

-  the  VN  reduction  factor. 

-  the  known  distance  between  weapon  i  and  installation  j. 

-  the  scaled  height  of  burst  for  weapon  i  in  feet. 

-  the  square  root  of  the  quantity,  (1  -ov^). 

-  the  distance  damage  sigma  for  installation  j. 

-  the  scaled  weapon  radius. 

-  the  T-factor  for  installation  j. 
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TLAT 

U 

UU 

V[PDAM] 

V[PDR] 

V[WRADS] 

VAL 

VALUE(j) 

VIND 

VN(j) 

VNA 

VNI 

W(5)[PDAM] 

W[WRADS] 

WR(l,j) 

WW(M12) 

WWW(MSQ) 

X(i) 

X(i  +  m) 
XX(j) 

XX(j  +  n) 

XXX 

X4 


the  sum  of  all  installation  latitude  coordinates  that  is 
used  to  determine  the  number  of  feet  per  minute  of  longitude. 

the  argument  of  the  error  function. 

an  intermediate  exponent  that  is  used  to  calculate  the 
gradient  of  the  CEDF. 

the  durniy  argument  of  BESB  and  BESC. 

the  duirny  argument  of  ERF. 

the  change  in  an  installation's  VN  number  when  it  is 
subjected  to  yields  other  than  20-kt. 

the  current  value  of  the  highest  valued  installation. 

the  value  of  installation  j. 

the  subscript  of  the  current  highest  value/hardest 
installation. 

the  integer  VN  number  for  installation  j. 
the  adjusted  VN  number. 

the  current  VN  number  of  the  hardest  installation. 

the  Gauss -Legendre  quadrature  coefficients. 

the  low  and  high  SWR  that  are  linearly  interpolated 
between  to  determine  the  actual  SWR. 

the  weapon  radius  for  the  weapon  i-installation  j 
interaction. 

a  ZXOGR  work  array. 

a  PWMIN  work  array. 

the  X  coordinate  of  weapon  i  in  feet. 

the  Y  coordinate  of  weapon  i  in  feet. 

the  X  coordinate  of  installation  j  in  feet. 

the  Y  coordinate  of  installation  j  in  feet. 

same  as  X(i)  and  X(i  +  m). 

the  (X^,Y^)  coordinates  of  the  m  weapons  after  ZXGGR  maxi¬ 
mization,  but  before  PWMIN  maximization.  Used  with  the 
mixed  CEDF  maximization  algorithm. 
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YIELD(i)  -  the  yield  of  weapon  i  in  kilotons. 

Z5  -the  Gauss -Legendre  quadrature  base  points. 

ZZ  -  a  standard  normal  random  variable. 

Source  Code 

The  next  27  pages  list  the  FORTRAN  V  code  of  the  CEDF  maximization 
algorithm. 


119 


m>v»  oonOM  o  n  o  n  ononn 


PROGRAM  OPTMZ 


C»**rt****vf  **«•>•.  '  ►  *  *  ft  ,#**•**♦**  4-  •*»*%•*••**»*****•«**#• 

OPT  MZ  IS  TH<„  DRIVER  MODULE  FOK  T'.iE  COMPLEX  EXPECTED 
OAMAGE  FUNCTION  <CEOF)  MAXIMIZATION  ALGORITHM,  THE 
ALGORITHM  DETERMINES  THE  OPTIMAL  rtfz  LOCATIONS  FO*  A 
FINITE  NUM8ER  OF  NUCLEAR  WEAPONS  AGAINST  INSTALLATIONS 
IN  A  TARGET  COMPLEX  BY  MAXIMIZING  THE  CEOF. 


PAR  AMEtER<MXM=J.O  •HXNslO  #MTM=20  «  MTNS20  » MSQaAZG  «  Ml  2  =  1 20  > 

EXT  IRINAL  GFUNCT 

INTEGER  M*N»VN<  MXN) (KCMXN)*UR CMXM*MXN ) *N2 *MAXFN • I ER *MAXlT»NU 
REAL  VALUE! MXN) (YIELOCMXM )*CEP( MXM) ,HOB< MXM) ,PA < MXM ) « SI 6MA< MXN ) 
REAL  X(MTM) #XX(MTN) *GR ( MTM) * FV* BETA ( MXN) » UW(M12 ) «  ACC,DFPREO 
REAL  ECMTM) «UUW( MSO )« XXXCMT M) t XM MT M) , ESCALEi EF » CEPS ( MXM ) ,ACC2 
LOGICAL  ONCE 
CHARACTER  T (MXN ) 

COMMON/INSTLN/  N ♦  VALUE  t VN *K  t  XX 
COMMON/PQIND/  T 

COMMON^UPNS /  M*YIELO»CEP*HOBfPA 
COMMON/PARA  MS /  UR  «S I GMA (BETA  *IT£RS 
COMMON/CNTRL/  N2  t AC C i ACC2 (DFPREQ ( ESC  ALE (E 

URI  TE ( 6  « 1 0  i  ) 

"0  F0RMATUXf3:  «•*•)  ) 

URITEC&tO  *  CEOF  MAXIMIZATION  ALGORITHM* 

URI  TEC61IOJ  ) 

ONCE*. TRUE. 

CALL  INITLZCX) 

CALL  URAOS 

00  10  C*1 «N2 
XXXCI  )»X(  I) 

S  CONTINUE 

CONJUGATE  GRAOIENT  OPTIMIZATION  OF  THE  CEP-EXCLUDEO 
CEOF  MODEL 

URI  TE  <  6  (1  33  ) 

WRI  TE<6(*  )  »  ZXCGR  MAXIMIZATION* 

ur:  tecs (i 03 ) 

ITERS*. 

MAX FN* l 00 
00  20  I al  »M 

CEPSCI  )  »CEP  Cl  ) 

C  EPC I ) *0 • • 

■J  CONTINUE 

3  CONTINUE 

WRITECSt*  )  »  ACC*  *  (ACC 

CALL  ZXCGRC GFUNCT (N2 (ACC fMAXFN(DFPREO*X(GR(FV(UU (IE R ) 

URI TC (6 f * )  •  IERa  *  d ER 

URI  TE<6(*  )  *  FUNCTICN=  •  *-F  V 

URI TECS  »* )  *  FUNCTION  EVALUATIONS:  *tITERS 
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1  ^ ; j f jr r«x v ,v ,-j ,cj  v.  *  j 


onoooo  u  *  oooo  u  ooon 


3'.' 


00  30  I  “l  »N2 

URITECfi.'  »  •  XC*,I,">=  *,XC  > 

WRITE<3**»  •  GRA0<V«1«*V=  *  t  GR  (I  ) 

CONTINUE 
CALL  OUTOGZCX) 

CONJUGATE  GPADIENT  OPTIMIZATION  OF  THE  CEP-EXCLUDED 
CEOF  M00EL  USING  A  REDUCED  CONVERGENCE  CRITERIA 

I FC  ONCE  )  THEN 
DO  33  I *1 «N2 
X4 ( I ) =X (I  ) 

CONTINUE 
WRI  T£  <6  *1  00  > 

WR1TE<6«»>  •  ZXCGR  WITH  ACC  REDUCED" 

WRITE 

ONCE*. FALSE. 

A  CC*ACC2 
I TERSsl 
GO  TO  23 
EMU  IF 

POWELL'S  CONJUGATE  DIRECTIONS  OPTIMIZATION  OF  THE 
CEP-INCLUDED  CEOF  MODEL 

URI  TECSfiar.  I 

WRITEfGfO  •  PWMIN  MAXIMIZATION" 

ur:  tecs  »io:  > 

MAX  IT*3 . 3 
DO  40  I»l«M 

CEP<I  »*CEPS<I  » 

CONTINUE 
MW*  N2HN2+S) 

CALL  PWMINC  XXX  »E»N2 iEF« ES CALEiMAXIT »UWU*NW> 

UR!TE<6*«>  •  FUNCTION*  "*-EF 
DO  30  t a  1  •  N 2 

WRITEC6  #* )•  XXX<"*I  t* )*  * tXXXCI » 

CONTINUE 

CALL  OUTOGZ(XXX) 

MIXED  OPTIMIZATION  OF  THE  CEDF  MODEL-  T HE  OOZ 
COORDINATES  FROM  THE  FUST  CONJUGATE  GRADIENT  OPTIMI¬ 
ZATION  BECOME  THE  INITIAL  OGZ  COORDINATES  FOR  P DWELL  "S 
CONJUGATE  DIRECTIONS  OPTIMIZATION. 

URI  TECfitlJC  ) 

URI TE »  MIXED  TECHNIQUE  MAXIMIZATION’ 

WRt  T£  (£  *1  0 :  > 

CALL  PUMINf  X4«E*N2»EFtE3CALE»MAXlT«UUUfNU) 

URITE<3**»  "  FUNCTION*  "»-£F 
DO  73  I  *1  #N2 

WRITER**!  *  X4C"fZi">s  »fX4<I) 

7  ‘  CONTINUE 

CALL  OU  TOG?  (X4 ) 

P3  CONTINUE 

END 
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SU8R0UTINE  GFUNCTCN2»X*F,GRAD> 

1#  CALCULATES  THE  CEOF(X).  THE  COMPLEX  EXPECTED 
DAMAGE  FUNCTION*  FOP  M  WEAPONS  AND  N  INSTALLATIONS 
USING  THE  CEP-EXCLUDED  CEDr  MODEL. 

2.  CALCULATES  the  2M  ELEMENTS  OF  THE  6RADIENT  OF 
THE  CEOF(X). 


*■**.  .*•*»***#•#*#»**,*  *».**■».  ★*#**« 1 


PARAMETER (MXMalC f  MX,N»iP *MTMs20  * MTNs^O } 

INTEGER  N.N*VMMXN)  »K(MXN>  «WR  (NX M .MXN  >  *N2 

REAL  VALUE! MXN)* YIELD (HXM).CEP( MX M) *HOB< MXM) *  PA C MXM) 

REAL  XCMTM)  »XX( MVN I  ,EC*E V (MXN > *P$ CMXN  > .PDAM. $  I Sfl A (MXN > 

REAL  POMG(NXM»MXN  I »  HE'TA  (MXN  I  *F 

REAL  FACTOR  (MXN  > »  CC  (MXN  >  »  OR  AO  (  Mr  M>  *  A  A  ,BH  «RR  ,UU  t  KK  *  MOM 
CHARACTER  T (MXN > 

CJHMON/INSTLN/  N. VALUE .VN.K* XX 
COMMON/POTNO/  T 

C OM MON/ W^NS /  M. ¥1  ELD* CEP. HO 8. PA 
COMMON/PARA  MS/  WR  .3I3MA.BETA # ITERS 

FUR  MAT(  *  WEAPON  XT  COORDINATES:  C  •  *F7. 3  »  •  ♦  •  *F7.  .  *•  »  *  > 

EC»  ? •? 

ITERS*ITERS*?. 

00  60  Lal.M 

WRITE ( 6.13  C  >  i.X<L»vX(L+M> 

M  CONTINUE 

CALCULATED  THE  CEOF(X). 

00  5-0  Jnl.N 
PS  ’  Jlal.O 

do  k  z-n.M 

POMT'f  tf  JlaPOAMCI  *  J»X> 

PSiJ?'PS(U»-li.  ‘-PACZMPONGCX  *J>) 

1?  contis m:m 

ev(J)*<  "P$<j)>*  value(j> 

URI  TE  !  6 .*  >  *  P3(**J»M  =  *.PS«U» 

E  C»EC*EV<  J) 

FACTOR^)**  VALUE!  Jl 

CC(  J!>*2.5f662R2'DETA(J> 

2  CONTINUE 

writer.*:  •  ec *  •»£;: 

Fa-EC 


J2.7 


L. 


u  u  u 


I 


CALCULATES  "HE  GRADIENT  OF  THE  CEDF(X). 

00  40  I *1 *M 
ORAOms:.-  .0 
GR  AD(  !♦  M)  =0  .0 
□  0  3'  J«1  *N 
AAaXXCUI-XC,  ) 

0B=XX CJ^NI-X(IoM) 

RRaSQRT<AA**2^BB*#2> 

IFiRR  .LT.C.COl  >  RR-J1C.J 

UU“ABS(  <<1.  /3ETA(JM*L0GCWR(I»J)/RR  l-BET  A  (  J )  )  /I  •  41  42 133  > 
NOMsFACTOR(J>*PA<I )#EXP(~UU**2> 

KK- NOM/ ( <1 •  -PACn*P0HGCX«d))*CCtd>*RK»»2) 

ORAL'CI )  *SPAOCI  >  ♦KK*  A  A 
GR  AOv  I*P>*GR  ADi:-*M)*KK*BB 
3*  CONTIN  IE 

4'.  CONTINUE 

DO  30  L»3lifJ2 

GRAD'D*  -GRAO<L» 

«L  CONTINUE 
END 
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£******«*  *♦**>**«  t**fr->**»-«**w  i  *.**.**•"*  *•*****.' »***<r*r**«rft»*k*«»'»V 


SU3R0UTINE  FUNC T  <  N2  »X  »F  I 

C  U  CALCULATES  THE  CEDF(X)  FOR  M  WEAPONS  ANO  N 

C  INSTALLATIONS  USING  THE  CEP-INCLUDEO  MODEL. 


P ARAH£TER<MXM=1C #MXN=lu  *MTM=2Q  *MTN=20  > 

INTEGER  M.N.VNCMXN*  ,K<MXN)»UR (MXM  »MXN  >.N2 
REAL  VAUIECMXN) » YIELOC MXM) .CEPCMX M> »HOB( MXM > • PA < MXM > 
REAL  XCMTM) *XX< MTN ) t EC . EV (MxN) . PS <MXN >  .POAM.F *BET AC MXN \ 
REAL  PDM3<MXM*MXN).SIGMACMXN> 

CHARACTER  TCMXN) 

C  OM  MON/I NSTLN/  N. VALUE. VN.K .XX 
C CM  MON/PQ IND/  T 

COMMON/WPNS /  M. YI EL  O.CEF .HOB  *PA 
COM  MON/PARA  MS/  WR . S I GMA .BETA . ITERS 
C 

EC=? .3 

DO  20  J=1*N 

?  sc  j  )=i  «c 
oo  u  i=i. m 

PDMGCI. JlsPOAMCI »J«X> 

PSCJ)=PS<J>M1.:  -PA(I)*PDMGCI  «J> ) 

V  CONTINUE 

EVC JlaCi.  -PSCdl  )  > VALUECJ1 
E  C=EC ♦EVC  J) 

2'!  CONTINUE 
C 

F*-£C 

RETURN 

EM0 


I 


I 
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Cw  *  *  *  *  *  *  ,  **»*..*•  (**■«**  *  *  * 


SUBROUTINE  INITLZCX) 

C  1.  P E AOS  USER-SPECIFIED  WEAPON  ANO  INSTALLATION 

C  PARAMETERS  FROM  THE  EXTERNAL  FILE.  INDATA. 

C  2.  ASSIGNS  INITIAL  OGZ  COORDINATES  ACCORDING  T 0 

C  THE  USER  OPTION  VARIABLE »  INC. 

C  3.  TRANSFORMS  ALL  WEAPON  AND  INSTALLATION  COORDINATES 

C  INTO  FEET  RELATIVE  TO  A  COMMON  ORIGIN  IN  A  XT 

C  COORDINATE  SYSTEM. 

C  INITIALIZES  ACCURACY  AND  CONVERGENCE  CRITERIA 

C  FOR  THE  OPTIMIZATION  SUBROUTINES,  ZXCGR  AND  PWMIN. 

(>**»  ***»***»**«,  .*».*»**^*»*— ******  **.****»*•**»,*  »**»»»»*«,, 

PAR  AME7ER  (MXMS10  .MXN-13  ,MTMs20  * MTNs20  ) 

INTEGER  0GLN*MXM,3>  ,OGLT(MXM,3)  , I SLN * MXN ,3 > , I SLT * MXN » 3 ) 

INTEGER  VN<  MXN)  ,K*MXN),M,N,INC»'V2»VIND,VNI 

REAL  DNCD<MXMI*DTCO(MXM>*INCDCMXN)»ITCOCMXM> ,ORLN,ORLT 

REAL  YIELD* MXM) ,CEP*MXM> ♦HOB(MXM) .PA C MXM> , VALUE <MxN » 

REA  L  X(MTM) , XX* MTN ) .ACC  ,ACC2 .DFPRED , ESCALE , E C MT M ) , VAL 
REAL  NOGLN*NOGLY ,NISLN,NISLT  »TLAT .LNGMN 
LOGICAL  INASG(MXN) 

CHARACTER  DLNC(MXM) .DLTCCMxM > ♦ I LNC < MXN I « ILTC < MXN ) »T (MXN ) 
CHARACTER  PH«SS 

COMMON/INSTLN/  N .VALUE * VN*K » X X 
COMMON/PQINO/  T 

CO;iMON/WPNS/  M. YIELD. CEP. HOB. PA 
COMMON/CNTRL/  N2 .ACC. ACC2 .OFPREO .ESCALE .E 
COMMON/ORIGIN/  CRLN . ORLT .LNGMN 
10.  0  F0RMAT(2X.I*.l2tl2,Al  .4X,  13 , 1 2 , 1 2  .A  1  ,4X, 

•  F%  0.2X.F5.  T.2X.F6.  ,  ,1X,F*.2  > 

10  3  FOR  MAT*  5X«I  2,5X  »F5*  ' »2X.F5.r  , 2X » F6. 0 . 1 X .Ft .2 > 

110  FOR  MAT*  2X.  I  A.  12.12,  A1  .4X.I3. 12.12  .A1.4X.I2.A1  ,11,2  X»F6.,'  ) 

120  FORMAT*/, MNITLZ.*,/) 

130  FORMAT**  THIS  PROBLEM  USES  *,12.*  WEAPONS*,/) 

1 A  0  FORMAT**  WEAPON  YIELD  CEP  HOB  PA*) 

13P  FORMAT*/,*  THIS  COMPLEX  CONTAINS  *,1 2»*  INSTALLATIONS*,/) 

160  FOR  MAT* •  LONGITUDE  LATITUDE  VNTK  VALUE*) 

170  FORMAT*/,*  THE  XY  COORDINATES  OF  THE  INSTALLATIONS  IN  FEET*,/) 

1?0  FORMAT*/,*  INITIAL  OGZ  LOCATIONS  ARE  *,A28»/) 

C 

W»I  TE*6  .12:  ) 

OPEN(15»F I LE- *  I NDAT  A  •  ) 

REWIND  15 
C 

CC  READ  THE  USER  INITIAL  DGZ  COORDINATE  OPTION  VARIABLE, I \C  • 

C  IF  INC  =  1,  THEN  USER-SPECIFIED  COORDINATES 

C  IF  INC  =  2,  THEN  HIGHEST  VALUE  INSTALLATION  COORDINATES 

C  IF  INC  =  3,  "HEN  HARDEST  INSTALLATION  COORDINATES 

L 

READU3,*)  INC 
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o  o  o  r»  n  fo  noon  ooo 


READ  USER-SPECIFIED  weapon  parameters* 

READC15**)  M 
WPT  TEC  6  tl  3C  )  M 
WRI  TEC  6  tl  40  I 
00  13  I=i#M 

READC15  *1C3 )  CDGLNCI*L>*L=l*3l,DLNCCX)tC0GLTCItL>*L=I  *3>t 
*0LTCCI>,YIELDCI  >*CEPfI >*HOB<I>tPACI  > 

WRITEC&«135  )  I*YIELO<I>tCEPCI)*HOBCI JtPACI) 

CONTINUE 

READ  USER-SPECIFIEO  INSTALLATION  PARAMETERS* 

R  EA  DC  15  »*  }  tl 
ur:  TEC 6*1 50  li 
MP.I  TEC 6  »l Sv  > 

DO  23  Jsl#N 

A  E AD  <15  #1 13  I  (ISLN(J,L)*L  =  1»3),ILNCCJ),<ISLTCJ,L>,L=1 *3) t 
*ILrC<J>  *VNC  J)  •MJ)«KC  JltVALUE(J) 

WRITE  <6  ,113  )CISLNCJ*L>*L=ltT)*ILNCC  J >♦ C ISLTC J,L > »L=1  »3  )• 
*i LTCCJI tVNC J> ,TCJ)tKC J)* VALUE CJ) 

CONTINUE 

$HNC0*iaa  •; 

SMTCOslSO  *o 
TLATsO.l 

TRANSLATES  INSTALLATION  4  DEGREE-MI NUTE-SECCNO 
COORDINATES  INTO  DEGREES* 

DO  30  J=1  #N 

I FCILNCC J)«EO**E*  >  THEN 

INCDC J)  =REALCISLNCJflH^CREAL<TSLNCj*2n/6?.5  >♦ 
*CPEAICISLNC 4*31 »/3S  0*? I 
ELSE 

NISLN=-REALCISLNC Jfl > I 

INCDC  J)  -NISLN-CREALCISLNC  J*2  »  )/6C.  ‘i  )- 
•CREALCISLNC  J»3)  )/36  0*3  > 

ENDIF 

I  FCINCDC4)*LT«SMNCD)  SMNCDsINCDC  J  ) 

IF<IL7CC4)*EQ*‘M  >  THEN 

ITC0C4>=REALCISLTC4*1  MMREALCI SLT i J ,2 > > /GO  •  3  >♦ 

♦  C  RE ALC 1 3L TC J*3 > )/33  : j*0  > 

ELSE 

NISLT  =- REAL  Cl SLTCJ«1>) 

ITCDC 4>=NISLT-CREALCISLTC4*2» »/60*0 )- 
* C RE ALC I SLTC J.3> >/36  lUO ) 

ENDIF 

I  FCI  TCDCJKLT.SMTCDI  3MTC0=I  TCD  C  J  > 

TLATsTLAWTCDCJ) 

3'.  CONTINUE 
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\  Jfc*  \  ■Si  1  N  I  <V  ;  ■>  *V  f*m  l“«  . 


•#  V.V  V.\iL»i.kUaWli  V.  » 


uuuu  #  u  <t  uuuu  uuuu 


LNGMN=COSO  .017453232*TLAT/NI*6'80.G 
ORLNsIMT(SMMCD> 

I  F<  ORLN.LT*:.'.  >  ORLN=XNT<SMNCD-l.  » 

ORLT=INTCSMTCO> 

IF<  ORLT  «LT*C*C)  ORLT-INTC  SMTCD-1 *0 > 

TRANSFORMS  INSTALLATION  J  DEGREE  COORDINATES  INTO 
FEET  RELATIVE  TO  A  COMMON  ORIGIN* 

WRI  TE  <  6  tl  7"i  J 

00  40  0=1 tN 

XK«J»aC INCDC JI-ORLNI *Sa.O -LNGMN 
XXC  J«N)a<ITCD<JI-ORLT>*364'»0  3  *G 
WRITER*)  •  XX<*«d*M  =  •»XX(J>»*  XX< • * J+N t • » =  »*XX<J*N> 
CONTINUE 

IF( CINC.EQ.SI.OP. <INC.E0.3> >  THEN 
00  45  Jal 'N 

INASGIJJa. FALSE. 

CONTINUE 
END  IF 

ASSIGNS  WEAPON  I  INITIAL  OGZ  COORDINATES  ACCORDING 
TO  THE  USER  OPTION  VARIABLE  t  INC. 

I  FI (INC.NE.2) .AND. (INC.NE.3 >  >  THEN 

IF  INC  s  li>  THEM  TRANSLATE  USER-SPECIFIED  WEAPON  I  DGZ 
OEG REE- MI  NUT E-SECOND  COORDINATES  INTO  DEGREES. 

PH=«USER  SPECIFIED  • 

DO  5’  I*: «M 

IFIDLNC (I>.EQ.*E*>  THEN 

ONCO<I)aREAL<OGLNCIfl) > ♦ C REAL < DGL* < I  * 2 » > /G 3 • ? )♦ 

*  <REAL<OGLN< 1*31 >/3£  C *0 > 

ELSE 

NDGLNs-REAL<OGLN<T»l> ) 

ONCDf I >=NDGLN~CREAL<DGLN  <  I  »2 ) )/6Q.0>- 
*<REALC0GLN<I*3) >/36 3.0  ) 

ENOIF 

I F  (  □  LTC  <I).£Q.*N*J  THEN 

OTCD<I)=REALCDGLTCI.I  I  >♦  (  f.  EAL  (  DGLT 1 1  «  2 )  ) /SO  )♦ 

-  <REAL<OGLT<  1 13) 1/36. 0 .3  > 

ELSE 

NOGLT  =-REAL<DG|.T<I*l>  > 

DTCU<t>=NOGLT- <REAL<DGLT<I  ,£>>/ 61.0  1- 
«-  <REAL<DGLT<  1.3)  1/367  •  J  ) 

ENOIF 

CONTINUE 
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c 

C  TRANSFORM  WEAPON  I  DEGREE  COORDINATES  INTO  FEET 

C  RELATIVE  TO  A  COMMON  ORIGIN, 

C 

0  0  60  IalfM 

X<n  =  CONCOCI>-0»l.NJ*6-3,  , *LNGMN 
X<I*M>=<DTCO(I)-ORLT>* 36460  ,1 
6'.  CONTINUE 

C 

ELS £t F< I NC. EQ ,2  )  THEN 
C 

C  IF  INC  =  Zf  THEM  ASSIGN  THE  COORDINATES  OF  ’’HE  M 
C  HIGHEST  VALUED  INSTALLATIONS  AS  THE  INITIAL  DGZ 
C  COORDINATES  OF  THE  M  WEAPONS. 

C 

PHs'HIGHEST  VALUED  INSTALLATIONS* 

0  0  7.1  1*1  vM 
VALaQ.r 
VIND=0 
DC  75  J=1*N 

IF(IHASG(J)>  GO  To  75 
IFCVALUE(J>.LT.VAL>  GO  TO  75 
VAL*  VA  LUE(  J> 

VIND=J 

75  CONTINUE 

XCI >*XX(VIND> 

X(I^M)*XX<VINO+N> 

INASG<VINDI=.TRUE. 

?'  CONTINUE 

C 

ELSE 

C 

C  IF  INC  s  3«  THEN  ASSIGN  THE  COORDINATES  OF  THE  M 
C  HARDEST  INSTALLATIONS  AS  THE  INITIAL  DGZ  COORDINATES 

C  OF  THE  M  WEAPONS. 

C 

PHa*HAROEST  INSTALLATIONS  * 

00  &n  I*l*M 
VNI=0 
VINOsO 
00  85  JsltN 

IF< INAS  G<  J) >  GO  TO  85 
IF(VNfJ).LT.VNI)  GO  TO  85 
VNI *VN<  J) 

VINO* J 

85  CONTINUE 

X«n*XX(VIND> 

XCI«-N>  =  XXCVINO*N» 

INASG<VIN0)=.TRUE. 

8'  CONTINUE 

END  IF 

WPI  TE  <  6  *1 8C  >  PH 

DO  95  1*1** 

WRITE  (6t*  )  *  X<*»I**)s  **X<I>**  X<**I*M*#>=  **X«I» 

95  CONTINUE 
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INITIALIZE  THE  CONVERGENCE  PARAMETERS 
OF  THE  SUBROUTINES*  ZXCGR  AND  PUMIN. 

ACCsfUH 
REA0(15f»>  ACC2 
REA  DC  15  *•  )  DFPRED 
N2=2*M 
escale=5ooc 
DO  90  I=liN2 
Ed)=:.i 
CON  TINUE 
CL0SEC13) 

END 
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C 


SUBROUTINE  OUTOGZ(X) 

1*  TRANSLATES  THE  FINAL  06Z  COORDINATES  FROM  FEET 
INTO  0E6REE»MI MUTE -SECOND  LONGITUDE  AND  LATITUDE 
COORDINATES. 


PARAMETER  (MXNslC  ,MTM*20  ) 

REAL  X(MTM) ,ORLN* ORLT, REM. LNGMN 

REAL  YXELD(MXM) *CEP(MXM)*HOB(MXM),PACMXM) 

INTEGER  DGLN(MXM,3> *DGLT(MXM*3) 

CHARACTER  DLNC( MXM ) »DLT C (MXM  > 

COMMON/UPNS/  M» YIELD* CEP. HOB *PA 
COM  MON/ ORIGINS  ORLNiORLT  *LNGMN 
1'.  "  F0RMAT<5X»I2*5X*IR»I2fI£»Al»AX*I3*I2*I2»Al  ) 

110  FOR*ATC/o«  UEA»ON  LONGITUDE  LATITUDE*) 

120  FOR  MAT <  /) 

C 

WRI  TE(S*110  > 

DO  10  Xal*M 

I FCORLN.GE.O • T )  THEN 
DGLNC III) a  I NT  COR LN) 

DLNC< I) *■ E* 

DGLNCI*2)aINT(X(I  ) /LNGMN  > 

REMsXCI)-LNGMN*OGLN(I ,2> 

ELSE 

OGLNCIil )sABS(INT<0PLN*1.0)  ) 

OLMC(I)s*M* 

DGLNC I*2)*INT(6'.0-X<I> /LNGMN  > 

REM*60 .3 • LNGMN-X(  I l-LNGMN*06LN( I *2  > 

ENOIF 

DGLNC 1*3) =1  NT (REM* GO. ./LNGMN) 
rFCORLT.GE.O. " )  THEN 
DGLTCIil IsINTCORLT) 

DITCC  I ) s * N • 

OGLTC  I*2)=I  NT  (X  ( I +M ) /&>)  80  ) 

REMaX  (I+M)-6DPi)uu*DGLTCI  *2) 

ELSE 

OGLTCI*l)sABSCINTCOPLTTl.r )) 

OLTCC  I)a«S* 

DGLTCI*2)3IMTC6  .0  •*(!♦«  >/6  !  BJ.C-  ) 
REM»364aOO.o-X<:*-M)-6O03  .0  * DGLT C I *2  ) 

END  IF 

D GLT (1*3) =1  NT CREM/1 01.333 333) 

WRI TE(6*10:  )  I,(OGLNCI.L),L  =  l*3)fOLNC(I)*(OGLT(I*L) *Ls 
1’  CONTINUE 

up: tecs, 12: ) 

RETURN 

END 


1 *3) ,OLTCCI) 
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SUBROUTINE  WRADS 


C 

C 

c 


:•  calculates  the  weapon 

WEAPON  I  -  INSTALLATION 
2*  CALCULATES  THE  B ETA ( J 


R  A  0 1  i  *  UR< ! #  J) #  FC  3  EACH 
0  INTERACTION. 

>  FOR  EACH  INSTALLATION. 


C**«*  •»***••«•****••*****••»•****..*•*«*••*•***«**»  *******  ****** 

PARAMETER CHXM*1 C * MXNsld ,MTMs2Q » MTN=2C > 

REAL  KK»YIEL0(MXM)*H0BCMXM)#E*EX#A#3IGMA(MXN>*SIG#R1#R2 
REAL  V#VNA* SHOBCMXM ) «6(2)*U(2 >#SUR#K0  »K1 *K2 *K3 « K4 »K5# K6 « K7 
REAL  VALUE(MXN>  *CEP CHXM > ,PA CMXM > #BETA < MXN) »XX< MTN > 

INTEGER  VN(MXN) #KCMXN)#ISH0BCMXM*2> «UR(MXM*MXN>*I «J«L«M#N 
CHARACTER  TCMXN)#TT*5 
COMMON/ZNSTLN/  N# VALUE* VN*K , XX 
COMMON/PQINO/  T 

COM MON/UPNS /  M#YIELD#CEP#HOB»PA 
C OM MON/ PA RAMS /  WR #SIGMA#BETA»ITERS 
10  FORMAT!/# •WRADS*#/) 

C 

WRI TEC4  #1 J) 

C 

C  SPECIFIES  THE  DISTANCE  DAMAGE  SIGMA  ANO  THE  BtTA  ffOR 

C  EACH  INSTALLATION  J* 

C 

J*? 

20  Jsj  *1 

I  FC  CTC  J).  EG.  «L«>.OR.CT<.J).EQ*  •)»•>>  THEN 
S I QMA ( J  IsC  *1 

EL3EIF(  <T  (  J  )  •  EQ  •  *  P  *  >. OR  * CTC J >.EQ. *S • t  >  THEN 
SIGMAC Jla;*2 

EL5EIFC  CTCUI*EQ.*M»  >•  OR •  C  TC  J ).E Q.  *Q*  )  I  THEN 
S  ZGMACJ>sC*3 

ELSEIFC  CTCU).£Q.«N«  > • OR • CTC J>* EQ.  •TM)  THEN 
3  I QMACU  Ja ' • A 
ELSE 

3 IGMACU>*0.5 
ENDIF 

I  FC CTC J>*EQ.*L*  >*0R .CTC J>*EQ.*P • >.0R*  <T  I. OR. 

*  <  TC J).EQ. *N»).0R.<T(J).FQ.«0,>>  THEN 
T  T»f PTYPE  * 

ELSE 

T  T=*QTYPE • 

ENDIF 

I  FC  TT,EG.  * PTYPE  *  )  THEN 
EXsO.5 
3  I G®'  *96 
ftsS*485 
R 1=2.3 
ELSE 

EX=Cl.fl/3*) ) 
srGs:.9i 
AS2.742 
Rl»3.  .1 
ENDIF 

BET  A  C  J>  sS  QRT  C-LCG  C 1 • " -SI GMA 
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n  o  o  o  o  *-*  *-*  o  o  o  h  h  nonoH  o  ooo  m  noon 


I  =0 

3r  i  *m 

W  Cl  )-Z .  3 
U( 2  >~J • 0 
S  Cl  ISO* 0 
G  <  2  >  s  0  •  0 

CALCULATES  THE  VN  REDUCTION  FACTOR*  R2*  AND  THE 
ADJUSTED  VN  NUMBER*  VNA. 

KKsREAL (K  (  J  >  > 

0  R2=1«Q-CKK/1Q.3  H-CKK/1Q  .0  >*  (20.C  /YI  ELD  (I  >  >  *  *  C 1  •  0 /3  *  3  )*R1**EX 

IF! ABS(R2-Pi>.GE. 0.0011  THEN 
R1  =  R2 

so  to  io: 

END  IF 

VsA*L0G<R2) 

VNA  =REAL(  VN  (  J )  >  ■*■  V 

CALCULATES  THE  SCALED  HEIGHT  OF  BURST  CSHOB). 

S  HQ  B(  I  )  =  HOB  <1  J/TIELOCI  >**  Cl.  0/3*0  > 

I  F(  SHOB  (Z  )•  GT  «9P  D  «*i)  >  THEN 
WRITE <6**  )  *H0B  T CO  BIG • 

SHOBCIlsSPC  .w 
ENDIF 

ISHOBCI  *1 )  =  INT CSHOB (I  1/110 .Q  l*109 
ISHOBCI  *2  )- ISHOB  (  I  *  1  >  ♦  1 1 , 0 
Es(  S  HOB  Cl  >- REAL  (ISHOBCI  *11)  )/lrj  '•  .0 

L*3 

.0  L  =L  n 

DETERMINES  THE  POLYNOMIAL  COEFFICIENTS  TO  CALCULATE 
G (V NA « ISHOB  ) • 

IFC TT.EQ.*PTVPE*>  GO  TO  2C0 
GO  TO  210 

20  GCL  >sK3  «-Kl*  VNA*K2*VNA**2*K3*VNA**3«K4*VNA**4 
GO  TO  140 

30  GCL>sK':.’*Kl*  VNA»K2*VNA**2*K3*  VNA*  *  3»  K4*  VNA*  *  4*. 

•  K5*  VNA**5*KS*VNA«*6*K7»VNA*  *  Y 

CALCULATES  THE  SCALED  WEAPON  RADIUS  CSUR). 

40  WCL  >*EXPCGCL>> 

IF< E.LT. 0.00011  GO  TO  150 
IFCL.LT.2)  GO  113 
52  SUR=U(1)+E*(U(2)~W(1 )  > 

INVERSE  YIELD  SCALES  THE  SUR  TO  DETERMINE  THE  WEAPON 
RADIUS*  WRC I  * J) *  FOR  THE  WEAPON  I  -INSTALLATION  J 
INTERACTION. 

WRC  I*  J>sINT<CSWR.*YIELDCI  >*<«  <  1  .0  /3 .0  )  *SI  G/  C 1  •  '  -  SI  G«A  C  J  )*  •  2  )  >  +  0 . 5  ) 
URITEC6,*)  •  UM*»  :*«*•*  J*«  )=  »»WRCI*J> 

SO  TO  170 
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163  URITEC6,*)  *  VN  TOO  LARGE  FOR  HOO» 

UPC  I* J)=0 

170  IFCI.L7.M)  GO  TO  33 
I FC  J.LT .N  >  GO  TO  20 
GO  TO  220 

'.'0  IFC  tSHOB«I*L>.LT.10  >  7 HEN 

COEFFICIENTS  OF  GCVNA.ISHOB)  FOR  OVERPRESSURE  IP-TYPE) 
TARGETS  FOR  ISHOB  FROM  FEET  TO  903  FEET, 

I F( VNA . LE. 7 ,5 )  THEN 
K0=8.206936 
Kls-9 .8  6 622 2 2 E-  2 
K2a-4„27D5319E-03 
K3  =  44. 67361’*  OGE-O  3 
K4=n.O 
ELSE 

K0=e.263243 
Kls-l  ,2109524E-n 
K2=12 • 74266E-04 
K3*-9.2il  63496E-'/  6 
K4*1.C 
ENOIF 
GO  TO  120 

ELSEIFC  ISH0BCI*L),LT,2D'.  I  THEN 
rFCVNA.GT.31,.'  >  GO  TO  160 
I  FCVNA.  LE.7.S  I  THEN 
KO 38.29123 
K  la-1  .132939E-U 
K2a31.19908E-03 
K3s0.0 
K4a0.0 
ELSE 

KO  , 29  ?59 
Kl=-i  .1  343338E-C1 
K2S-4.6494C85E- . 4 
K3=63.63Q1E-G6 
K43-9  .1 680  378E-'1 7 
ENOIF 
GO  TO  12' 

ELSEIFCISHOBCItl.  >,LT.3C‘:  )  THEN 
I FCVNA. GT. 41. C l  GO  TO  160 
KOse.395223 
Kls-l  .4717856E-  -1 
K2S12.74489E-03 
K 3 3-2.0  632771 E-  '.3 
K4sl6.67391E-r.5 
K5=-6*89342E-06 
K6sl4.2  3714f.-'.' - 
K7  =  -l  .1  675C  1  St! - '  9 
GC  TO  13;-. 
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ELSEIF(ISHGB<IfL>.LT.v:-‘-.  >  THEN 
I FCVNA.GT.34.C )  GO  TO  160 
KG =6. 41 9304 
Kl=-9  .9027016F-i  2 
K2a-A.18727«7E-: 3 
K3  =  S4.49C'84E-05 
K4a-3.758352E-15 
K5sl4.1C969E-C7 
K6*-2.;:i7Q989E-CB 
K  7s  0 • J 
GO  TO  13° 

EL3  EIF(  I  5H0B(I  * L )  •  L T.  30 >  THEN 
IFCVNA.GT.30.0  I  80  TO  16u 
*5  =  8.499469 
Kls-1 •0963211 E-* 1 
K2=-3.4443747E-;3 
K3=72.617C6E-05 
K4  =  -7 .1 C  9C3E-C3 
K3  =  33 .1 9013E-1  7 
K6S-5.6635  j37E-"'3 
K7=1.C 
GO  TO  13.) 

ECSEIF(ISH0B<l*L>*LT.6i  >  THEN 
I  FCVNA.3T.27.  )  >  OC  TO  160 
K0  =  8«  323983 
Kl*-6.312G5S2E-"2 
K2a-2.5622191E-,'2 
K3*34.26447E-Q4 
K4S-3.926339E-14 
K3  =  34  .8  330  4E-0  S 
K&S-1.C228646E-  6 
K7sll.4432C-09 
GO  TO  130 

EL3EXFCXSH0B(ZtU«LT*70(->  THEN 
I F( VNA.  GT  .23. C )  GO  TO  160 
Kuxd.SB  6222 
Kls-1  •I?  0  2  711E-.il 
K2--9.9171739E-'.  3 
K3s26.P  232E-04 
.  K4=-3.6(,28224E-,:,4 
K3»20.C2315E-C6 
K6s-l«t 826364E-"6 
K7sl3.41337E-C-3 
GO  TO  13: 

CLSSIFCl3HOBUfU.LT.8G.:  I  THEN 
1FCWNA.GT.22.:  )  GO  TO  160 
KJa  3. 653962 
Kls-1  .367988  6E-1 1 
K2sl 4 .2  6201 E-G  3 
K3"*-4 929993E-0J 
K4a30»2H23E-C3 
K3S-2.5712239E-05 
K6S43.79003E-C8 
K7*0  # J 
5  0  TO  13.: 
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ELSEIFC  ISH08<I*D.LT.9C  )  THEN 
r  FIVNA*  0T.21*;:  )  GO  TO  160 
I F(VNA*LE.7*5 )  THEN 
K0 =3.681285 
K1*-1*1432A64E-U 
K2*“l  *78866-  6E-  '<3 
K3=13*9590?£«CS 
K4*0 • 0 
ELSE 

KO  *12  <*51342 
Kl*-1 *516344 
K2*17.69944E-02 
K3»-8  *9U0633E-03 
K4»14*C  0736E-G3 
ENOIF 
GO  TO  12f 
ELSE 

IFIVNA. QT*2C.  )  GO  TO  16C 
IF<VNA*LE.7*3>  THEN 
KO  >8*  71 9634 
Kl»-1  •  213ft326E-'ll 
K2«12.G  36C4E-04 
K3»-1*3«632S1E-14 
K4»0*  0 
ELSE 

KQ«13  *4  7289 
Kla-1  *971983 
K2*23*47267E-Q2 
K3»-l*4323113E-:2 
K4s26*4C371E~C5 
ENOIF 
00  TO  120 
ENOIF 

COEFFICIENTS  OF  G<VNA»ISMOB>  FOR  DYNAMIC  PRESSU»E 
<  Q«  TYPE »  TARGETS  FOR  ISHOB  FROM  !.  FEET  TO  95  C  FEE 

10  IF< ISHOB<I*L>*LT*l).  )  THEN 
I  F(WNA*  GT.33.r  )  SO  TO  16" 

KOaa.313139 
Kl=-0  *1  C 60 A 68 
K2»G.OO C 3224 
K3»-I3  *f,  ID  31  J 
K4«3.2264fE-05 
KV=-l*23227E>*C6 
K6»l  •  96  7 7E*  J  8 
K7a-t ;BeoE-n 
GO  TO  130 


EL3EIF<ISH0B(I ♦L>.LT.2-  1  THEN 
I  F<VN4.GT  .35. C  >  GO  TO  ISO 
K3*3.3763B2 
Kla-Q  ,1042945 
K2  =  -0. 0012014 
K'35-3.U136E-0  5 
K4=i.2&757£-05 
K3*-4.97579E-07 
K6*5. 77257E-3 9 
K7a0.0 
GO  TO  13C 

ELSElF<ISHC8<I»L>.LT,3i:  '  >  THEN 
I  F<VNA,3T,35.:;  >  GO  TO  160 
K0»8,42t'24 
Kls-1  ,P  °473E"C  1 
K2»14.6226flE-04 
K3S-3.969792E-14 
K4366,97002i-aS 
K3»-3.: 149459E-  )6 
K6=61.8622flt>C9 
K7s-4,K66633E-1 
SO  TO  13C. 

EL3EXFCl«JHOBCl»L).Lr,4r’  ‘  )  THEN 
l  P<  VNA,  ST  ,3S,  >  50  TO  160 
KQ  *8,  49  53 15 
Kla-». 1,031393 
K2»-0  ,D  ~  34114 
K3«0.0:  0  3087 
K4*-t.07267E-05 
K3*3.13662E-Q7 
K6*-3.3*646E-03 
K7  «■:).(} 

0  3  TO  ir 

CI,SEIP«:SH0B«t#L>*LT,9C  .• »  THEN 
I  F<VNA,GT,3U.  >  GO  TO  16’) 

K3 #3. 37 SO  03 
Kla-0  ,1C3«PP3 
K»«-0  .“''637M9 
K 3a 0,00 12362 
K4s-3  #0  “ C 1333 
K3»a*Ol397E-06 
K6*-2  .34684F-0  7 
«7«2.51293E-3V 
5  0  TO  1 3>: 

ELS  EXP<  ISHGS<  I  ft  l«tLT*60  •  >  THEN 
I FC VNA, GT ,2P,  '  )  GO  TO  160 
Kfl *4,643504 
Kl«-9  ,1110364 
K2*-0  #Sfl4».9t'4 
K3»3  •  OC  4  6644 
K4s**7  ,76848E«f(5 
K9«3*9e69ffE-36 
K6»«2,27C79i>f'T 
KT=3,:i0626F-09 
GO  TO  13" 
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ELSE  IF<  ISHCBf!  »L>.Lr.  7v  )  THEM 
I  FCVNA.GT.26...  )  GO  TO  160 

«:=£• 636697 
Kis-0  .1164822 
K2*3*  3C  ‘3634 
K3=-0  .0  0061  69 
K4*?l.57541E-S3 
K3S-4.L  7263E-06 
K6=5. 66  4C2E-1 8 
K7srt.O 
GO  TO  130 

ELS EXF(ZSHOB(Z«L).LT.8Qj>  THEN 
I  FCVNA.GT.23.'  J  GO  TO  160 
KO  *8*  70  7449 
Kls-0  .1 1753C  2 
K2*iU 00  23  48  3 
K3»-0  •<fl3034 
K4a0.00rl9C« 
K3*-i«152C0E-03 
KC  *2* 330  79E-0  7 
K7b-2.447Ci4£-09 
30  TO  133 

ELS  E£F( I S HCB(  I  *1  > .LT,  ?(’ .  »  THEN 
I  F<VNA.  GT.23.ij  )  GC  TO  160 
K038. 736328 
Kl*-0.U316j3 
K2«3. 2021175 
K33-3. 0013218 
K4*  0 ■ QC  0  2634 
K9*-l«967B9E-w5 
K6aft.18P13F.-07 
K7s-7  *20  362E-C  9 
00  TO  1  3: 

ELS  E 

I  F< VNA.  GT  *22.i  )  G  3  TO  163 
K0a8. 7930  42 
Kla-3  .1134*3!’ 5 
K2*1) .02  01871 
K3»-3«0 C  1105 8 
K4=0.3C 02337 
K5S"2  •■’1 1 362E'-,t  5 
K6=6. 97  320  E-C 7 
K7S-8.74866E-C9 
GO  TO  13) 

F.NOUF 

220  CONTINUE 
250  FOft  MAT < / /  1 

Mfit  TE<6*2  JO  ) 

END 
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C5***ri»***  f*v**«**«***-**'ft*-t'*»»  ••****#*•*»**  ***»4n-*»ar****t»**t*»* 


REAL  FUNCTION  PDAM<I*dtX> 

U  CALCULATES  THE  PROBABILITY  OF  ACHIEVING  A 

SPECIFIED  LEVEL  OF  DAMAGE  TO  INSTALLATION  d  FROM 
WEAPON  I* 


PAR AMETER <MXMsi:,MXNsi: »MTM=20 * MTN=20 ) 

PARAMETER<R95sO.O  > 

COEFFICIENTS  OF  THE  POLYNOMIAL  APPROX  IMATI IN  OF  THE 
MODIFIED  ZEROTH  OPOER  BESSEL  FUNCTION, 

PAR  AMETER<Hl*3.*I56  223tB2s3.!)a9«424,B35l.2  )6T4?2t 
♦B  Mi.  ,2639  73  2*05  =  0,0  36C76S  *B6=0  J  4561  3,  CO  ='.  .39894226  « 

•  C  la  r,,G  13285  *2  ,C2=G  •  *•  J 22 5319  ,C3=  ,CC  1  57563  *04  =  3  •  3  3  916281  * 

*C  5=0  ,0  20  577G6  ,C6*&  •  •.2*353  37  *C7=:  ,  v  l  647633  »C3  =0  .  0  3  392377  > 
INTEGER  U«<MXM,MXN>  *2  t J *K2*L2  »M »N  *VN  CMXN)  » K <MXN  > 

REAL  BESBtBESCt V*Z<5 ) * U <5 ) ♦ S ♦ X < MTM) , XX C MTN ) * R ,BET A C MXN ) , R95 
REAL  S I OMA< MXN ) ,CEP <MXM I iCEPA tPDM *PP ,PDR ,NUR *NS , VALUE C MXN > 
REAL  A,B,FN,Hfdd,F<5»2>  * YIELD CMXM) ♦ HOB CMXM> *PA( MX  Ml 
COMMON/INSTLN/  M . VALUE ♦ VN *K • X X 
COMMON/WPNS/  M»YIELO*CEP#HQB#PA 
COMMON/PARAMS/  UR  ,S I GMA *BETA ♦ I TERS 

QUA  ORATURE  BASE  POINTS  Z<2»  AND  COEFFICIENTS  UCI> 

DATA  Z( i )  »Z(2 ) ,2  <3 ) ,Z<4  )  ,Z<3  )  *U(1 > ,U<2  )»U<3  >  »U<4  )  ,U<5  )/ 

*3,1  43874339  ♦0.4  333933941  »C •  6 794  ;  9566 3  ,i  .6b5i  633  667, 

*0  .9  7  39  j  6526  5 , 0.2955242247 tC  .  2692667193 .2193  863625 , 

.1  494313 49 2,0.  '  666  713433/ 

BE3B(V>  =  1 .: +V* I Bl» V ' < 02 ♦V* < B3*V» ( 04 ♦ V* ( B5* V* 06 >  > > )) 
9E3C<V>=C0  +  V«  <Ci*V*  <C2-V*<C3-Vv  <C4-V*  <C3-VMC6-V* 

*<CT-V«C9>  HID) 

S=3QRT<<X<I )-XX<d))"*2+<X<I+M)-XX<J+N))**2) 

IF<$.LT.0.rci)  sal.' 

CEPAaSQRT<CEP<I  )  **2  ♦*  .231*9  95**  2  > 

IF<WR<lfd>,EQ."  »  THEN 
POMs^.O 
GO  TO  12C 
END  IF 
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THE  PO(I »J>  I S  THE 


CEP-EXCi.tJOEO  CEDF  MODEL. 

DISTANCE  DAMAGE  FUNCTION. 

IFC  CEPA.LT.i-.001  )  THEN 
IFCS.GT.l  >  THEN 

PDMsPOR  CREALC UR (I ♦ J II « S i BET  A C J> ) 

ELSE 
PQH=l.Q 
ENOIF 

ELSE 

CEP -INCLUDED  CEOF  MODEL.  DETERMINATION  OF  THE 
INTEGRATION  LIMITS  A  AND  8. 

NURsl«lT7**REALCUR<I«J) i/cepa 
NS=1.1774*S/C£PA 
d  -1 .  fi«NUR*EXP(2«86*SIaMA(  J) ) 

1FCB.3T.C NS+4.0 >>  B=NS+4.0 
A=N$-4.G 

TFCA.LT*:  .0  >  A  =  0. 

PNsO.O 

GAUS3-LEGENDRE  TEN  POINT  QUADRATURE  TO 
DETERMINE  PDCI,J). 

DO  110  K2=l*5 
DO  100  L2=l*2 

R=Q  .5M  <B-A>*  Zf  K2>*  <-l  >*♦  L2+A  +  8 ) 
PP=PDPCNWR,R*BETACjn 
H=NS*R 

IFCH.EQ.O.J)  THEN 

F«K2»L2 >=PP*R*EXP<«-R**2>/2.0  ) 

£LSEIFCH.LC«3.75)  THEN 
J Js (H/3 .7 5  >  » >2 

F  (K2  »  L2  >  =  PP «  R  *EXP  <-(  NS* 2  >/ 2.:  >*BESB(JJ> 

ELSE 

JJ=%75/H 

FtK2*L2  >-pP*R*EXPC-<NS-R>**2/2..  ) - BESCC JJ ) /5QRT  C H  ) 
END  IF 

FN“FN+U<K2)*F(K2*L2  > 

-  ..  CONTINUE 

110  CONTINUE 

POM*0 .5*<B-A»*FN 
END  IF 

120  P DA  M* INT ( 1C  C J . i * P DM-*3 . « ) /I : 

END 
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RCA L  FUNCTION  PDR ( WR *R »BETA l 

1.  CALCULATES  THE  DISTANCE  DAMAGE  FUNCTION--  THE 
PROBABILITY  OF  ACHIEVING  A  SPECIFIED  LEVEL  OF 
DAMAGE  TO  INSTALLATION  J  FROM  WEAPON  I  WHEN  THEY 
ARE  SEPARATED  A  KNOWN  01  STANCE  R. 


P A3  AMEYER  <£l  =  C*r  7u  5230  7tJ4*E 2=0..:  422320123  *E3s?  .  0  3  9270  52  72* 

*E 4*0. 0)0 152  143  ,E5=  j.OC 0  2765672  *  E6=0  .0 0  0  0 43? 633 > 

REAL  ERF *22  iBET  A*UR#R*U*SIGN*V 

ERF  (Visl«0-l»C/<l»0  ♦V* IE1»V*  CE2*V*<E3*V,J<£4+V*<E5«,V*E6>  »>))**:  6 

CALCULATES  Z2 »  A  STANDARD  NORMAL  PANOOM  VARIABLE*  AND 
TESTS  27  TO  DETERMINE  THE  EXTREMES  OF  THE  PROBABILITY 
OF  DAMAGE  FUNCTION. 

ZZ- 1 1.  /BETA)*L0G<<WR*EXPC-HETA*f*2  )  )/R) 

I F(ZZ .GT.3.67 *  THEN 
PORsl.O 
RETURN 
ENDIF 

XF<ABS(Z2 ).LT.3E-06»  THEN 
PDRrC  .5 
RETURN 
ENDIF 

IF(ZZ.LT.-3.aT>  THEN 
PORs?  .0 
RETURN 
ENDIF 

U=ABS(ZZ)  /SQRT(2.0  ) 

SlGNsl.G 

IFtZZ.LT.C.  )  S I  GN=-1  •  ? 

CALCULATES  THE  (DISTANCE  DAMAGE  FUNCTION  USING  A 
POLYNOMIAL  APPROXIMATION  OF  THE  ERROR  FUNCTION. 

POR“0.5-»SIGN*0.5*ERF<U) 

RETURN 

ENO 
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C********  If  ******** 


*********  ***********************  ****** 


SUBROUTINE  PWMI  h  (  X,  E  ♦  N  *EF  *E  SCALE  •  HA  XlT  *  W«  NU  ) 

i*  P QUELL* S  METHOD  OF  CONJUGATE  DIRECTIONS 
DETERMINES  THE  MINIMUM  OF  A  FUNCTION  USING 
ONLY  FUNCTION  EVALUATIONS. 


REAL  XCN> «E<N  >«U<NU) 

DDM  AG=T  .1  *£SCAl.E 
SCE  4*0  «05/E SCALE 
J  Js  N«  <N*1  ) 

JJJsJJ+N 
K*N  *1 
NFC C=1 
INDsl 
I  NN  “1 

WRI TE(6**  )•  PUM IN* 

DO  4  I-1*N 
U  tl  >=ESCALE 
DO  4  J*i*N 
U  CK 1  =  0  *  0 

I F< I  *  £Q« J  )  W<  K  >  =  ADS (E (I ) ) 
KsK^l 
ITERC=1 
£$GRAD»2 

CALL  FUNCT<  N  »X*F ) 
FKEEP«2.0*ABSCF» 


START  THE  NEXT 

ITERATION. 

ITONEsl 

DO  200  Isl.N 

URITEC6**  >  • 

X(*<I»*)*  *»X(I> 

CONTINUE 

URITECS**)  • 

EF  =  • *-F 

FPsF 

SUMsO*! 

IXPsJJ 

DO  6  I s 1  * N 

I  XPs  r  XP»1 

uc:  xpjsxci) 

I  OI  RN=N  +  l 

I  LINE=  I 

START  THE  NEXT  ONE  DIMENSIONAL  SEARCH. 

DMA  X=«(ILINE> 

OAC  C^DMAX *S CER 
QMAG=AMINl(ODMAG*I  • i *  DMAX) 

0 MA G=AMAXl<OMAGt23. ' *DACC  > 

0  CM  A  X  s  l  ■)  *  0  *  DM  A  G 

GO  TO  (70  *7  *  71  )  »  I  “ONE 
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7'  DL=  j .3 
0=0  MAG 
fpre/=f 
:s=5 

FA=FPREV 

DA=DL 

3  D0=  D-QL 

0  L=  0 

SELECT  THE  NEXT  SEARCH  DIRECTION  FOR  ITERATION  I. 

5f<  K=I  DIRN 

DO  9  X=ltN 

XCI  )  =  XCI)+DD*H<K> 

9  K*X*1 

CALL  FUNCT <  N  *X  »F  I 
NFC  C=NFCC*i 

30  TO  (19«ll»12«13il4t9S)»  IS 
14  IF(  F«EQfiFA)  THEN 

r  F(  AdS  <D>  .LE.DHAX  >  THEN 
D=D»0 
00  TO  8 
,  EUDIF 

URIT£(6«*>*  MAX  CHANOE  DOES  NOT  ALTER  FUNCTION* 
00  ro  2 C 

ELSEIF(F.LT.FA>  THEN 
F  B*F 
09=0 
ELSE 
FB=FA 
09*0  A 
FAaF 
0  A=0 
END  IF 
C 

00  TO  <33*23>»  ISGRAD 
23  QxOB^OB-QA 

IS*  I 
00  TO  <3 

83  D*0  .5*COA*OB-CFA-F8I/COB-OAI » 

IS*  ♦ 

X FC  <D  A-D  >  *<D*OB  I  *GE  *0  >  GO  TO  8 

25  X  S* I 

IF! ABS(D-DB>.LE.OOMAX»  GO  TO  8 
ZA  D=DB*SI3N(D0MAX  »DB-OA I 

IS=1 

DOW  AX  =  DON AX+DOMAX 

OONAGsODMAG+ODMAO 

I  FC  ODMAQ.GE.l.C  E*6?  »  OOMAG=1  *0  E«-60 

IFC DDMAX.LE«DMAX »  GO  TO  8 

DDMAX  =  OMAX 

GO  TO  3 
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13  IFCF.Gt.FAl  GO  TO  23 

2  H  FC=FB 

00*08 

29  FB=F 

08=  Q 

so  To  3n 

12  IFfF.LE.FB)  60  TO  21 

FA»F 
0A*0 
SO  TO  33 

11  IFCF«S£«FB>  60  TO  1J 

FA*F8 
0A=08 
80  TO  29 
C 

71  01=1.1 

DOM AX=3. 0 
FA*FP 
DA*~1  •? 

F8«  FHOLO 
0B«  },•) 

0=1  .) 
i:  FC»  F 

00=0 

33  A=(OB-OC>MFA-FC> 

B«<  OC-OAI*«FB-FCI 

TFC  fA*B>*  (OA-OCI.GT  .3  .  ’  >  GO  TO  34 
FA  =  FB 
OA*OB 
FB*FC 
08=00 
SO  TO  26 

34  0=3  .3*<A*  IDR+OC)*B*  (OA  +  OC  ) ) / (A  +  B ) 

01*08 

F  [a  FB 

IFKFB.LE.FC)  GO  TO  44 
01  =  00 
f:=fc 

44  SO  TO  <66*3&*33>t  I  TONE 

£5  IT0NE*2 

80  TO  43 
C 

0  CHECK  THE  ONE  DIMENSIONAL  MINIMIZATION  SEARCH 

0  FOR  CONVERGENCE* 

C 

36  IF(  A8$<0**DI  ).LE.DACC)  GO  TO  41 

IFC  ABSCO-OI  t.LE.Ct  •  3*ABS<D>M  PC  TO  41 
43  IF< <OA-OC)« CDC-0)*SE.O.: )  THEN 

F  A=F8 
0A=08 
FB*FC 
D8=0C 
GO  TO  23 
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ELSE 
I  S=2 

IF(<0B-O)*(D-0C)*GE*fl*C>  'GO  TO  8 
I  S»3 
GO  TO  8 
ENOIF 
41  FaF! 

0  =0 1 -DL 

RE* (0C-08 >* CDC-DAI* (OA-DB>/(A*B) 

IFC SE.LE.G.r )  THEN 

WRXTECS  »*  >•  ACCURACY  LIMITED  BY  THE  FUNCTION1* 
RETURN 
ENOIF 

DD»SQRT<RE> 

COMPLETES  ONE  OF  THE  N  ONE  DIMENSIONAL  SEARCHES 
FOR  ITERATION  I.  UPDATE  X(l). 

DO  49  IsltN 

Xfl  )=X( I ) ♦D*M<I  DI  RN ) 

W  (  I  OIRN  »  =  OD*U<  I  01  RN) 

9  IDIRNaIOIRN*l 

M(ILINE)aU<ILINE)/OD 
X LI  NEalLINE«l 

I F(  ITONE* EQ *2  I  GO  TO  38 
IF<  CFPREW-F-SUM)*LT  «U  •£  )  GO  TO  94 
SUN«FPREV-F 
JtLalLlNE 

*  IPC XOXRN«LE« JJ)  80  TO  7 

ALL  ONE  DIMENSIONAL  SEARCHES  COMPLETED 

GO  TO  (92*72) «  1N0 
2  FHOLDaF 

ISs6 

IXPaJJ 

00  39  I»l*N 
I XPal XP«-1 

9  w(txp>axm-uaxp> 

0  D»  1.0 

CALCULATES  THE  EXPANDED . POI NT • 

GO  TO  33 
96  GO  TO  (1X2*87) *  IND 
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THE  MODIFICATION  TEST 


C 
C 

c 

112  IFCFP.LE.F)  GO  TO  37 

0-2  .0*<FP»F-2.0>FHOLO)/<FP-F>*  •£ 

IFC  <O*{FP-FHOLD-SUM)<m>2-SUM).GE*0,:  )  GO  T0  37 
8  7  d=JIL*N*l 

IF(J.GT.Jd)  GO  TO  61 
00  62  IsJ*JJ 
K*I-M 

62  U(K>*U(X) 

00  97  IsJILfN 
97  U(I-l)s«(t) 

C 

C  SEARCH  IN  DIRECTION  OF  THE  EXPANDED  POINT. 

C 

61  I DX  RN*I DI RN-N 

IT0NE*3 
K-I 01 RN 

I  XP  *,ld 
AAA  *0  •  j 
DO  67  :*1,M 
I XPaIKP+1 
U(K)sUdXP) 

!F(AAA.LT.ABS(U(K)/E(I  )))  THEN 
AAAsABS (U<K  )/E<I  )  > 

ENOIF 

67  K  =K  ♦! 

0  DM  AGs  l.*i 
U(N)aESCALE/AAA 
ILlNEsN 
GO  TO  7 
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» r*+r, jk  :ji  *.*."«*  .*•*.  .*  .v  * 


i“w»  r  a  ;  *t  r«A 


SOI 


AMB mu 


.  ■  J  •  •  A  i  m/i,  .  t  ^  i  J 


n  n  n 


UPDATE  X  C 1 1  AMD  USD  THE  PREVIOUS  SEARCH  DIRECTIONS. 

37  IXP=JJ 
AAA  sO 
F*FHOLD 
□0  99  1=1 *  N 
IXP=IXP*1 
X<I>=XCI)-U<IXPI 

I F C (AAA*ABS  CE ( I > > ) .LT . A8S ( UCI XP )) )  THEN 
AAA=ABS<U<XXPI/ECI>  > 

ENQIF 

99  CONTINUE 

GO  TO  72 

33  AAA=AAA*<l.r*DI > 

GO  TO  <72*106),  INO 
72  GO  TO  <10  9*88).  INO 

10  9  IF<  AAA.LE*?  <*1  )  GO  TO  20 
C 

1FCF.LT.FP)  SO  TO  35 

URI  TE  <  6  .* ) *  ACCURACY  LIMITED  BY  THE  FUNCTION* 

go  to  2: 

3R  IN0*1 

35  0  DM  AOs*1  •♦*SQRT(A8S<FP-F|) 

IFCODMAQ.se. 1.0E*60 >  0QMAG=1 .CE+60 
I SQR A0=1 
C 

ITERC*ITERC»l 

IFIITERC.LE.MAXIT)  GO  TO  5 
IFCF.LE.FKEEP1  GO  to  20 
F  *FKEEP 
DO  111  1  =  1. N 
4«JJ*  JJJ+1 
111  X(I)  =  UCJJJ) 

21  WRtTECS.O*  ITERATI  0NS  =  *.ITERC 
EF=F 
RETURN 

114  IFC  AAA.LE.:.  .1  S  GO  TO  21 

I NN  =1 
SO  TO  35 
C 

END 
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Appendix  E:  User  Guidelines  and  a  Sample  Problem 

This  appendix  provides  basic  instructions  for  inputting  user- 
specified  weapon  and  installation  parameters  to  the  CEDF  maximization 
algorithm.  These  instructions  are  presented  using  an  example. 

Initially,  the  user  must  determine  the  values  of  the  convergence 
control  parameters,  ACC,  DFFRED,  E,  and  ESCALE.  These  values  depend 
on  installation  values  end  the  number  of  weapons  and  installations. 

User  guidelines  in  Chapter  V  discuss  specific  considerations.  The 
source  code  Initializes  the  PWMLN  convergence  parameters  to  E(i)  "  0.1 
end  ESCALE  “  5000.  Two  minor  code  changes  would  be  necessary  to  change 
either  of  these  values.  The  user  must  input  values  for  the  ZXGGR  con¬ 
vergence  parameters,  ACC  and  DFPRED. 

Next,  the  user  must  decide  on  the  type  of  initial  DGZ  coordinates 
to  use.  The  user  has  three  options.  The  user  communicates  the  desired 
DGZ  location  option  to  the  algorithm  through  the  input  variable  —  INC. 

If  INC  -  1,  then  the  algorithm  usas  tha  user-specified  estimates  of 
the  initial  DGZ  locations.  If  INC  ■  2,  then  the  algorithm  assigns  the. 
coordinates  of  the  m  highest  valued  installations  to  be  tha  initial 
coordinates  of  the  in  weapons  in  decreasing  order:  of  yield.  Finally,  if 
INC  -  3,  then  the  algorithm  assigns  the  coordinates  of  the  m  hardest 
installations  to  be  the  initial  coordinates  of  the  m  weapons  in  de¬ 
creasing  order  of  yield. 

Then  the  user  needs  to  compile  the  necessary  input  data  in  a  FORTRAN 
external  file,  INDATA.  Figure  E~1  is  the  input  data  file,  INDATA,  for  a 
two  weapon- four  installation  example,  The  first  line  is  the  decision 

variable,  INC.  For  this  example,  INC  «  2  and  the  initial  DGZ 
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Fig'.ire  E-l.  The  CEDF  maximization  algorithm  input  file,  INDATA,  for 
a  two  weapon- four  installation  complex. 


coordinates  were  the  coordinates  of  the  two  highest  valued  installations. 
This  second  lj.no  in  the  file  ’NTiATA  is  in,  the  number  of  weapons  for  the 
^complex;  for  this  example,  m  -  2.  The  next  m  lines  contain  the  user- 

i 

specified  weapon  parameters .  The  FORTRAN  input  format  for  these  para¬ 
meters  is  staitamcuit  100  in  subroutine  INITLZ  of  the  source  code.  The 
ordpr  and  units  of  the  weapon  rvirameters  ares  the  longitude  and  lati- 

i  ’ 

tude  coordinates  (these  mint  be  initialized  to  some  value  even  if 
'INC  *  2  or  3),  the.  yield  in  kilotons,  the  CEP  in  feet,  the  HOB  in  feet, 
and  the  Fa. 

The  line  after  the  last  weapon's  parameters  is  n,  the  number  of 
installations  in  the  complex;  for  this  example,  n  -  4.  The  next  n 
lines  contain  the  user-specified  installation  parameters.  Has  FORTRAN 
Input  format  statement  for  these  parameters  is  statement  110  in  sub¬ 
routine  INITLZ  of.  the  source  code.  The  order  of  the  installation 
parameters  is:  the  longitude  and  latitude  coordinates,  a  VNTK  code, 
and  a  value  (a  real  number  less  than  loOOOO.O). 

Finally,  the  ZXGGR  convergence  control  parameters  complete  the 
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external  file,  INDATA.  The  line  after  the  last  installation's  para¬ 
meters  contains  the  value  of  ACC  for  the  second  stage  of  ZXCGR  maxi- 
miation.  The  last  line  of  the  file  contains  the  value  of  DFPRED.  For 
this  example,  ACC  -  0.001  and  DFPRED  ■  1000. 

The  CEDF  maximization  algorithm  outputs  the  results  of  a  problem 
to  another  external  file,  TAPE6.  The  results  of  four  maximizations 
are:  (1)  a  ZX0GR  maximization  for  ACC  -  0.01;  (2)  a  ZXCGR  maximi¬ 
zation  for  a  user-specified  value  of  ACC;  (3)  a  PWMIN  maximization; 
and  (4)  a  mixed  maximization. 

The  next  six  pages  list  the  results  of  the  two  weapon-four 
installation  problem. 
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.#*#■»#!**  4)  B  ••♦****•%#  U  <t  »<*•  1»  *  .**.*■*.*.•<>*<**».*.  4<  *  * 

CEUF  MAXI  MI  ZA  TI-QN  ALGORITHM 

«*»***  A*..**.*.*..****. ...A***  ■**»**.  «***.**♦*..*  »  * 


INITL7 


THIS  PROBLEM  USES  2  WEAPONS 


UEAPON  YIELD  CEP 

HOB 

PA 

1 

70.  2.53. 

IC  C  . 

•  99 

2 

70.  250. 

10  0  0  • 

.99 

IS  COMPLEX 

CONTAINS  ♦ 

INSTALLATIONS 

longitude 

LATITUDE 

VNTK 

VALUE 

*310  ■’  E 

46  33u  N 

16P2 

3350. 

♦  310  30  E 

*6  *  ON 

22P2 

250  .1. 

♦  511  OE 

*6  33i.'N 

21  P* 

500  ■:• 

♦31130E 

*6  3  5N 

1903 

70  J  1  • 

THE  X Y  COORDINATES  GF  THE  INSTALLATIONS  IN  FEET 

XXCIla  *2199  .35367441  XX(5»a  21280.0  0  On  C  Ct  2 
XX (2  >3  4*298.8213981 5  XX<6)s  2*319.99999998 
XX (3  >£  *8*08 .2090*18  *  XX«7)=  2330  6.66666666 

XX<*»«  *8317.75672959  XXC8IS  187*6.6666666* 

INITIAL  DGZ  LCCATIONS  ARE  HIGHEST  VALUED  INSTALLATIONS 

X(lli  *8317.75672559  X«3>  =  187*6.6666666* 

XC2»«  *6*09.2890*18*  '  XC*)s  23306.66666666 

UP  AOS 

UR  f  1  «l  >a  3272 
UR  <2*1  »3  3272 
UR  <1  »2>s  19*7 
UR  (2  «2  I*  19*7 
UR  (1  *3  )  =  2200 
UR (2  t3 )=  2200 
UR  (1  •*  >a  2863 
URI2t*>=  2863 


150 


ft 


ZXC6R  MAXIMIZATION 

t>  *  'I  .1*4. ft.*#..**.*.*. 6. *.< 


ACCs  .U 

WEAPON  1  XY  coordinates: 
WEAPON  2  XY  COORDINATES! 
P$  <1  )s  ,97624 
PS  (2  )*  <,86833 
PS <3  I*  ,009999999999993 
PS <4  )*  *00  9641599994993 
ECs  12293*4438 
WEAPON  1  XY  coordinates: 
WEAPON  2  XY  COORDINATES! 
PS<1)  =  .97228 
PS <2 )s  .62373 
PS  <3  )s  *30  999999999999  8 
PS  <4  )  =  *309851499999993 
EC»  12406*2345 
WEAPON  1  XY  coordinates: 
WEAPON  2  XY  COORDINATES: 
PSClls  .91585 
PS(2>s  *13  315 
PS  (3  )  =  .00  3999099999998 
PSC4)s  *{J0  9940  59999999  8 
ECs  14 1 92*3  653 
WEAPON  1  XY  coordinates: 
WEAPON  2  XY  COORDINATES: 
PSU)s  .65345 
PS  <2  >s  .009999999999999 
PS <3  is  . 64500 
PS  <4)s  .00  9990  C  99999998 
ECs  12125.0943 
WEAPON  1  XY  coordinates: 
WEAPON  2  XY  COORDINATES: 
PS<1  )3  .86239 
PS  <2  )=  .22  :  09 
PS (3  )s  .01196 
PS  <4  )s  *0 199 6r  39999999  3 
ECs  14799.7872 
WEAPON  l  XY  coordinates: 
WEAPON  2  XY  COORDINATES: 
P$U)s  .83566 
PS  (2  )*  .01091 
PS  <3  )s  .0238  60  COOC 0  031 
PS  (4  )s  .0  0  99  7?  29999999  3 
ECs  14358. 6229 
UEAPON  1  XV  coordinates: 
WEAPON  2  XY  COORDINATES: 
PS<1>=  .17632 
PS  <2  )=  .05257 
PS  <3  )=  .63  02  3 
PS  (4  )s  .00937?  299595993 
ECs  13  750.5129 


(  48518.* 
(  46  4  u  8  •  * 


C  43518*  * 
(  46333.* 


C  48517.1 
<  45657.* 


(  48517.* 
<  44168.. 


C  48517.* 
<  4523 1*  * 


C  43517., 
(  45124.* 


<  43516.* 
(  44110., 


16747.  ) 

23307.  I 


107  47.  ) 
23331.) 


16747.  ) 
23554.  ) 


18747.  ) 
24037.  ) 


18747.  ) 
23678.  ) 


18747.  ) 
23729.  ) 


1 67 4  7.  ) 
230  12.  ) 
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PS<2>-  ,0  2  63  3 
PSC3)  =  .13226 
PSC4)  =  ,03  9960  3999999?  8 


EC~  13253.7322 
WEAPON  1  XY  COORDINATES:  ( 
WEAPON  2  XY  COORDINATES:  < 
PSC1)=  .51787 
PSC2>  =  .38321 
PSC31=  .09613 

3  PS  C4  >=  .0?994v599999998 

|  EC*  15416.6958 


WEAPON 

i  xy  coordinates: 

C 

WEAPON 

2  XY  coordinates: 

c 

PSC1>  = 

.51391 

PS  (2  )  = 

.07831 

PS  C3  )* 

.13306 

PS  <4  )  = 

.00  9940  599999998 

EC*  13420.6358 

WEAPON 

i  xy  coordinates: 

c 

WEAPON 

2  XY  coordinates: 

c 

PS  Cl  )= 

.53767 

PSC2  )s 

.06742000000001 

PS  <3  )* 

•o 5319: noooom 

PS  <4  )s 

.0  3994C  599999  992 

EC*  13429.3708 

WEAPON 

i  xy  coordinates: 

c 

WEAPON 

2  XY  coordinates: 

t 

PS  Cl  )= 

.72773 

PS  (2  )* 

•  0228  7 

PS  C3  )«= 

.02584 

PS  <4  )* 

•  0  0  99  3G  49999999  8 

EC*  15196.8465 

WEAPON 

l  xy  coordinates: 

c 

WEAPON 

2  XY  coordinates: 

c 

PS  <1  )* 

.35351 

PS  C2  )= 

.06143 

PS  C3)* 

.0812300900-001 

PS  C4  )* 

.QC  9940599999993 

EC*  15433.0338 

4351 7. ♦  16747. » 
44629.*  23319.) 


4S5 1 7. «  1B747.) 
44312.*  23334.) 


48516.*  18747.) 
448  35. «  23C69.) 


48516.*  18747.) 
45040.*  23336.) 


48516.*  1H747.  ) 
44851. »  23394.  ) 


iep.3  a 

FUNCTION*  13433. 03C6 
FUNCTION  EVALUATIONS:  IB 
X  Cl  )*  43516. 460374C  7 
GR  AO (1 )  *  .000338669791 4471 
XC2>=  44350  .73B  76378 
GRA0C2)*  .01948  70170  614 
XC3I*  1874T. 19690065 
GR  AO  (3  )  *  -.00f  14745020  3560  6 


X  C4  )*  230  93.65196727 
<Jfl  AD  C4  )  *  .1193911576447 


WEAPON 

1 


LONGITUDE 

451129E 

451C37E 


LATITUDE 
46  5  3N 
46  34 7N 


I 
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Fimi^OTfi'mi'imnTrrrrTi'foT  &  mw  r  k’a  tvA  '•  a> LV  •*  iv  ^  ^  l>  ap. 


ZXCGR  WITH  ACC  REDUCED 


ACC  =  ,3  1 

WEAPON  1  XY  COORDINATES:  C  48516,*  18747,1 

WEAPON  2  XY  COORDINATES:  (  44851,*  23394,) 

PS  Ills  ,53351 

PS<2)  =  ,06143 

PS  <3  1=  ,0312900  0000061 
PS  <4  1=  ,0  099  4C  599999998 
EC^  15433,0338 
IE*  =  0 

FUNCTIONS  15433. 03G8 
FUNCTION  EVALUATIONS:  2 
XC1  IS  43516.46037407 
GRAD  Cl.  >*  .0003386697914471 
X (2  Is  44353 .73876378 
GR  AD  (2  I  =  .’194870170614 
X 13  Is  IS  747  ,19690  f*  65 
GRAD  <3  I s  -,00014743023  3560  6 
X  <4  Is  23093.65196727 
GRAD  (4  Is  ,0193°1157644  7 

WEAPON  LONGITUDE  LATITUDE 

1  451129E  46  3  5N 

2  451037E  46  347N 
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***#*.#*  *.*!*«*■*.*****.**  **■**»******■•»■*•*«***•*»«**** 


PUMIN  MAXIMIZATION 


PUMIN 

XU):  43517.75672559 
X  <2  >  =  464 08.2890 4184 
X<3)=  13746.66666664 
X  <4 )=  2333 6. 6 66 66 666 
EF-  12353.3581 
X(l)=  43517.75672559 
X<2  )=  44888 .9  44979  61 
X  C3 Is  18746.66666664 
X  <4  )s  23169.22996394 
EF=  15137.2158 
XXI )s  48517.75672559 
X<2>s  44838.630  71814 
X (3 )=  13746.66666664 
X(4)s  23165. 79903C85 
£F -  15141.9708 

ACCURACY  LIMITED  BY  THE  FUNCTION 
ITERATIONS*  3 
FUNCTION*  15141.9706 

XXHCDs  43517.75672559 
XXXC2>=  44391.57331154 
XXX  < 3  )  =  13746.66666664 
X  XX (4 1*  2 3159 .53991 0  93 

WEAPON  LON 61 TUDE  LA7I TUDE 

1  451 130 E  46  3  5N 

2  45.1C36E  46  348 N 
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*.  V  wV  .%•  UW  L>-  u'W  w  u  Ul 


i  njk .  jfi  ;\4,  .  „*//<,  » V  Wj  •*.  ^  4  4  «£wtV  »•  VN^******.  4  • 


MIXED  TECHNIQUE  MAXIMIZATION 

*♦♦**.*. *».*»****.*.**.*  A************************* 

PWMIN 

X<1 )=  43516.46037407 
XC2  )=  448  30.73876376 
XC3 1=  18747. 19690065 
X<4)=  23093.65196727 
EF«  15133.0608 
XClIs  48316.46037407 
X<2ls  44851.32790633 
X<3>  =  18747. 19690063 
X<4>=  23127.72924026 
EF  s  15136.0308 
XU>  =  48516.46037407 
X<2I=  44851.74210346 
X<3U  18747. 19690065 
X<4J  =  23123*46811256 
EF  s  15138. U108 
X(1  1  =  48516.46037407 
X (2  >“  44853.60958023 
XC3>=  18747.19690  :65 
X<4)=  23123.46811256 
EF  s  13139.4958 

ACCURACY  LIMITED  BY  THE  FUNCTION 

ITERATIONS::  4 
FUNCTION*  13132. 07CG 

X4U  is  48516.460  3740  7 
X4(2I=  4*835:, 35435102 
X4 <3  Is  16747.19690065 
X4  <4  I s  23105.5184597 

NEAP ON 
1 
2 


LONGITUOE  LATITUDE 
451129E  46  3  5N 

451C37E  46  348N 
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Appendix  F:  Verification  of  the  Gradient 
of  the  CEP-Excluded  CEDF  Model' 

The  results  of  two  example  problems  verified  that  the  subroutine 
GFUNCT  correctly  calculates  the  gradient  of  the  CEP-Excluded  model.  The 
pencil  and  paper  results  for  each  example  were  compared  with  the  results 
from  GFUNCT. 

The  first  example  included  one  weapon  and  two  installations.  The 
weapon  and  installation  parameters  are  presented  below.  This  verifi¬ 
cation  example  used  a  graph  of  CEDF(x)  versus  x.  For  40  equally  spaced 
DGZ  locations,  values  of  CEDF(x)  were  calculated.  The  x  direction  was 
along  the  line  connecting  the  two  installations.  Table  F-l  lists  the  40 
values  of  x  and  the  corresponding  function  and  gradient  values.  Figure  9 
in  Chapter  IV  is  a  plot  of  this  data.  A  DGZ  between  the  two  installa¬ 
tions  was  selected  (x  ■  63500),  and  the  gradient  was  calculated  using 
two  methods.  In  this  example,  the  gradient  had  only  one  element  because 
the  y  variable  was  constant;  only  the  x  variable  was  allowed  to  vary. 

The  gradient  values  for  the  two  calculation  methods  were  compared  with 
the  gradient  value  from  GFUNCT. 

A  CEDF 

Tl.~  first  method  used  a  difference  equation,  --  ,  to  approximate 
the  gradient.  The  slope  of  the  line  segment  connecting  the  CEDF  values 
for  the  two  DGZa  on  either  side  of  the  selected  DGZ  was  an  approximate 
gradient  value.  From  Table  F-l, 

CEDF(x  -  63000)  -  6257 
CEDF(x  -  64000)  -  10196 
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TABLE  F-l 


The  CEDF(x)  and  the  Gradient  of  the  CEDF(x)  for 
a  One  Weapon-Two  Installation  Complex 


************ ********* 

GRADIENT  VERIFICATION 


Xfl  > 

Xf  21 

EC 

GRAD(l) 

BR ftD(2  I 

5400  0 • 

20  3  00  « 

2336. 

1.6274 

u .1)900 

54500. 

2  30  3'  . 

321  A. 

1.6494 

0.0001 

5500(5. 

2100: , 

3935. 

1.3530 

0 .0000 

5550  0. 

238  30  . 

4534. 

•  3335 

I.OOCfi 

96003. 

20000  • 

4326. 

.3564 

3.03CG 

5650  C  # 

21 0  O'.  • 

4930. 

.0  895 

3.0  GO  9 

57000. 

20  3  3  r.  . 

495:. 

•  J  1 0  6 

0.0GG C 

57500. 

2  v  3  o ::  • 

495  :  • 

•  -j  0  G  4 

C  .  0  r  2  G 

54000. 

20  3  30  • 

4951. 

•  0  000 

0  •  0  j  O  v 

39  50  0. 

20  3  3':  • 

495*. 

•  0  C 1 2 

.000  3 

390GQ. 

213 OC  * 

4951. 

•  0  005 

o.o: o : 

3a50  0. 

2:3  33  « 

495  0. 

.0017 

0  .00  9  0 

6000 P. 

210  3',  . 

495? . 

.1053 

0.0091 

6133  0. 

20  0  39  • 

4962. 

.1168 

8  .003  0 

61CJ0. 

20  0  0'.  • 

4974. 

•  >513 

3.0090 

61390. 

2 1 0  0  C  • 

5021. 

-1493 

0.0000 

62033. 

2130:  • 

5152. 

•  4064 

V  .0  2  0  0 

6250  C . 

2'j0  0i.;  • 

5485. 

1.0093 

a .0300 

63000. 

203  00  • 

6257. 

2. 20  36 

: .  o  o  o  e 

6350C. 

2100C  . 

7793. 

3.9791 

0.0013 

64030. 

2  u  0  0  C  • 

10196. 

5.4652 

1  .0100 

64530. 

2  ,io  at . 

12910. 

4.5543 

1  •  0  C  C  9 

63000. 

23  3  C-C. 

14724. 

1.9935 

■:  .ooc3 

65330. 

21009  • 

14931. 

-.6246 

1 . 0  ■  > •;  1 

660 3  3. 

2  13  0  9. 

14266. 

-1.5783 

:.oc:  i 

66530. 

2 GO')'.-.  • 

13514. 

-1.3643 

1 .  O  0  trf  >j 

67000. 

29  99C  . 

12921 • 

-1 .  :  0  76 

0.00  CO 

6750 C. 

2CGOO . 

12499. 

-.6735 

7.0.  ■?  •: 

6* 030* 

23309  • 
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The  difference  equation  approximation  of  the  gradient  at  x  -  63500  was 

ACEDF  „  10196  -  6257  .  3  939 
Ax  1000 

The  second  method  was  pencil  and  paper  calculations  of  all  the 
steps  necessary  to  determine  the  gradient.  Chapter  J.I  presented  these 
steps.  Only  a  sunnary  of  the  calculations  are  presented  here. 


Given:  Weapon 

Yield 

CEP  HOB 

Pa 

(x,y)  in  feet 

1 

100-kt 

0  feet  1000  feet 

0.99 

(63500,20000) 

Target 

VNTK 

Value(v) 

WR 

(xx,yy)  in  feet 

1 

11P2 

5000 

6194 

(60000,20000) 

2 

15P2 

12000 

4066 

(68000,20000) 

ftCEDF 

Note:  In  this  example,  y  ■  x(2)  -  constant.  Hence,  ^  ■  ■  ■  0  and 
CEDF(g)  *  CEDF(x). 

Fran  Eq  (1), 

CEDF(x)  -  v1*P*1*Fd(l,l)  +  v2*Pa  *Pd(l,2) 

However,  and  Pa^  are  constants,  so 

CEDF(x)  -  4950*Pd(l,l)  +  11880*Pd(l,2) 


and 


where 


d  CEDF(x) 
dx 


4950* 


+  11880* 


d  Mill) 

dx 


-  x) 


and 


u 


1 

jr 


*  abs 
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For  overpressure  (P-type)  targets,  “  0.2  and 


P-  J-lnd-O 

(F-l) 

-  0.202045 

Lat 

AA  -  xxj  -  Xj 

(F-2) 

88  ■  “jtn  '  Xl-Hn 

(F— 3) 

r  -  M2  +  BB2 

(F-4) 

For  this  verification  example,  BB  "  0  for  both  targets  and  r  -  |AA|, 
the  absolute  difference  between  the  x  coordinates  of  the  weapon  and 
installation  j. 

For  target  1:  AA  ■  60000  -  63500  ■  -3500 

r  -  3500 

u  -  1.8548527 

d  -  -1.8080608*10“5 

dx 

For  target  2:  AA  -  68000  -  63500  -  4500 

r  -  4500 
u  -  0,49780381 

»  3,4247393*10“4 
dx 

Therefore, 

»■  ■— (*>  -  4950 . 0*( -1 , 8080608*10*  )  +  11880.0*(3.4247393*10~4) 
dx 

-  3.97909 
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The  value  of  the  gradient  of  the  CEDF(x)  from  GFUNCT  for  the  DGZ 
selected  (x  ■  63500)  was  3.9791.  The  gradient  results  from  the  dif¬ 
ference  equation  approximation  and  the  pencil  and  paper  calculations 
were  compared  with  the  value  from  GFUNCT.  These  two  comparisons  indi¬ 
cated  the  subprogram  GFUNCT  was  properly  calculating  the  gradient  of 
the  CEDF(x) . 

The  second  gradient  verification  example  included  two  weapons  and 
three  installations.  The  gradient  of  the  CEDF(x)  had  2m  or  four 
elements.  However,  only  one  element  was  completely  checked  by  pencil 
and  paper  calculations. 


Given: 


I 

i 

i 


Weapon 

Yield 

CEP 

1 

100  kt 

0  feet 

2 

100  kt 

0  feet 

Target 

VNTK 

Value(v) 

1 

12P2 

5000 

2 

14P2 

8000 

3 

12P2 

4000 

HOB 

Pa 

(x,y)  in  feet 

1000  feet 

0.99 

(61000,21500) 

1000  feet 

0.99 

(62000,17500) 

WR 

(x,y)  in  feet 

5550 

(60000,19500) 

4495 

(68000,20500) 

5550 

(63000,23500) 

From  Eq  (1) 


«v,*Pa9*  (1  -  Pa,;*Pd(l,l))* 

3x0  11  L  9x9 


+  v *Pa9  (1  -  Pa0*Fd(l,2))* 

*  1  L  3x^ 

+  v  *Pa9  (1  -  Pa0*Pd(l,3)* 

J  1  L  3x2 


(F— 5) 
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From  the  subprogram  FUNCT, 


Pd(l,l)  -  1.000 
Fd(l,2)  -  0.007 
Bd(l,3)  -  0.999 

Hanca,  Vl*pa2  "  4950 

V2*Pa2  "  7920 
V,*Pa2  »  3960 

•nd  Pa2*Pd(l,l)  -  0.99000 

Pa2*H(i,2)  -  0.00693 
Pa2*Fd<l,3)  -  Ov 98901 

Eq  (F-5)  becomes 

-  49.50*  +  7865.1144*^!^^ 

9*2  9X2 

+  43.5204* 

3x2 

(xXj  -  x2) 

(-S2J1)] 

Eqs  (r  *1)  through  (F-4). 


where 


3M2J)  - 


-u 


■Kpr2 


and 


1 

F 


u  ■  tsp  *  abs 


@ ,  AA,  BB,  and  r  are  calculated  from 
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For  target  1:  AA  «  60000  -  62000  -  -2000 

BB  -  19500  -  17500  -  2000 
r  -  2828.427 

-  -3.633135*10”^ 

3x2 

For  target  2:  AA  ■  68000  -  62000  ■  6000 

BB  -  20500  -  17500  -  3000 
r  -  6708.204 

»  2,426656*10"'’ 

0x2 

For  target  3:  AA  -  63000  -  62000  -  1000 

BB  -  23500  -  17500  -  6000 
r  «  6082.763 

iM.2^,)  «  4.304235*10”'’ 

^x2 

Therefore, 

&$9®Sx)  m  49 . 5Q*(  -3 . 633135*10“6 ) 

®X2 

+  7865.1144*(2.426656*10”5) 

+  43.5204*(4.304235*10"5) 

-  0.19254823 

lha  value  of  the  gradient  element  of  the  CEDF(x)  for  x0  from  GFUNCT 
waa  0.19255136.  This  comparison  also  indicated  that  the  routine  GFUNCT 
was  correctly  forming  the  gradient  of  the  CEDF(x). 
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